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This  report  describes  error  models  for 
gravity  gradiometers  in  survey  applications.  The 
major  emphasis  is  on  the  baseline  Gravity  Gradio- 
raeter  Instrument  No.  1  developed  by  the  Bell  Aero¬ 
space  Division  of  Textron,  Inc.,  for  which  compre¬ 
hensive  test  data  were  available.  Preliminary  re¬ 
sults  are  included  for  Bell's  Ball-Bearing  Gradio- 
meter  and  the  Floated  Gradiometer  Instrument  devel¬ 
oped  by  the  Charles  Stark  Draper  Laboratory. 

Models  for  the  vibration  environment  of 
an  airborne  gradiometer  system  are  also  described; 
these  are  based  on  extrapolations  of  available  air¬ 
craft  flight  test  data.  The  report  concludes  with 
examples  of  gradiometric  survey  simulations  based 
on  these  error  models. 
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1. 


INTRODUCTION 


1.1  AIRBORNE  GRADIOMETRIC  GRAVITY  SURVEYS 


The  potential  future  use  of  gravity  gradiomelers  in 
airborne  surveys,  as  presented  in  Fig.  1.1-1,  may  significant¬ 
ly  augment  present-day  survey  techniques.  Two  key  advantages 
to  be  gained  by  using  airborne  gradiometers  are 


•  More  accurate  measurements  of  the  short- 
wavelength  features  in  the  earth's  gravity 
field 

•  Denser  surveys  for  areas  where  conventional 
survey  methods  are  relatively  slow  and  costly. 


R — 51324 
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Gradiometric  data,  when  combined  with  long-wavelength 
data  obtained  from  other  sources  such  as  satellite  altimetry 
and  satellite  tracking,  can  potentially  provide  more  accurate 
estimates  of  gravity  over  larger  areas  than  is  practical  with 
current  survey  methods.  This  report  describes  error  models 
for  gravity  gradiometers ,  which  in  combination  with  the  TASC 

provide  estimates  of 
ted  with  the  inclu¬ 
sion  of  airborne  gradiometry. 

1.2  MAJOR  GRAD IOMETER  ERROR  SOURGES 

Gradient  measurement  errors  are  caused  by 

•  Self  noise  originating  inside  the 

gradiometer 

•  Environmental  disturbances  to  which  the 

gradiometer  is  sensitive 

•  Systematic  errors. 

The  self  noise  depends  only  on  the  gradiometer  and 
sets  a  fundamental  limit  on  the  ultimate  accuracy  of  its  meas¬ 
urements,  independently  of  the  operating  environment.  In  con¬ 
trast,  the  errors  due  to  environmental  disturbances  depend  on 
both  the  gradiometer  and  the  severity  of  the  thermal,  magnetic, 
mechanical,  and  gravitational  disturbances  acting  on  it. 

There  are  standard  methods  for  isolating  the  gradi¬ 
ometer  from  environmental  disturbances  during  moving-base 
operations.  Control  of  the  thermal  environment  requires  the 
use  of  thermal  barriers  and  temperature  control  systems.  Mag¬ 
netic  shields  are  built  into  the  gradiometer  to  attentuate 
both  external  magnetic  fields  and  those  induced  by  current- 
carrying  elements  in  the  gradiometer  itself.  The  gradiometer 


Multisensor  Survey  Simulation  software, 
survey  accuracy  improvements  to  be  expec 


can  be  partially  isolated  from  linear  and  angular  mechanical 
vibrations  by  the  combined  use  of  a  passive  vibration  isola¬ 
tion  system  and  an  inertial  platform. 

The  gradiometer  is  also  naturally  susceptible  to  the 
time-varying  gravity  gradient  fields  caused  by  relative  motions 
of  nearby  masses  such  as  gimbals  in  the  inertial  platform. 
Corrections  can  be  applied  to  the  gradiometer  output  signals 
to  compensate  for  the  gradient  fields  caused  by  disturbing 
masses  whose  positions  can  be  monitored. 


1.3  AVAILABLE  TEST  DATA 

Prototype  gravity  gradiometers  have  been  developed 
for  moving-base  applications  at  both  the  Charles  Stark  Draper 
Laboratory  (Draper)  and  the  Bell  Aerospace  Division  of  Textron, 
Inc.  (Bell).  The  comprehensive  body  of  formal  test  data  on 
Bell's  baseline  Gravity  Gradiometer  Instrument  (GGI  No.  1)  was 
the  primary  data  base  used  for  developing  the  error  models 
described  in  this  report.  These  data  were  supplemented  by 
recently  available  data  on  a  modified  version  of  the  Bell  in¬ 
strument  known  as  the  Ball-Bearing  GGI.  Recent  data  for  the 
Draper  gradiometer  were  also  analyzed.  Error  models  for  the 
Draper  gradiometer  have  not  been  as  fully  developed  as  those 
for  the  Bell  instrument,  but  preliminary  results  are  related 
to  survey  performance. 

Bell  is  currently  building  and  testing  smaller  ver¬ 
sions  of  their  baseline  gradiometer  that  employ  ball  bearings 
for  the  rotating  assembly.  The  design  of  the  ball-bearing 
gradiometers  is  similar  to  the  baseline  GGI,  and  the  error 
models  described  in  this  report  are  therefore  applicable  to 
these  newer  instruments.  Only  preliminary  test  data  on  the 
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ball-bearing  prototype  are  currently  available.  Error  model 
refinements  associated  with  additional  test  data  from  the  new 
instruments  are  likely  to  consist  of  parameter  adjustments. 

The  dominant  instrument  errors  during  an  airborne 
survey  are  expected  to  be  caused  by  the  self  noise  and  the 
mechanical  vibration  environment  of  the  gradiometer.  Ample 
test  data  on  the  Bell  GG1  No.  1  were  available  for  modeling 
these  error  sources  (Refs.  1-11).  Temperature  sensitivity 
data  (Refs.  3,  11,  12)  indicated  that  thermally-induced  errors 
should  be  small  enough  to  ignore  compared  to  other  types  of 
errors  under  survey  conditions.  This  conclusion  is  based  on 
Bell's  measured  thermal  sensitivity  coefficients,  which  are 
less  than,  or  equal  to,  49  eotvos  per  degree  Celsius  change 
in  the  temperature  of  the  gradiometer ' s  instrument  block. 

1.4  ERROR  MODELS  FOR  SURVEY  SIMULATIONS 

The  accuracy  with  which  gravimetric  quantities,  such 
as  anomalies  and  vertical  deflections,  can  be  recovered  from 
gradiometric  survey  data  depends  on  the  characteristics  of  the 
gravity  field,  the  geometry  of  the  survey,  the  data  available 
from  other  sensors,  and  the  errors  in  the  gradiometer  measure¬ 
ments.  The  TASC  Multisensor  Survey  Simulation  software  uses 
models  for  the  gravity  fields,  survey  geometries,  and  avail¬ 
able  sensors  to  estimate  the  rms  errors  inherent  in  a  statis¬ 
tically  optimal  recovery  of  gravity  products  through  multi¬ 
dimensional  Wiener  filtering.  These  calculations  are  carried 
out  in  the  frequency  domain.  Therefore,  the  gradiometer  error 
models  described  in  this  report  are  represented  by  power  spec¬ 
tral  densities.  These  densities  provide  an  easily  understand¬ 
able  characterization  of  the  random  errors  due  to  gradiometer 
self  noise  and  aircraft  vibrations. 


1.5 


OVERVIEW  OF  REPORT 


Models  for  the  self  noise  of  the  Bell  and  Draper  gradi- 
ometers  are  described  in  Chapter  2;  these  models  were  derived 
from  test  data  by  using  two  different  methods  of  spectral  esti¬ 
mation.  Chapter  3  describes  models  for  vibration-induced  errors 
in  the  Bell  baseline  instrument.  Also  in  Chapter  3,  models 
for  the  vibration  environment  of  an  airborne  gradiometer  are 
developed  from  aircraft  flight  test  data.  This  chapter  ends 
with  an  error  model  for  a  complete  survey  system  consisting  of 
three  gradiometer  instruments  mounted  on  one  inertial  plat¬ 
form.  Examples  of  gradiometer  survey  simulation  results  ob¬ 
tained  with  the  survey-system  error  model  are  described  in 


Chapter  4.  Chapter  5  describes  the  status  of  ongoing  gradi 


ometer  development  and  discusses  key  issues  for  future  air¬ 


borne  gradiometer  test  plans.  The  report  ends  with  conclu¬ 


sions  and  recommendations  in  Chapter  6. 
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2. 


GRAD IOMETER  SELF  NOISE 


2.1  INTRODUCTION 

This  chapter  describes  techniques  for  deriving  gradi- 
ometer  self-noise  models  from  test  data.  The  analysis  of  data 
for  the  Bell  baseline  Gravity  Gradiometer  Instrument  (GG1)  is 
described  in  detail.  Recently  obtained  self-noise  spectra  for 
the  Bell  ball-bearing  and  Draper  floated  gradiometers  are  also 
presented . 


A  series  of  experiments  (Refs.  1,  3,  4,  10,  11)  were 
performed  at  Bell  to  measure  the  self  noise  of  their  baseline 
GG1.  These  wei  e  conducted  in  accordance  with  the  test  procedure 
TP1  (Ref.  A).  The  data  sets  discussed  in  this  report  were 
obtained  in  two  experiments:  one  with  the  gradiometer  spin 
axis  vertical  (SAV)  and  the  other  with  the  spin  axis  hori¬ 
zontal  ( SAH ) .  These  experiments  spanned  13  days  for  the  SAH 
orientation  and  almost  11  days  for  SAV.  An  overview  of  these 
data  is  shown  in  Figs.  2.1-1  and  2.1-2.  r.ach  figure  shows 
both  the  inline  and  cross-channel  outputs  of  the  GGI  as  a 
function  of  lime.  The  time  histories  in  Fig.  2.1-2  clearly 
show  start-up  transients,  while  the  plots  in  Fig.  2.1-1  are 
free  of  start-up  transients  because  the  GGI  had  been  running 
several  days  before  these  data  were  taken. 


-This  terminology  reflects  the  fact  that  the  Bell  gradiometer 
simultaneously  measures  "inline"  and  "cross"  elements  of  the 
gravity  gradient  tensor.  A  description  of  these  measurements 
is  given  in  Appendix  B. 
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Figure  2.1-1  Bell  GGI  No.  1  Self-Noise  Data 
(SAV,  April  Data,  Ref.  24) 
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Figure  2.1-2 


Bell  GGI  No.  1  Self-Noise  Data 
(SAH,  March  Data,  Ref.  11) 
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2.2 


PREPROCESSING  OF  DATA 


The  gradiometer  noise  data  were  analyzed  in  three 
stages.  The  first  stage,  termed  preprocessing,  is  described 
in  this  section.  The  laboratory  data  are  sample  values  of  the 
outputs  of  the  cross  and  inline  channels  of  the  gradiometer. 

The  sampling  rate  for  each  channel  was  0.5  samples  per  second, 
and  the  data  were  scaled  to  read  directly  in  eotvos  (E) .  The 
self  noise  consisted  of  small  fluctuations  about  relatively 
large  mean  values  once  the  initial  start-up  transients  died 
out.  These  mean  values  were  arbitrary  offsets  and  were  there¬ 
fore  subtracted  from  the  data  before  the  power  spectral  den¬ 
sities  of  the  noise  were  estimated.  This  removal  of  the  means 
reduced  low-frequency  bias  in  the  spectral  estimates. 

To  facilitate  low-frequency  analysis,  the  data  were 
lowpass  filtered  and  resampled  at  1/64  the  original  sampling 
rate.  This  "decimation"  was  implemented  in  two  steps  with  a 
finite-impulse-response  (FIR)  digital  filter.  The  first  step 
was  to  filter  contiguous  blocks  of  laboratory  data.  Each  of 
these  blocks  contained  8192  sample  pairs  (cross  and  inline 
data),  which  at  the  output  of  the  FIR  filter  produced  1024 
sample  pairs  with  a  sampling  rate  of  0.0625  samples  per  sec¬ 
ond.  The  second  step  was  to  form  contiguous  blocks  from  these 
filtered  samples  and  pass  them  again  through  the  FIR  filter  to 
produce  samples  having  the  final  desired  rate.  In  this  way 
seven  blocks  of  decimated  data,  each  containing  1024  sample 
pairs,  were  generated  for  the  SAH  orientation,  and  five  blocks 
were  obtained  for  SAV.  The  FIR  filter  was  designed  by  using  a 
Kaiser  windowing  technique,  and  the  proper  operation  of  the 
filter  was  verified  through  the  use  of  sinusoidal  test  signals. 

The  two  data  tapes  were  dated  March  9  and  April  25, 
197(8.  Plots  of  these  data  were  previously  shown  in  Figs.  2.1-1 


and  2.1-2.  During  the  preprocessing  of  the  April  25  data,  the 
gap  due  to  tape  recorder  failure  was  avoided  since  it  would 
otherwise  have  introduced  bias  into  the  spectral  estimates. 

The  first  day  of  data  in  the  March  9  tape  was  also  avoided 
because  of  the  start-up  transient. 

2.3  SPECTRAL  ESTIMATION  METHODS  EMPLOYED 

The  analysis  of  the  gradiometer  noise  data  was  con¬ 
ducted  so  as  to  balance  the  conflicting  goals  of  obtaining 
spectral  estimates  that  were  both  precise  and  unbiased.  There 
was  also  the  need  to  verify  that  the  data  were  consistent  with 
a  stationary  random  model  for  the  self  noise.  Without  evidence 
of  stationarity  the  entire  procedure  of  modeling  the  self  noise 
with  a  power  spectrum  could  be  misleading. 

The  preprocessed  data  were  analyzed  in  blocks  con¬ 
taining  1024  sample  pairs  (inline  and  cross-channel  data). 

These  blocks,  in  the  case  of  decimated  data,  each  spanned  about 
1.5  days,  which  was  long  enough  to  provide  spectral  estimates 
down  to  the  lowest  frequencies  of  interest  in  gravity  survey 
simulations.  On  the  other  hand,  the  blocks  were  short  enough 
to  test  for  nonstationary  behavior  over  the  durations  of  the 
data  sets.  The  block  length  also  provided  spectral  estimates 
from  the  raw  and  decimated  data  that  overlapped  each  other  for 
a  frequency  range  of  four  octaves.  This  overlap  provided  a 
consistancy  check  for  the  estimates. 

As  a  further  consistancy  check,  two  methods  of  spec¬ 
tral  estimation  were  used:  smoothed  fast  Fourier  transforms 
(FFT)  of  Kaiser-windowed  data,  and  least -squares  autoregressive 
(AR)  models.  Both  methods  were  used  to  estimate  autospectra, 
c ross- spec t ra ,  and  squared  coherencies  for  the  two  channels  of 
gradiometer  noise. 


These  spectral  quantities  are  defined  in  terms  of  the 
correlation  functions  (R^(t),  ( T  )  *  RXc(*))  t^e  inli-ne 

and  cross-channel  noise  signals,  which  are  respectively  denoted 
fj(t)  and  rc(t) .  The  correlation  functions  are 


Il(x  )  =  E 

[r x ( t+x )  r j ( t )] 

Cc ( x )  -  E 

Tc<tn>  rc(t>] 

XC(T)  =  E 

rI(t+x )  rc(t)] 

(2.3-1) 


where  "E"  denotes  the  expected  value,  t  is  time,  and  t  is  the 
time  shift. 


The  autospectra  of  the  inline  and  cross-channel  noises 
are  the  Fourier  transforms  of  their  respective  correlation 
functions 


Su(f)  -- 

00 

L 

Rn(x)e'l27TfTdx 

scc(f)  = 

00 

/. 

RCC(x)e  l2nftdx 

spectrum  is 

slC(f>  = 

00 

L 

RIC(t )e'l2nfxdt 

(2.3-2) 


(2.3-3) 


where  f  is  the  frequency,  -<*  <  f  <  ®.  These  power  spectra  are 
"two-sided,"  which  means  that  the  mean-square  value  of  each 
noise  signal  is  computed  by  integrating  its  spectrum  over  both 
positive  and  negative  value  of  frequency,  e.g. 
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E 


(2.3-4) 


[r]}--  £sn(f)dt 

The  squared  coherency  between  the  inline  and  cross-channel 
noises  is 


C(f) 


lsIC(  ^1 

SCC^ 


(2.3-5) 


This  function  is  useful  in  analyzing  gradiometer  data  because 
it  indicates  whether  correlations  between  the  noise  signals 
Tj(t)  and  l~c(t)  need  to  be  included  in  gradiometer  error  models 
for  survey  simulations.  The  values  of  C(f)  fall  in  the  range 
from  0  to  1  and  indicate  the  fraction  of  power  that  is  common 
to  fj(t)  and  rc(t)  as  a  function  of  frequency.  If  the  inline 
and  cross-channel  noises  are  uncorrelated,  then  C(f)  =  0, 
while  C(f)  =  1  indicates  that  the  two  noises  are  determinis¬ 
tically  related  to  each  other  by  a  linear  transformation. 

If  C(f)  <<  1,  then  the  cross  and  inline-channel  noises  can 
be  modeled  as  uncorrelated  processes. 


2.3.1  Fast  Fourier  Transform  (FFT) 

The  FFT  method  of  spectral  estimation  consists  in 
multiplying  the  data  time  series  by  a  data  window  and  then 
fast  Fourier  transforming  the  result.  The  Kaiser  window  was 
used;  it  has  the  values  (Ref.  14) 

wk  =  lo  («Vi  ■(~rli))/i0M  (2-3-m 


where 
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k  -  1  ,2 , . .  .  ,n 


W(k)  =  window  value  at  the  k^  sample 

I  =  modified  Bessel  function  of  order  zero 

n  =  number  of  samples  in  the  time  series  block 
being  analyzed  (1024  for  this  study) 

a  =  resolution  parameter  (7  for  this  study). 

The  value  of  the  resolution  parameter  was  chosen  to 
provide  large  side-lobe  attenuation  in  the  frequency  domain 
(72  dB)  while  preserving  adequate  frequency  resolution  with  a 
major-lobe  width  of  nine  frequency  bins  (one  frequency  bin  is 
equal  to  the  sampling  frequency  divided  by  n). 

For  each  block  of  windowed  data  an  unsmoothed  spec¬ 
tral  estimate  was  obtained  by  fast.  Fourier  transforming  the 
block  and  then  computing  the  periodogram,  which  for  windowed 
data  is  the  squared  modulus  of  the  spectrum  divided  by  the  sum 
of  the  squares  of  the  window  values.  This  periodogram  is  an 
unbiased  estimate  of  spectral  density  when  the  original  data 
are  a  sample  of  stationary  white  noise. 

The  unsmoothed  spectral  estimate  has  low  precision 
(large  variance).  Smoothing  was  performed  on  the  decimated 
data  by  computing  unsmoothed  estimates  from  as  many  blocks  as 
were  available  and  then  averaging  them  together  at  each  fre¬ 
quency  in  the  FFT.  Five  of  these  data  blocks  were  available 
for  SAV  and  seven  were  available  for  SAH.  Five  to  seven  block 
averages  were  also  used  to  smooth  the  high-frequency  estimates 
from  raw  data. 

The  resulting  spectral  estimates,  when  plotted  with 
respect  to  logarithmic  scales,  showed  that  additional  smoothing 
was  warranted  at  mid  and  high  frequencies.  This  was  implemented 
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by  performing  a  running  mean  on  the  spectral  estimates.  The 
number  of  frequency  bins  in  this  mean  was  adjusted  so  that  each 
average  was  performed  over  an  interval  of  frequencies  having  a 
width  equal  to  30%  of  its  center  frequency.  This  "constant 
percentage  bandwidth"  running  mean  provided  final  spectral 
estimates  having  a  nearly  constant  frequency  resolution  when 
plotted  against  a  logarithmic  frequency  axis.  The  center  fre¬ 
quencies  were  chosen  to  avoid  overlap  of  adjacent  means. 

The  resulting  estimates  for  autospectra  and  magnitudes 
of  cross-spectra  are  shown  in  Figs.  2.3-1  to  2.3-6,  where  the 
dotted  lines  show  the  sample  standard  errors  (one  sigmas)  of 
the  estimates.  The  highest  frequency  in  these  estimates  is  24 
millihertz  (mHz),  which  is  the  cut-off  frequency  of  the  lowpass 
filter  in  the  Bell  gradiometer.  Above  this  frequency  the  self¬ 
noise  spectra  fall  sharply  due  to  the  filter  cut-off  rate  of  24 
decibels  (dB)  per  octave.  The  lowest  frequency  in  the  estimates 
is  the  reciprical  of  the  block  length  of  the  decimated  data 
(1.5  days ) . 

Although  the  cross-spectral  estimates  are  statistically 
significant  at  most  frequencies,  their  magnitudes  are  small  com¬ 
pared  to  the  corresponding  autospectra  for  the  two  output  chan¬ 
nels  of  the  gradiometer.  This  implies  that  the  inline  and  cross¬ 
channel  noises  can  be  modeled  as  uncorrelated  for  the  purposes 
of  estimating  gravity  gradients  from  gradiometer  measurements. 

2.3.2  Autoregressive  Model  Fitting 

The  second  method  of  spectral  estimation  that  was 
used  in  this  study  is  called  least-squares  autoregressive  (AR) 
model  fitting  (Ref.  15-20).  The  AR  method  has  the  merit  of 
yielding  spectral  estimates  with  greater  accuracy  than  the  FFT 
technique  when  t  he  data  ('an  be  modeled  by  an  au  t  o  regress  i  ve 
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Figure  2.3-5  Self-Noise  Autospectrum  (Cross  Channel , 

SAV,  25  April  Data,  FFT  Estimates) 
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Figure  2.3-6 


Self-Noise  Cross-Spectrum  (Inline  ard  Cross 
Channels,  SAV,  25  April  Data,  FFT  Estimates) 


model  of  sufficiently  low  order.  Other  names  associated  with 
this  approach  to  spectral  estimation  are  "maximum  entropy," 
"linear  prediction,"  and  "all-pole  modeling." 

The  AR  method  used  in  this  study  is  explained  in  the 
following.  Suppose  the  vector  x(k)  is  the  sample  of  a  vec¬ 
tor  time  series  of  data,  where  k  =  l,2,...,n.  (In  analyzing 
gradiometer  self-noise  data  for  this  report,  n  =  1024  and  x(k) 
was  two-dimensional  because  the  cross  and  inline-channel  out¬ 
puts  were  simultaneously  measured.)  Suppose  further  that  x(k) 
is  an  estimate  of  x(k).  The  error  of  this  estimate  is  then 
defined  as  e(k)  =  x(k)  -  x(k).  Therefore 

x ( k )  =  x(k)  +  e(k)  (2.3-5) 

The  AR  method  is  based  on  the  assumption  that  the 
errors  e(k),  k=l,2,...,n,  are  a  sample  of  zero-mean  stationary 
white  noise  when  each  estimate  x(k)  is  chosen  to  be  the  best 
linear  one-step-ahead  prediction. 

x(k)  =  A1x(k-1)  +  A2x(k-2)  +  ...  +  Apx(k-p), 

k  =  p+1,  p+2,...,n  (2.3-6) 

In  this  equation  the  positive  number  p  is  called  the 
order  of  the  predictor,  while  the  square  matrices  A^  ,  A2i-*»A  , 
which  contain  real  numbers,  are  called  parameter  arrays.  Us¬ 
ing  (2.3-6)  in  (2.3-5)  yields  the  following  AR  model  for  x(k) 

x(k)  =  A1x(k-1)  +  A2x(k-2)  +  ...  +  Apx(k-p)  +  ef(p,k), 

k  =  p+1,. . .  ,n  (2.3-7) 

where  the  one-step-ahead  or  forward  prediction  error  is  e^(p,k). 
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The  data  time  series  x(k),  k=l,2,...,n,  is  modeled  as 
a  sample  from  a  stationary  random  process.  Therefore,  it  also 
makes  sense  to  define  a  one-step-back  predictor 

x(k)  =  B^xCk+l)  +  B2x(k+2)  +  ...  +  BpX(k+p) 

k  =  1,2,.. . ,n-p  (2.3-8) 

Using  Eq.  2. 3. -8  in  Eq.  2.3-5  yields  a  second  AR  model 

for  x(k) 

x(k)  =  BjX(k+l)  +  B2x(k-+2)  +  ...  +  BpX(k+p)  +  e^pjk), 

k  =  1,2, . . .n-p  (2.3-9) 

where  the  one-step-back  or  backward  prediction  error  is  e^(p,k). 
The  AR  method  is  based  on  the  assumption  that  this  backward 
prediction  error  sequence  is  also  a  sample  of  zero-mean  sta¬ 
tionary  white  noise  when  the  parameter  arrays  Bp,  B2,...Bp  are 
optimally  chosen. 

By  specifying  the  order  p,  the  parameter  arrays,  and 
the  covariance  matrix  of  the  prediction  errors  for  an  AR  model, 
one  determines  an  estimate  of  the  power  spectrum  of  the  process 
generating  the  time  series.  The  forward  AR  model  yields  the 
estimated  spectral  matrix 

Sp(F)  =  A_1(F)CFA"(F),  |F |  <  1/2 

F  =  f/f  =  frequency/sampling  frequency 

b  a 

JL 

A  =  complex  conjugate  transpose  of  A 
n 

A(F)  =  I  -  £  A, exp ( i 2nkF )  (2.3-10) 

k-p+1  K 

I  =  ident i t y  matrix 
Cp  =  covariance  matrix  of  ep(p,k) 
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The  matrix  f  ^S-ff/f  ),  where  f  =  sampling  frequency  of  the 
data  series,  contains  estimates  of  autospectra  along  its  diagon 
al  and  estimates  of  cross-spectra  in  its  off-diagonal  locations 

A  similar  formula  yields  the  spectral  estimates  pro¬ 
duced  by  backward  predictors 

Sb(F)  =  B_1 (F )CgB*(F ) ,  j  F |  <  1/2  (2. 3. -11) 

Several  methods  for  optimally  fitting  AR  models  to 
vector  time  series  have  been  discussed  in  the  literature.  The 
method  used  in  the  present  study  satisfies  the  following  re¬ 
quirements,  which  are  appropriate  for  the  purposes  of  spectral 
estimation:  For  specified  order  p,  sample  size  n,  and  time 

series  x(k),  the  sum  of  the  squared  forward  and  backward  pre¬ 
diction  errors,  SE(p),  is  minimized  subject  to  the  requirement 
that  the  forward  and  backward  estimates  of  the  spectrum  must 
be  identical.  In  other  words,  the  parameter  matrices  Ai ,  B^, 
i=l,2,...,p,  and  covariance  matrices  Cp  and  Cg  are  chosen  to 
minimize 


n-f 


SE(P)  =  £  |ef(p,k)|  ^  +  2-  |eh(p.k) 

k=p+l  r  1  k=l  1  D 


(2.3-12) 


-f,b^p,kM  =  -f  b^p,k^-f  b^p,k^  =  inner  Pr°duct 


subject  to  the  constraints  that  Eqs.  2. 3. 2-3  to  2. 3. 2-7  are 
satisfied  and  the  two  spectral  estimates  are  equal: 


Sf(F>  =  Sb(F),  J  F |  <  1/2 


(2.3-13) 


The  algorithm  used  for  computing  the  solution  to  this  problem 
recursively  was  presented  in  Ref.  15. 
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The  AR  spectral  estimate  obtained  for  a  particular 
time  series  depends  on  the  order  p  of  the  model.  In  this  study 
the  Akaike  information  criterion  (AIC  statistic)  was  used  to 
select  the  value  of  p  for  each  time  series  analyzed  (Refs.  18 
through  20).  AR  models  were  fitted  for  p  =  1,2,... 20,  and  the 
value  of  AIC  was  computed  in  each  case.  The  value  of  p  that 
yielded  the  minimum  value  of  AIC  was  then  used  to  compute  the 
spectral  estimate.  This  method  of  order  selection  was  employed 
because  it  can  be  derived  from  basic  information- theoretic 
principles  and  has  been  found  to  be  suitable  for  a  wide  vari¬ 
ety  of  estimation  problems  (Refs.  18-20). 

Standard  errors  (standard  deviations)  were  computed 
for  selected  AR  spectra  to  assess  their  precisions.  The  param- 
eter  covariance  routines  in  PARAIDE1  were  used  for  this  pur¬ 
pose.  These  routines  compute  the  inverse  Fisher  information 
matrix  corresponding  to  the  particular  AR  model  under  study. 
This  information  matrix  equals  the  covariance  matrix  of  the 
parameters  in  the  parameter  arrays  of  the  AR  model  when  the 
one-step  prediction  errors  of  the  model  are  Gaussian  and  the 
sample  size  is  large. 

Small  perturbations  in  the  AR  parameter  values  pro¬ 
duce  small  perturbations  in  the  AR  spectra  of  the  gradiometer 
data  because  the  spectra  are  free  of  lightly  damped  resonances 
in  the  frequency  range  of  interest.  These  perturbations  are 
related  to  each  other  through  sensitivity  functions  that  can 
be  computed  from  the  AR  model.  Therefore,  one  can  estimate 
the  standard  deviations  of  the  spectrum  estimates  due  to  the 
sampling  variability  of  the  AR  parameters  by  computing  these 


PARAIDE,  a  trademark  of  The  Analytic  Sciences  Corporation, 
refers  to  a  set  of  computer  routines  for  parameter  identi¬ 
fication  and  stochastic  modeling. 
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sensitivity  functions  together  with  the  covariance  matrix  of 
the  AR  parameters  (inverse  Fisher  information  matrix). 

Gradiometer  self-noise  data  were  analyzed  in  blocks 
containing  1024  pairs  of  measurements  of  the  inline  and  cross¬ 
channel  outputs.  The  resulting  AR  spectral  estimates  were 
then  plotted  at  512  frequencies,  which  correspond  to  the  fre¬ 
quency  bins  of  a  1024  sample  FFT.  Typical  examples  of  the 
results  obtained  in  this  way  for  the  SAH  orientation  are  shown 
in  Figs.  2.3-7  to  2.3-11,  where  the  dotted  lines  are  the  com¬ 
puted  standard  errors  (one  sigmas)  for  the  estimates.  Similar 
results  were  obtained  for  the  SAV  orientation  except  that  the 
low-frequency  rise  in  the  auto  and  cross-spectra  occured  at  a 
lower  frequency. 

A  few  comments  on  the  interpretation  of  AR  spectral 
estimates  are  useful  at  this  point.  In  Figs.  2.3-7  through 
2.3-9,  the  estimates  from  decimated  data  have  ripples  above 


Figure  2.3-7 


Self-Noise  Autospectrum  (Inline  Channel, 
SAH,  9  March  Data,  AR  Ksti mates) 
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Figure  2.3-8  Self-Noise  Autospeetrum  (Cross  Channel 
SAH,  9  March  Data,  AR  Estimates) 


Figure  2.3-9 


Self-Noise  Cross-Spectrum  (Inline  and 
Cross  Channels,  SAH,  9  March  Data, 

AK  Estimates) 


0.5  niHz .  This  is  typical  behavior  for  AR  spectra  derived  from 
data  containing  both  low-frequency  random  drift  and  white  noise. 
The  spectral  estimates  from  raw  data,  in  contrast,  are  free  of 
ripples  because  they  are  derived  from  data  batches  having  short 
durations  (34  minutes).  This  duration  is  appropriate  for  high- 
frequency  analysis  but  is  too  short  to  yield  an  accurate  esti¬ 
mate  of  the  spectrum  at  low  frequencies,  which  characterizes 
the  long-term  random  drift.  Therefore,  the  spectra  from  deci¬ 
mated  and  raw  data  are  expected  to  have  different  qualities  in 
their  frequency  range  of  overlap  from  0.5  mHz  to  3.9  mHz.  The 
differences  between  decimated  and  raw  estimates  are  judged  to 
be  statistically  insignificant  at  a  two-sigma  (95%)  confi¬ 
dence  level . 

Two  conclusions  can  be  drawn  from  these  AR  spectral 
estimates.  The  first  is  that  the  self  noises  in  the  two  output 
channels  are  uncorrelated  for  the  purposes  of  gravity  gradient 
measurements.  This  conclusion  derives  from  the  AR  coherency 
plot  in  Fig.  2.3-11,  which  indicates  less  than  10%  coherency  in 
the  frequency  bands  where  the  estimates  are  statistically 
significant . 

The  second  conclusion  is  that  the  AR  spectral  esti¬ 
mates  based  on  decimated  data  are  statistically  consistent 
with  the  high-frequency  estimates  derived  from  raw  data,  i.e. 
in  the  frequency  range  of  their  overlap  the  t wo- sigma  (95%) 
confidence  intervals  of  the  estimates  intersect  each  other. 

(The  dotted  curves  in  the-  spectral  plots  (Figs.  2.3-7  to  2.3-11) 
art-  only  one-sigma  ((>8%)  conf  idence  intervals).  The  one  excep¬ 
tion  to  this  consistency  occurs  in  the  phase  of  the  cross-spectrum 
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(Fig.  2.3-10).  Fortunately,  precise  estimates  of  this  phase 
function  are  not  needed  for  developing  gradiometer  error  models 
because  the  two  output  channels  contain  essentially  incoherent 
noises . 


Three  sets  of  AR  spectral  estimates,  such  as  illus¬ 
trated  in  Figs.  2.3-7  to  2.3-11,  were  compared  with  each  other 
to  validate  the  use  of  a  stationary  model  for  the  gradiometer ' s 
self  noise.  These  estimates  were  generated  from  non-overlapping 
blocks  of  both  SAH  and  SAV  decimated  and  raw  data.  The  com¬ 
parisons  showed  no  statistically  significant  evidence  of  non¬ 
stationary  behavior  in  the  frequency  range  below  24  mHz,  i.e., 
the  two-sigma  confidence  intervals  of  paired  estimates  over¬ 
lapped  each  other.  At  higher  frequencies,  above  50  mHz,  non¬ 
stationary  behavior  was  observed  in  the  form  of  time-varying 
narrow-band  spectral  peaks  in  both  the  spectrum  and  coherency 
plots  for  SAH  and  SAV.  These  components  were  above  the  in¬ 
tended  24  mHz  cutoff  frequency  of  the  Butterworth  low-pass 
filters  in  the  gradiometer  and  did  not  contribute  significant 
power  to  the  self  noise.  For  survey  applications  requiring 
one  measurement  every  10  seconds,  these  resonances  would  not 
be  a  practical  concern  because  the  required  bandwidth  is  50 
mHz . 

2.4  SELF-NOISE  MODELS 
2.4.1  Spectral  Model 

The  spectral  estimates  for  the  Bell  baseline  GGI  ob¬ 
tained  by  the  FFT  and  AR  methods  support  each  other  and  lead 
to  three  important  conclusions. 
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•  At  lower  frequencies  (  |f|  <  0.1  mHz)  the 
autospectra  are  inversely  proportional 
to  the  frequency  squared  (trend  =  -20 
dB/decade ) 

•  At  high  frequencies  (2  mHz  <  |f|  <  24 

mHz)  the  autospectra  are  independent  of 
frequency 

•  Over  the  range  of  frequencies  relevant 

to  airborne  gradiometric  surveys  0.01 
mHz  <  | f |  <50  mHz,  the  self  noises  in 

the  inline  and  cross  channels  are  inco¬ 
herent  and  have  nearly  the  same  auto- 
spec  t  ra . 


In  view  of  these  results,  an  appropriate  analytic 
noise  model,  which  summarizes  all  of  the  FFT  and  AR  spectral 
estimates  for  the  Bell  baseline  GG1  ,  consists  of  two  uncorre 
lated  processes  with  power  spectra 


SSAV 

SSAH 


(f) 

(f) 


R 


SAV  +  w 
f2  WSAV 


R 


SAH  +  w 
c2  SAH 


0.01  mHz  <  |  f |  <  24  mHz 

(2.4-1) 

0.01  mHz  <  |  f |  <  24  mHz 


These  models  correspond,  respectively,  to  spin-axis-vertical 

(SAV)  and  spin-axi s-horizontal  (SAH)  orientations  of  the  grad- 

iometers.  (Results  presented  in  Sections  2.5  and  2.6  show 

that  similar  formulas  also  apply  to  the  Bell  Ball-Bearing  GGI 

and  the  Draper  Floated  Gravity  Gradiometer  No.  2.)  Both  the 

inline  and  cross-channel  noise  have  these  spectra.  The  posi- 

2 

t  ive  constants  R^^  and  R^.^  have  the  units  of  E  •  Hz  and 
determine  the  strength  of  the  "red"  noise*  at  lower  frequencies. 


The  term  "red”  noise  refers  to  noise  with  a  power  spectrum 

proportional  to  l/i  .  This  is  a  natural  extension  of  the 
term  "pink"  noise,  which  is  established  nomenclature  for 
acoustical  signals  having  spectra  proportional  to  1/f. 


36 


The  positive  constants  W0...  and  W.,...  have  tne  units  E^/Hz  and 
determine  the  strength  of  the  white  noise  at  higher  frequencies. 
The  lower  frequency  limit  of  0.01  mHz  for  this  model  corresponds 
to  a  period  of  28  hours  and  is  determined  by  the  manner  in 
which  the  data  were  analyzed.  This  frequency  limit  is  consis¬ 
tent  with  the  fact  that  simulations  of  airborne  surveys  do  not 
require  lower  frequencies.  The  upper  frequency  limit  of  24 
mHz,  on  the  other  hand,  was  the  result  of  the  Butterworth  fil¬ 
ters  in  the  gradiometer,  which  cause  the  spectra  to  decrease 
sharply  above  this  frequency. 

2.4.2  Time-Domain  Markov  Model 

A  Markov  model  consistent  with  the  spectra  in  Eqs. 

2.4-1  is  shown  in  Fig.  2.4-1,  where  the  symbols  R  and  W  denote 
the  red  and  white  noise  parameters.  The  small  parameter  f  is 
required  for  this  model  to  yield  stationary  self  noise.  Its 
value  within  the  range  0  Hz  <  fQ  <  10  ^  Hz  is  arbitrary  for 


R  *>099 1 


Figure  2.4-1  Markov  Model  of  Self  Noise,  Validated  for 

Frequencies  Ranging  from  0.01  mHz  to  24  mHz 


t: 
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the  purpose  of  modeling  self  noise  in  the  frequency  range  above 
0.01  mHz  (periods  less  than  28  hours).  The  value  f  =  10  ^  Hz 
yields  a  minimum-variance  model  consistent  with  the  spectral 
estimates.  The  white  noise  driving  this  model  has  a  power 
spectral  density  equal  to  unity. 


Numerical  values  for  the  noise  parameters  R  and  W 
were  computed  by  individually  fitting  R  and  W  values  to  all  of 
the  TASC  spectral  estimates  (both  FFT  and  AR)  and  then  comput¬ 
ing  the  mean  values  of  R  and  W  for  both  SAV  and  SAH  orienta¬ 
tions.  The  results  for  the  Bell  baseline  GG1  are 


SAV  =  2'°  "  10'6 

2 

E  -Hz 

SAH  =  16  *  10 

2 

E  •  Hz 

'SAV  =  81  E2/Hz 

SAH  =  86  e2'/Hz 

The  resulting  power  spectral  densities  (FSD's)  computed  from 
Eqs .  2.4-1  are  shown  in  Fig.  2.4-2.  These  models  apply  equal¬ 
ly  t o  the  grad i omet er ' s  inline  and  cross  channels.  The  cross¬ 
spectrum  between  these  channels  is  modeled  to  be  zero  because 
of  the  small  coherency  ,n  the  FFT  and  AR  spectral  estimates. 


2.4.3  Estimated  Random  Drift 


A  model  for  random  drift  in  the  gradiometer  outputs 
is  implied  by  the  self-noise  power  spectra.  If  x(t)  denotes 
the  output  in  ebtvbs  at  time  t  of  either  the  inline  or  cross 
channel,  then  the  random  drift  over  an  interval  T  can  be  de- 
f  l ned  a s 


d  ( T )  =  x  (  t  +  T )  -  x  (  t  ) 


(2.4-3) 
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Figure  2.4-2  ARMA  Models  of  Bell  Gradiometer  Self-Noise, 
Based  on  Smoothed  FFT's  and  Least -Squares 
Autoregressive  Models 

Note  that  d(T)  is  the  accumulated  end-to-end  drift  over  the 
interval  T,  as  opposed  to  a  drift  rate .  It  can  be  shown  that 
the  expected  value  of  this  drift  is  zero,  and  its  variance 
ad2(T)  is 


°/(T)  =  2[R _ (0)  -  R _ (T)  ] 


(2.4-4) 


where  RXX(T)  -  E[x(t  +  T)  x(t)]  is  the  autocorrelation  func¬ 
tion  of  the  self  noise  x(t). 

When  this  expression  is  evaluated  for  the  Markov  model 
it  is  found  that  the  following  approximation  is  valid  for  T 
<  T  <  105  s 


(2.4-5) 


V(T>  '  ''"2t  K  +  2  ^'‘ave 

T  =  measurement  averaging  time  (20  s  for 
24  mHz  bandwidth) 

2 

where  R  is  the  red-noise  parameter  (E  *riz)  and  W  is  the  white- 
noise  floor  (EvHz);  see  Eqs.  2.4-2.  The  rms  random  drift 
a^(T)  is  plotted  as  a  function  of  the  drift  time  T  in  Fig. 
2.4-3.  The  larger  rms  drift  for  the  horizontal,  as  compared 
to  the  vertical  spin  axis,  corresponds  to  the  higher  level  of 
red  self  noise  produced  by  the  horizontal  orientation.  This 
figure  indicates  that  the  accumulated  rms  drift  over  the  dura¬ 
tion  of  a  5-hour  gravity  survey  would  be  3.0  E  for  SAV  and  4.4 
E  for  SAH.  If  the  bandwidth  of  the  gradiometer  were  increased 
to  50  mHz  (10  second  averaging  time),  then  these  values  would 
increase  to  4.2  E  for  SAV  and  5.3  E  for  SAH. 


P-50937 


Computed  RMS  Random  Drift  for  Bell  Baseline 
DC  I  Based  on  Sell -Noise  Models  for  |f|  <  24  mHz 
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F i gu  re  2.4-3 


2.4.4  Comparison  of  TASC,  BELL,  and  NSWC/DL 
Spectral  Estimates 


Bell  presented  averaged  and  smoothed  FFT  autospectra 
on  the  SAV  and  SAH  orientation  of  their  baseline  GG1  in  Ref.  1. 
In  addition,  Naval  Surface  Weapons  Center/  Dahlgren  (NSWC/DL) 
generated  an  averaged  FFT  autospectrum  for  SAH.  These  spectral 
estimates  and  the  TASC  estimate  agree  as  shown  in  Fig.  2.4-4. 
The  slightly  reduced  level  of  the  Bell  estimate  below  0.2  mHz 
may  be  due  to  their  use  of  a  different  data  set  as  compared 
to  NSWC/DL  and  TASC.  The  differences  between  estimates  below 
0.02  mHz  are  attributed  to  the  finite  length  of  the  data  time 
series  and  the  different  methods  that  were  used  in  processing 
the  data.  Since  these  differences  occur  at  frequencies  cor¬ 
responding  to  1-day  time  periods,  they  are  not  important  for 
airborne  survey  modeling. 
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Figure  2.4-4  Comparison  of  Spectral  Estimates 
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2.5 


NOISE  SPECTRA  FOR  THE  BELL  BALL-BEARING  GG1 


Self-noise  data  for  the  Bell  prototype  Ball-Bearing 
gradioraeter  recently  became  available.  These  data  were  spec¬ 
trally  analyzed  using  the  AR  technique,  and  the  resulting  self 
noise  autospectra  are  shown  in  Fig.  2.5-1  along  with  the  ana¬ 
lytic  model  for  the  Bell  baseline  GGI . 


PERIOD 


Figure  2.5-1  Comparison  of  Self-Noise  Autospectra  for  Bell 
Baseline  and  Ball-Bearing  GGI's 


The  Bell  Ball-Bearing  GGI  (BB  GGI)  was  operating  with 
its  spin  axis  inclined  55  degrees  with  respect  to  the  vertical. 
This  is  the  intended  orientation  for  future  moving-base  appli¬ 
cations.  The  BB  GGI  data  spanned  six  days  with  a  sampling 
frequency  of  0.5  Hz.  To  provide  a  low-frequency  noise  spec¬ 
trum,  these  data  were  decimated  128:1,  which  yielded  2074  data 
pairs  at  a  reduced  sampling  frequency  of  3.91  mHz.  A  ramp  was 
removed  from  each  channel  of  data  (0.92  E/hr  in  the  inline 
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channel  and  0.30  E/hr  in  the  cross-channel),  and  an  AR  model 
was  fitted  to  the  residuals  by  minimizing  the  sum  of  the  squares 
of  the  one-step  prediction  errors  (Ref.  15).  The  AIC  statistic 
(Refs.  18  through  20)  dictated  an  order  9  model  (out  of  a  maxi¬ 
mum  permitted  order  of  24).  In  Fig.  2.5-1,  the  resulting  auto¬ 
spectra  for  inline  and  cross-channel  noises  of  the  BB  GG1  show 
small  ripples  at  high  frequencies  that  are  typical  of  high-order 
AR  spectra  with  white-noise  floors.  The  BB  GGI  has  a  white-noise 
floor  that  is  about  four  times  higher  than  that  of  the  baseline 
GGI,  which  corresponds  to  twice  as  much  short-term  variability 
on  an  rms  basis.  The  rising  low-frequency  parts  of  the  spectra 
indicate  more  long-term  drift  in  the  inline  channel  as  compared 
to  the  cross  channel  of  the  BB  GGI.  At  frequencies  between 
0.01  mHz  (1  day  period)  and  0.3  mHz  (1  hour  period),  the  in¬ 
line  and  cross-channel  noise  spectra  of  the  BB  GGI  fall  between 
the  baseline  spectra  for  vertical  and  horizontal  spin  axes. 

For  reference  purposes,  the  spectrum  level  at  the 
bottom  of  Fig.  2.5-1  corresponds  to  the  white-noise  level  of  a 
fictitious  "goal"  gradiometer  that  is  free  of  long-term  drift 
and  has  an  rms  self  noise  of  1.0  eotvos  with  a  10  second  aver¬ 
aging  time  between  measurements. 


2.6 


SELF-NOISE  MODEL  FOR  THE  DRAPER  FLOATED  GRADIOMETER 


Noise  data  were  spectrally  analyzed  for  the  Draper 
Floated  Gravity  Gradiometer  (FGG  No.  2).  An  overview  of  the 
largest  data  set  is  shown  in  Fig.  2.6-1,  where  the  output  sig¬ 
nals  of  the  FGG  are  resolved  into  independent  cross  gradients, 
T yZ  and  these  are  defined  by  the  gradiometer ' s  axis  of 

weights  and  two  orthogonal  axes,  denoted  Z  and  L  in  Draper's 
documentation  (Ref.  21).  Fig.  2.6-1  shows  the  residuals  after 
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Figure  2.6-1  Draper  FGG  No.  2  Self-Noise  Data 
(9  July  Tape) 


ramps  were  subtracted  from  the  data:  0.39  E/hr  in  the 
channel  and  0.58  E/hr  in  the  F^j  channel. 

The  data  span  13  days,  and  the  obvious  periodic  com¬ 
ponents,  evidenced  by  peaks  and  troughs,  correlate  with  normal 
work  hours,  weekends,  and  holidays  over  the  period  from  June 
23  through  July  4,  1979.  Larger  spikes  in  the  data  are  at¬ 
tributed  by  Draper  to  known  mass  disturbances  in  the  labora¬ 
tory  environment.  Draper  has  also  indicated  that  the  sensi¬ 
tivity  of  the  FGG  No.  2  to  environmental  temperature  cycles 
may  have  been  impacted  by  the  fact  that  two  gradiometer  tem¬ 
perature  controllers  were  inoperative. 

Analyses  by  means  of  FFT  and  AH  techniques  indicated 
that  the  self-noise  of  the  FGG  No.  2  can  be  represented  ana¬ 
lytically  with  the  following  model 
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SWL<f>  "  SWZ(f>  =  p-  *  W 

R  =  2.0  x  10-7  E2*Hz  (2.6-1) 

W  =  2.3  E2/K*. 

0.01  mH2<  J  f |  <50  mHz 

This  noise  spectrum  is  compared  with  analytic  models  for  the 
Bell  baseline  GGI  in  Fig.  2.6-2.  The  Draper  white-noise  floor 
is  about  three  percent  of  Bell's;  this  corresponds  to  about 
1/6  as  much  high-frequency  noise  (short-term  variability)  on 
an  rms  basis.  A  slightly  larger  difference  between  the  two 
instruments  occurs  at  lower  frequencies. 

2 

The  Draper  white-noise  floor  falls  beneath  the  10  E  / 
Hz  level  of  the  fictitious  "goal"  gradiometer.”  The  bottom  of 
Fig,  2.6-2  represents  the  theoretical  lower  limit  on  self  noise 
in  the  Draper  gradiometer  due  to  Brownian  motion  (Ref.  22). 

The  rms  accumulated  random  drift  in  the  FGG  output 
signals  (due  to  self  noise)  is  compuLod  by  means  of  the  theory 
described  previously  in  Section  2.4.3.  For  a  gradiometer  band¬ 
width  of  24  mHz,  which  corresponds  to  the  20-second  averaging 
time  previously  used  for  the  Bell  GGI,  the  drift  variance  over 
a  time  interval  T  is 


od2(T)  =  4  n2T  R  +  0.096  W 
20  s  <  T  <  105  s 


(2.6-2) 


*The  goal  gradiometer  is  arbitrarily  defined  to  have  an  rms 
error  of  1  E  when  the  output  is  averaged  for  10  seconds.  It 
is  described  further  in  Section  4.6. 
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Figure  2.6-2  Comparison  of  Self-Noise  Autospectra  for 

Draper  FGG  No.  2  and  the  Bell  Baseline  GGI 


R  Si  250 


Figure  2.6-!l  Computed  RMS  Random  Drift  for  Draper  FFG 
No.  2,  Based  on  Self -Noise  Model  for 
| f |  <  24  mHz 


The  rms  drift  a^d)  is  plotted  in  Fig.  2.6-3  for  the 
self-noise  model  in  Eq .  2.6-1.  This  plot  indicates  an  rms 
random  drift  of  0.60  E  over  a  5-hour  survey. 
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3 .  VIBRATION- INDUCED  ERRORS 

3.  1  INTRODUCTION 

This  chapter  describes  vibration  error  models  for  the 
Bell  Gravity  Gradiometer  Instrument  (GGI),  which  are  based  on 
analyses  of  test  data  that  extended  over  the  last  two  years. 
Vibration  test  data  for  the  Draper  Floated  gradiometer,  in 
contrast,  have  only  recently  been  available.  As  a  result, 
vibration  models  for  the  Draper  gradiometer  are  not  as  fully 
developed  as  those  for  the  Bell  instrument,  although  prelimin¬ 
ary  results  are  related  to  survey  performance  in  Chapter  4. 

This  chapter  is  organized  as  follows.  In  Section  3.2 
a  general  model  for  vi brat i on- i ndueed  errors  is  described  for 
the  Bell  GGI.  This  model  is  used  in  Sections  3.3  and  3.4  to 
interpret  the  Bell  data  on  vibration  sensitivities  and  to  de¬ 
velop  error  models  for  the  effects  of  aircraft  vibrations  such 
as  those  occuring  during  airborne  gravity  surveys.  Section 
3.5  describes  models  for  the  vibration  environment  of  an  air¬ 
borne  gradiometer  during  a  gravity  survey.  These  models  are 
ext rapolat ions  based  on  available  flight  test  data.  Upper 
bounds  for  vibrat ion- induced  errors  during  an  airborne  survey 
are  estimated  in  Section  3.6.  The  chapter  closes  with  Section 
3.7,  which  describes  an  error  model  for  an  airborne  survey 
system  consisting  of  three  gradiometers  mounted  on  a  single 
inert i a  1  plat  form . 
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3.2 


GENERAL.  ERROR  MODEL.  FOR  THE  BEI.I.  GRAVI  TY  GRAD  1  OME'J'ER 


An  ideal  Bell  gravity  gradiometer  is  completely  in¬ 
sensitive  to  linear  vibrations  as  well  as  to  angular  vih>ra- 
tions  about  its  spin  axis.  The  only  vi brat  ion- induced  errors 
are  those  caused  by  angular  vibrations  about  axes  normal  to 
the  spin  axis.  (These  vibrations  cause*  centripetal  gravity 
gradients  that  cannot  be  distinguished  by  the  gradiometer  from 
naturally  occuring  variations  in  the  earth's  gravity  field.) 

This  ideal  behavior  is  approximated  in  the  Bell  GGI 
by  the  use  of  on-line  error-correcting  feedback  loops  built 
into  the  gradiometer  (Refs.  3  and  26).  These  loops  provide 
corrections  for  accelerometer  scale  factor  unbalances,  mis¬ 
alignments  of  input  axes,  and  even-order  error  coefficients 
(bilinear  nonlinearities).  They  also  provide  correction  lor 
jitter  in  the  spin  rate. 

A  time-domain  model  of  the  input-output  behavior  of 
the  Bell  gradiometer  is  shown  in  Fig.  3.2-1.  The  part  of  this 
model  contained  within  the  dotted  box  represents  an  error- free 
gradiometer,  i.e.,  one  free  of  self  noise,  mi  sea  1 i brat  ion  or 
misalignment,  and  free  of  environmental  disturbances.  From  a 
mathematical  viewpoint,  the  error- free  system  has  the  same 
structure  as  a  communications  system  employing  coherent  ampli¬ 
tude  modulation  and  demodulation.  The  gravity  gradient  compo¬ 
nents  Fyy(t)  -  rxx(t  )  and  T  (t)  play  the  roles  of  modulation 
signals,  while  the  sinusoids  at  twice  the  spin  frequency  f 
play  the  roles  of  inphase  and  quadrature  carriers.  1’he  syn¬ 
chronous  detector  is  phase  locked  to  the  gradiometer  spin, 
which  has  a  frequency  of  about  0.28  Hz.  This  value  feu  f  is 
much  higher  than  the  effective  band  limits  of  either  the  gra¬ 
dient  signals  or  the  lowpass  output  filters  in  the  detector. 
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As  result,  the  inline-channel  output  Fj(t)  and  cross-channel 
output  l^(t  )  n(  the  error-free  instrument  are  lowpass  filtered 
■(■  rsions  of  the  input  gravity  gradient  signals: 
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Measurement  errors  induced  by  jitter  and  linear  vi¬ 
bration  of  the  gradiometer  are  also  indicated  in  Fig.  3.2-1. 

The  jitter  is  characterized  by  the  angular  velocity  vector 
Q(t)  measured  in  radian/second  (rad/s),  while  a(t)  is  an  ac¬ 
celeration  vector  that  characterizes  the  linear  vibration  in 
milli-g  (mg).  In  studies  of  airborne  gravity  surveys,  Q(t) 
and  a(t)  represent  the  vibration  environment  of  the  gradiometer. 
Normal]  Q ( t )  is  expected  to  be  dominated  by  jitter  in  the 
inertial  platform  that  supports  the  gradiometer,  and  a(t)  is 
primarily  driven  by  the  structural  vibrations  and  low-frequency 
motions  of  the  aircraft. 


A  model  for  the  vibration- induced  measurement  errors 
is  shown  in  Fig.  3.2-2.  The  measurement  error  e(t)  is  modeled 
as  the  sum  of  three  signals: 


•  Linear  jitter  error,  caused  by  rotational 
vibrations  near  harmonics  of  the  spin 
frequency.  The  gradiometer  output  error 
is  proportional  to  the  angular  velocity 
of  the  jitter. 

•  Quadratic  jitter  error,  caused  by  centri¬ 
petal  gradients  due  to  rotational  vibra¬ 
tions  about  axes  normal  to  the  spin  axis. 
The  gradiometer  output  error  is  propor¬ 
tional  to  the  square  of  the  angular  ve¬ 
locity  of  the  jitter. 

•  Linear-vibration  error,  caused  by  linear 

( translational )  vibrat ions  near  harmonics 
of  the  spin  frequency.  The  gradiometer 
output  error  is  proportional  to  the  linear 
acceleration . 


In  Fig.  3.2-2  the  linear  jitter  and 
sitivities  of  the  gradiometer  are  modeled  by 
gains  Gj  j  and  Gj ^ .  If  the  gradiometer  were 

*  1  g  =  9.81  m/s^ . 


v i brat i on  sen- 
t he  two  matrix 
free  of  linear 
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jitter  and  vibration  errors,  these  gain  matrices  would  be 
lilled  with  zeros,  and  the  induced  measurement  error  e(t) 
would  equal  only  the  quadratic  jitter  error  eQ ( t ) . 

The  linear  jitter  error  e^j(t)  is  a  linear  combina¬ 
tion  of  the  angular  velocities  in  the  jitter  vector  Q(t): 


eLj(  1  ^  M  Usl  }  °LJ 


(3.2-2) 


where  ui  =  2n (  ,  C,  ,  is  the  matrix  of  gains  that  model  the 
s  s  I..I  .j. 

grad i ome t e r ' s  linear  sensitivity  to  jitter,  and  M  (m^t  )  is  the 
transpose  of  th<-  modulation  victor 
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(3.2-3) 


This  modulation  vector  contains  inphase  and  quadrature 
harmonics  of  the  gradiometer  spin,  through  order  2.  The  phys¬ 
ical  basis  of  this  model  for  linear  jitter  error  has  been  do¬ 
cumented  in  Bell  technical  reports  (Refs.  3  and  26).  For  the 
purposes  of  characterizing  input-output  behavior,  specific  phys¬ 
ical  interpretations  of  the  individual  elements  of  ,  GjV, 
and  M(wgt)  are  not  required.  The  net  effect  on  the  performance 
of  the  gradiometer  is  that  certain  components  of  the  jitter 
and  vibration  are  modulated  by  the  spi n- frequency  harmonics 
because  the  accelerometers  in  the  gradiometer  are  attached  to 
an  instrument  block  that  rotates  at  the  spin  frequency  f  . 


An  equation  of  identical  form  describes  the  linear-* 
vibration  error: 


eLV(t)  =  MT(mst  )  Glv  a ( t  )  (3.2-4) 

In  this  equation,  a(t)  is  the  vibration  vector  of  linear  accel¬ 
erations  and  G^jj  is  the  matrix  of  gains  that  model  the  gradio- 
meter's  sensitivity  to  linear  vibration. 

The  scalar  output  e(t)  of  the  error  model  in  Fig.  3.2-2 
generally  has  a  broadband  spectrum,  but  only  the  spectral  compo¬ 
nents  of  e(t)  that  lie  within  24  mHz  of  2f.  (twice  the  gradi¬ 
ometer  spin  frequency)  contribute  significant  errors  to  the 
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outputs  of  the  synchronous  detector.  This  frequency  selec¬ 
tivity  is  characteristic  of  narrowband  AM  modulation  and  de¬ 
modulation.  The  selectivity  implies  that  the  input-output 
behavior  of  the  error  model  can  be  identified  from  laboratory 
data  obtained  by  shaking  and  rotating  the  gradiometer  with 
sinusoidal  test  signals  having  frequencies  equal  to  f  and  its 
first  four  harmonics.  Bell  performed  such  experiments  in  ac¬ 
cordance  with  formal  test  procedures  (Refs.  5  and  8).  A  method 
for  using  these  vibration  test  data  to  characterize  the  output 
of  the  error  model  for  random  inputs,  such  as  the  vibration  en¬ 
vironments  in  airborne  gravity  surveys,  is  described  in  the 
next  two  sections. 


3.3  ERRORS  INDUCED  BY  LINEAR  VIBRATIONS 

3.3.1  Interpretation  of  Bell  Test  Data 

The  Bell  tests  for  linear-vibration  sensitivity  are 
discussed  first.  During  these  experiments  the  gradiometer  in 
both  spin-axis-horizontal  (SAH)  and  sp i n- axi s-vert ical  (SAV) 
orientations  was  subjected  to  linear  vibrations  along  specified 
directions.  The  1 i nea r-accelera t ion  vector,  therefore,  ideally 
had  the  form 


a(t)  =  a  cos  (2nf  t  +  0)  a  (3.3-1) 

p  a 

where 


f  is  the  frequency  in  hertz 

is  the  peak  acceleration  in  milli-g 

a  is  the  unit  vector  that  defines  the  direction  of 
vibration  (along  one  of  the  principal  axes  of 
t  he  grad i omet er ) 

is  the  phase  angle  in  radians 
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The  vibration  frequency  f  was  chosen  to  equal  different  low- 

a 

order  harmonics  of  the  gradiometer  spin  frequency  fg.  In  ac¬ 
cordance  with  the  error  model,  the  resulting  steady-state  out¬ 
puts  of  the  detector  for  f  =  f  ,  2f_,  3f  ,  and  4f  can  be 
shown  to  be  of  the  form 


rT(t)  =  constant  =  C  cos  (t|i) 

1  (3.3-2) 

rc(t)  =  constant  =  ±C  sin  (<Jj) 

where  C  and  ip  are  real  numbers.  The  phase  angle  tj/  depends  on 
the  phase  0  of  the  applied  vibrations,  which  was  held  fixed 
for  each  frequency  used  in  the  experiment. 


The  root-mean-square  (rms)  values  of  the  two  output 
signals  are  equal  to  each  other  when  the  average  is  over  all 
equally  probable  values  for  the  two  phase  angles  0  and  <(>  of 
the  linear  vibration  and  the  gradiometer  spin: 


(r1(t) 


rms 


=  lrC“> 


rms 


=  r 


rms 


(3.3-3) 


where 

rrms  =  lcWZ  <33-4> 

A  per-channel  rms  sensitivity  to  linear  vibration, 

denoted  SrTxr(a,f  ),  can  be  defined  for  each  vibration  direc- 
L1N  a 

tion  a  and  frequency  f  : 

3 


ar. 

SLlN(a,fa)  =  da 


rms 


rms 


(3.3-5) 


where  the  rms  acceleration  is 


a  =  a  /V2 
rms  p7  N 


(3.3-6) 
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Numerical  values  for  these  sensitivities  were  computed 
from  the  Bell  test  data  by  fitting  least-squares  straight  lines 
to  estimate  the  slopes  3Fj(t)/8ap  and  8r£.(t)/8a  .  The  sensi¬ 
tivities  were  then  computed  as 


^LIN^’^a^  *s  t'ie  ®ain  which  the  gradiometer 

demodulates  (heterodynes)  the  input  accelerations  near  fre¬ 
quency  f  to  produce  low-frequency  errors  in  the  output  of  the 
detector.  The  significance  of  this  sensitivity  for  random 
acceleration  disturbances  is  described  next.  It  will  be  s  wn 
that  random  accelerat i ons ,  such  as  occur  in  airborne  gravi  ■ 
surveys,  cause  white-noise  errors  in  the  gradiometer  output?. 


3.3.2  Mode  1 _ for  Errors  Induced  by  Random 

Linear  AceeTeraTlon 


Let  the  acceleration  vector  a(t)  be  a  random  process 
of  the  form 


a(t)=a(t)a  (3.3-8) 

where 

a  is  a  fixed  unit  vector  specifying  the  direction  of 
the  acceleration 

a(t)  is  a  stationary  random  scalar  process  with  power 
spec  t  rum 

S  (f)  =  f  H  it )  e'i2nft  dt  (3.3-9) 

,M  -00  rlcl 

N  ,(t)  is  the  autocorrelation  function  of  a(t) 

cl  <1 

K  (t)  =  l'.|  a  ( t  +t  )a(  t  )  |  (3.3-10) 

<i  cl 
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Also  let  the  phase  angle  <j)  of  the  gradioiueter  spin  be  a  random 
variable,  uniformly  distributed  on  the  interval  (0,2n). 


The  random  acceleration  a(t)  induces  low-frequency 
random  error  outputs  from  the  synchronous  detector.  According 
to  the  error  model,  the  power  spectrum  of  the  error  signal  at 
the  output  of  the  inline  channel  is  the  same  as  the  error  spec¬ 
trum  for  the  cross-channel  output  and  is  denoted  Spp(f).  This 
spectrum  is  the  sum  of  heterodyned  versions  of  the  acceleration 
spectrum,  weighted  by  the  squared  sensitivity  coefficients  for 
each  spin-frequency  harmonic  and  weighted  also  by  the  squared 
frequency  response  of  the  output  filter. 


Spp(f)  = 


HD(f) 


n?0 


fn> 


lSaa<f‘fn>  + 


S  (f+f  )] 
aa  n  J 

(3.3-11 ) 


where  fn  =  nfs  =  n*1^1  harmonic  of  the  spin  frequency,  and  Hpj(f) 
is  the  Fourier  transform  of  h^U). 


The  pass  band  of  the  detection  filter  Hp^f)  is  only 

24  mHz,  which  is  much  smaller  than  the  harmonics  f  of  the 

n 

spin  frequency.  Therefore,  as  long  as  the  shifted  spectra 
S  (f±f  )  are  free  of  resonances  in  the  band  |f|  5=  24  mHz,  the 

a  a  n 

expression  for  the  output  error  spectrum  can  be  approximated  as 


Spp(f)  s 


HD(  f ) 


4 

n?0 


SLIN(a,fn) 


Saa<fn> 


(3.3-12) 


which  simply  represents  white  noise  in  the  output  bandwidth. 
This  means  that  the  errors  induced  by  random  linear  vibrations 
are  equivalent  to  an  increase  in  the  white-noise  level  of  the 
grad i omet er ' s  self  noise. 
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The  analysis  is  extended  to  3-diinensional  anisotropic 
random  accelerations  as  follows.  Assume  that  the  acceleration 
vec  tor 


a(t)  =  ax(t)x  +  ay(t)y  +  az(t)z  (3.3-13) 

where  a  (t),  a  (t),  and  a  (t)  are  uncorrelated  zero-mean  sta- 
x  y  z 

tionary  random  processes  that  have  the  power  spectra  S  (f), 

S  (f),  and  S  (f),  and  x,  y,  z  are  unit  vectors  pointing  along 
y  y  zz 

the  principal  axes  of  the  gradiometer.  Then  the  white  noise 
induced  in  the  two  detector  outputs  both  have  the  same  power 
spectrum: 


srr(f) 


nD(f) 


HDror 


2  n?0  lSxx(fn)  SLIN(X’  fn} 


+  Syy(tn)  SL1 N( ^ ’  1n)  +  Szz(fn*  SLlN(z’  fn)J 


(3.3-14) 


The  value  of  the  sensitivity  Sj  (a,f^)  depends  on 
the  matrix  gain  G(  v  in  the  error  model.  Fortunately,  each 
element  of  G does  not  have  to  be  known  for  the  purposes  of 
estimating  output  errors  due  to  stationary  random  inputs;  the 
available  test  data  are  sufficient  for  developing  a  full  3- 
dimensional  error  model.  However,  such  a  model  is  based  on 
the  assumption  that  the  gradiometer  sensitivity  to  linear  vi¬ 
bration  is  isotropic  about  the  spin  axis.  This  assumption  is 
appropriate  because  practical  limitations  of  the  test  stand 
used  in  the  Bell  tests  made  it  impossible  to  apply  vertically- 
directed  accelerations  to  the  gradiometer. 
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3.3.3  Estimated  Linear  Vibration  Sensitivities 


A  set  of  estimated  linear  vibration  sensitivities  for 
the  GG1  is  listed  in  Table  3.3-1.  (Data  were  not  available 
for  sensitivities  at  the  4th  spin- frequency  harmonic.)  Sensi¬ 
tivities  at  higher  frequencies  were  expected  to  be  small  com¬ 
pared  to  the  sensitivities  in  Table  3.3-1.  This  was  verified 
by  test  data  taken  at  4  Hz,  8  Hz,  and  16  Hz.  However,  at  30 
Hz  an  unexplained  larger  sensitivity  of  5  E/mg  was  observed 
for  the  SAH  orientation  when  the  vibration  was  normal  to  the 
spin  axis.  This  anomaly  is  attributed  by  Bell  to  a  structural 
resonance  in  the  test  stand  or  the  structure  of  the  gradiometer 
instrument.  With  regard  to  future  testing  of  ball-bearing 
gradiometers ,  this  suggests  that  the  fundamental  vibration 
modes  of  the  test  stand  (and  their  possible  interaction  with 
the  gradiometer)  be  identified  to  permit  the  clearest  practi¬ 
cal  understanding  of  gradiometer  vibration  sensitivities. 


TABLE  3.3-1 

ESTIMATED  LINEAR  VIBRATION  SENSITIVITIES 
FOR  BELL  GG1  NO.  1 


SPIN  AXIS 
ORIENTATION 

VIBRATION 
DIRECTION 
WITH  RESPECT 

LINEAR  VIBRATION  SENSITIVITY 
( E/mg ) 

VIBRATION  FREQUENCY 
(f  =  SPIN  FREQUENCY  ~  0.28  Hz) 

■ 

TO  SPIN  AXIS 

o 

f 

s 

2f  |  3f 

s  s 

Horizontal 

Paral lei 

mi 

0.13 

0.79 

1 . 68 

(SAH) 

Normal 

0.081 

1 

1  .17 

Vert ical 
(SAV) 

Norma  1 

0.4(»t 

4 . 30 

4.40 

Paral lei 

0.1  :C 

m 

1  .  68" 

"'Not  available  due  to  test  geometry  limitations  -  modeled  con 
si  stent  ly  with  SAI!  response. 

tl'tom  tilt  test  data  on  static  g  sensitivity. 
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LRRORS  INDUCED  BY  ANGULAR  VIBRATIONS 


P 


5 .  6 


<.6.1  1 n t e rp re t at i on  o f  Bell  Test  Data 


Bell  conducted  formal  tests  (Refs.  8,  27)  to  measure 
the  angular  vibration  sensitivity  of  their  baseline  GG1 .  For 
this  purpose  sinusoidal  angular  vibrations  (jitter)  were  ap¬ 
plied  about  the  principal  axes  of  the  gradiometer,  and  the 
result  i ng  error  signals  at  the  detector  outputs  were  measured. 
These  data  can  be  interpreted  with  the  help  of  the  error  model 
m  Fig.  <.1-1  bv  expressing  the  angu 1 ar-veloe i t y  vector  Q ( t ) 
in  the  1 1  >  rin 


o(t)  r  U  cos  ( 2/i  f  t  +  0)1/  (3.6-1) 

p  a 


where 


is  the  peak  angular  velocity  in  radians  per  second, 

is  the  i remnncv  in  hertz,  0  is  the  phase  angle  in  radians, 
a 


.mil 


is  <i  unit  vector  pointing  along  the  fixed  axis  of  the 
a p p lied  rotation.  1 ht 
tor  h . 


•  x,v,x  components  of  this  velocity  vec- 
ive  t  he  peak  va 1 ues 


PX 


py 


pz 


U  ■  X 

p. 

V'v 

u  u  ■  7. 

p 


(3.6-2) 


where  x,  y,  ’/.  are  unit  vectors  parallel  to  the  x,  y,  and  z 
axes.  (The  7.  axis  is  chosen  to  coincide  with  the  ^radiometer 
spin  axis.  )  In  the  formal  tests  considered  here,  Q  was  point¬ 
ed  along  principal  axes  of  the  gradiometer,  1  (  was  a  harmonic 
of  the  gradiometer  spin  lieipiency  I  ,  and  0  was  a  constant 
phase  angle  for  each  choice  ol  1 


The  error  model  includes  both  quadratic  and  linear 
errors  that  are  induced  by  the  jitter.  The  linear  measurement 
error  e^j(t)  computed  first,  in  the  same  way  the  linear 
vibration  error  e^y(t)  was  analyzed  in  the  lasL  section;  the 
quadratic  error  eQj(t),  which  is  due  to  centripetal  gradients, 
is  analyzed  in  Section  3.4.3. 


From  Fig.  3.2-2,  the  j i tter- induced  linear  measurement 

error  is 

eLJ(t)  3  «T<V>  GLJ  ^(t)  (3.4-3) 

The  resulting  steady-state  error  outputs  of  the  de¬ 
tector  are  constants  when  f  is  a  spin- frequency  harmonic,  and 

ct 

the  rms  values  of  these  constants  (averaged  over  all  equally 
probable  phase  angles  0  and  4> )  satisfy  equations  identical  in 
form  to  those  derived  for  linear  vibrations  in  the  previous 
section.  Specifically, 


If  j(t)  ] 


rms 


(rc(t) 


rms 


T 

rms 


(3.4-4) 


where  the  value  of  r  depends  on  the  jitter  frequency  f  and 
the  jitter  axis  fi.  A  per-channel  rms  sensitivity  to  jitter 
can  be  defined  in  the  same  fashion  as  was  done  for  linear  vi¬ 
bration  in  Eq.  3.3-5: 


ar 

s  ( o  f  )  =  - _ rgg 

aANG'w,Ia;  BO 

rms 


(3.4-5) 


where  the  rms  velocity  is 


Q 


rms 


=  VV2 


(3.4-6) 
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Numerical  values  for  these  angular  sensitivities  were 
computed  from  the  Bell  test  data  by  fitting  least -squares  straight 
lines  to  estimate  the  slopes  3rT(t)/3f)  and  3rr(t)/3Q  .  The 
sensitivities  were  then  computed  as 


3.4.2  Es  t ima t ed  Angu lar- Vibration  Sensitivities 

The  angular  sens i t i vi t i es  computed  from  the  available 
test  data  are  sufficient  for  developing  a  full  3-dimensional 
model  for  j i t t er- i nduced  gradient  measurement  errors.  In  build¬ 
ing  this  model  it  is  assumed  that  the  gradiometer ' s  response 
to  jitter  is  both  isotropic  about  the  spin  axis  and  independent 
of  spin  axis  orientation.  These  assumptions  are  appropriate 
because  the  design  of  Bell’s  test  stand  permitted  application 
of  rotations  only  about  a  horizontal  axis.  A  full  set  of  jit¬ 
ter  sensitivities  is  shown  in  Table  3.4-1. 


TABLE  3.4-1 

ESTIMATED  JITTER  SENS1TIVITES  FOR  BELL  GGI  NO.  1 


SPIN  AXIS 
ORIENTATION 

VIBRATION 
DIRECTION 
WITH  RESPECT 

TO  SPIN  AXIS 

JITTER  SENSITIVITIES 
(10A  E/( rad/ s ) ) 

f 

s 

2fs 

3L  

4fs 

Horizontal 

Paral lei 

9.2 

4.5 

8.3 

7.6 

(SAH) 

JL 

Norma  1 

9.7 

2.4 

66 

35 

Ve  rt i ca  1 

Paral lei t 

9.2 

4.5 

8.3 

7.6 

(SAY) 

No rma 1 

9.7 

2.4 

66 

35 

-'•Linear  drift  removed  from  raw  data. 

tNot  available  due  to  test  geometry  limitations  -  modeled  con¬ 
sistently  with  SAH  response. 

“•'Data  not  available  -  modeled  consistently  with  SAH  response. 


62 


Some  of  the  Bell  data  that  were  used  to  compute  the 
sensitivities  in  Table  3.4-1  exhibited  linear  drift,  which  was 
subtracted  from  the  data  before  estimating  the  sensitivities. 
Bell  attributed  this  drift  to  the  fact  that  resonance-hunting 
experiments  preceded  the  jitter  tests  at  spin- frequency  har¬ 
monics.  As  a  result,  one  of  the  error-correcting  feedback 
loops  in  the  gradiometer  was  over-loaded,  and  this  loop  had 
not  recovered  fully  when  the  jitter  tests  were  performed.  The 
jitter  data  therefore  exhibited  drift  corresponding  to  the 
transient  response  of  the  feedback  loop.  This  drift  would 
not  occur  during  normal  gradiometer  operation. 

3.4.3  Model  for  Errors  Induced  by  Random 
Angular  Velocities 


The  power  spectra  of  errors  caused  by  3-dimensional 
stationary  random  jitter,  which  is  used  later  to  model  air¬ 
craft  vibrations  during  airborne  gravity  surveys,  can  be  com¬ 
puted  by  using  the  same  reasoning  used  in  the  previous  section 
for  linear  vibration.  The  jitter  vector  is  represented  by 

«(t)  =  Qx(t)x  +  Oy(t)y  +  Qz( t)z  (3.4-8) 

where  fi  (t),  Q  (t),  and  fi  (t)  are  uncorrelated  zero-mean  sta- 
x  y  z 

tionary  random  processes  that  have  the  power  spectra  S  (f), 

XX 

S  (f),  and  S  (f).  Then  the  induced  errors  in  the  two  detee- 
yy  zz 

tor  outputs  both  have  the  same  power  spectrum: 


where  I  '  nt  1'h  i  s  is  !  hr  sum  ol  the  contributions  from  the 
II  s 

x,  v.  and  /.  components  taken  separately. 

The  one  remaining  component  of  error  in  Fig.  3.2-2  is 
the  quadratic  jitter  error  e^^(t)  caused  by  centripetal  gradi¬ 
ents.  This  error  source  was  negligible  during  the  Bell  tests 
conducted  at  spin- 1  requency  harmonics  and  was  therefore  omitted 
in  the  preceding  discussion.  However,  in  other  tests  at  high 
frequencies  (A  Hz,  10  Hz,  20  Hz,  30  Hz,  00  Hz,  and  50  Hz)  the 
quadratic  jitter  error  was  the  dominant  error  source. 


According  to  Figs.  3.2-1  and  3.2-2,  the  inline  and 
cross -channel  outputs  due  to  quadratic  jitter  error  during 
these  high-frequency  tests  are 


-  $r. 

f  ,(t  )  =  px  7t  p>  Hu(0) 
I  (;(t  )  =  Wp2*  py  llj/0) 


(3.0-10) 


.  u 

where  H.(0)  -  10  F/(rad/s)^.  The  Bell  data  correspond  to  the 

,  2  2 

case  where  0  '>  0  ~  0,  which  implies  that 

py  px  K 


l'  i<>)|  >>  | rc(.  >|  '  0 

l 1 1  < 1  >1  ; 


(3.0-11) 


These  relations  are  verified  by  the  Bell  high-frequency 
test  data.  For  example,  in  Fig.  3.0-1  the  magnitude  of  the 
i n 1 i ne-channe 1  output  1  j(t)  for  jitter  normal  to  a  vertical 
spin  axis  is  plotted  as  a  function  ol  t he  peak  angular  velo¬ 
city  0  .  lln-  dotted  line  is  the  upper  bound  for  these  data 

as  predicted  by  l.q .  d.O.-ll.  The  fact  that  the  data  are  grouped 
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Figure  3.4-1  Bell  Jitter  Data  t SAV ,  Bell  Tape  22) 

near  and  beneath  this  line  indicates  agreement  with  the  error 
model.  Figure  3.4-2  shows  similar  results  for  jitter  normal 
to  a  horizonal  spin  axis. 


For  gravity  gradiometer  survey  simulations  it  is  nec¬ 
essary  to  compute  the  power  spectrum  of  the  gradient  errors 
caused  by  random  quadratic  jitter.  To  do  this  the  Q  (t)  and 
fty(t)  jitter  components  are  modeled  as  independent  jointly- 
stationary  zero-mean  Gaussian  random  processes  with  power  spec 

tra  S .... (f)  and  S  At)  and  standard  deviations  o  and  a  .  Ac- 
xx  yy  x  y 

cording  to  Fig.  3.2-2,  the  output  error  signals  due  to  quad¬ 
ratic  jitter  are  accurately  approximated  bv 


scc(r>  -  isxx(f)*syy<m  inD(f)i2 

sIC(f)  =  0 

Hd(0)  =  JO9  £/( rad/s)2 
6(f)  =  delta  function 

The  formulas  for  the  autospectra  can  be  simplified 
with  a  negligible  loss  in  accuracy  when  the  spectra  S  (f)  and 
S  (f)  are  analytic  functions.  The  approximat i ons  are  derived 
by  expanding  the  integrands  of  the  convolutions  in  powers  of  f 
and  then  retaining  only  the  first  terms.  This  procedure  is 
appropriate  when  the  lowpass  filter  H^(f)  has  a  small  bandwidth 
compared  to  the  scale  of  the  spectra  Sxj.(f)  and  Syy(f).  For 
example 

00 

S  (f)*S  (f)  =  J  S  (f-f')S  (f')df' 

xx  yy  xx  yy 

90 

=  J  S  (f'-f)S  (f')df' 

xxv  yy v 

OO 

=  /  [S  (f')  -  S'  (f')f  +  S  (f  )f2/2  - .  . .  ]S  (f')df' 
i,  xx  xx  xx  yy 

—  00  J  J 

00  00 

=  J  s  (f')S  (f')df'  -  f  J  S'  (f')S  (f')df'  +  0(f2) 

_i  xx  yy  xxv  yy 

(3.4-14) 

2  2 

where  0(f  )  denotes  terms  proportional  to  f  and  higher  powers 

of  f.  For  f  sufficiently  small,  only  the  first  term  is  sig¬ 
nificant,  and  this  is  the  term  that  is  retained  in  the  approx¬ 
imation.  The  resulting  simplified  formulas  for  inline  and 
cross-channel  error  spectra  are 
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Su(f)  -  Hj*(0)  |o2  -  oy  J26 ( f )/4 

Xj4(f)  +  S2y(f))df  (3.4-15) 

Scc(f)  s  |HD(f)j  2_/\x(f)Syy(f)df 


HD(f) 


These  equations  state  that  the  quadratic  jitter  er¬ 
rors  produce  independent  lowpass  white-noise  processes  in  the 
two  output  channels  of  the  gradiometer. 

These  approx imat ions  are  valid  for  jitter  described 
by  continuous  power  spectra,  such  as  the  broadband  spectra 
used  later  in  this  report  to  model  the  vibration  environment 
of  the  gradiometer  in  an  airborne  gravity  survey.  For  narrow- 
band  jitter  containing  spectral  delta  functions,  the  above 
approximation  is  invalid.  This  case  can  actually  aris°  in 
practice  because  the  inertial  platform  that  supports  the  gradi¬ 
ometer  may  have  stietion,  which  is  expected  to  result  in  narrow- 
band  spectral  components  at  the  resonance  frequency  of  the 
platform  and  its  harmonics.  An  analysis  of  platform  jitter  is 
given  in  Appendix  A  where  its  practical  impact  on  gradiometer 
performance  is  discussed.  It  is  concluded  that  uncompensated 
platform  jitter  could  be  a  dominant  error  source  in  airborne 
surveys.  However,  compensation  of  the  gradiometer  outputs  is 
feasible  (via  gyro  measurements  on  the  platform)  to  reduce  the 
jitter  errors  to  negligible  levels.  Such  conipensat  i  on  has  not 
yet  been  demonstrated  by  Bell. 
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3.5  MODELS  FOR  THE  VIBRATION  ENVIRONMENT  OF  AN  AIRBORNE 
G RADIOMETER  DURLNG  GRAVITY  SURVEYS 

3.5.1  Introduction 

In  the  previous  sections  it  was  shown  that  both  linear 
and  angular  vibrations  induce  errors  in  the  gradiometer  output 
signals.  These  errors  were  shown  to  consist  of  biases  plus 
white  noise  when  the  vibrations  are  stationary  random  process¬ 
es.  In  Section  3.5.2  of  this  chapter,  models  are  developed 
for  random  aircraft  vibration.  These  models  are  based  on 
available  flight  test  data.  Section  3.5.3  describes  models 
for  the  passive  isolation  system  and  inertial  platform  that 
support  the  gradiometer  in  the  aircraft. 

3.5.2  Linear  and  Angular  Aircraft  Vi  brat i on s 

Flight  Test  Data  -  A  variety  of  flight  test  data  on 
the  linear  and  angular  vibrations  of  military  aircraft  were 
evaluated  for  the  purposes  of  modeling  the  vibration  environ¬ 
ment  of  a  gradiometer  during  an  airborne  gravity  survey. 

Ideally  these  data  would  provide  power  spectra  for  the  x,  y, 
and  z  components  of  both  linear  and  angular  vibrations  at  the 
location  where  the  gradiometer  system  and  its  supporting  hard¬ 
ware  would  be  mounted.  These  spectra  would  ideally  span  the 
frequency  range  from  0.1  mHz  to  1000  Hz  and  would  correspond 
to  typical  survey  flight  conditions.  The  available  test  data 
that  came  closest  to  these  ideals  were  in  Ref.  2,  which  pro¬ 
vides  vibration  power  spectra  for  an  NKC-135A  aircraft  in 
"moderate  turbulence"  at  airspeeds  ranging  from  0.5  mach  to 
0.8  mach.  The  data  are  worst-case  spectra  obtained  from  two 
test  flights.  The  altitude  ranged  from  5,000  feet  to  37.000 
feet,  but  was  not  uniquely  specified  in  the  test  data.  In 
addition  the  low-frequency  limit  of  the  power  spectra  was 
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approximately  0.1  Hz.  These  data  were  used  to  generate  vibra¬ 
tion  models  for  survey  simulations  by  fitting  parametric  spec¬ 
trum  models  as  explained  in  the  following. 


Vibration  Models  -  The  parametric  vibration  models 
were  structured  at  low  frequencies  to  extrapolate  the  test  data 
down  to  zero  frequency  in  the  simplest  possible  way  consistent 
with  the  physical  constraint  that  the  linear  and  angular  dis¬ 
placements  should  have  finite  variances.  The  extrapolations 
provided  constant  spectral  density  values  in  the  neighborhood 
of  zero  frequency  for  the  linear-displacement  and  angular- 
displacement  spectra.  These  extrapolations  correspond  to  a 
low-frequency  slope  of  12  dB/octave  in  the  linear  acceleration 
spectrum  and  6  dB/octave  in  the  angular-velocity  spectrum. 

Both  the  linear  (vertical  acceleration)  and  angular  (pitch, 
roll,  and  yaw  rates)  test-data  spectra  had  trends  that  could 
be  modeled  well  by  the  following  rational  functions 

Sfift(f)  =  8  X  10'4  K  (f/°-2)2/l (f/0.2)4  +  1)  (rad/s)2/Hz 

S*  (f)  =  4  xio"2  K  (f/0.15)4/|(f/0.15)4  +  1  ]  |  (f/0. 1 5 ) 2  +  1)  g2/Hz 

a  3 

2 

K  =  1/30  =  scale  factor 

(3.5-1) 


A  A 

where  S^(f)  and  S  (f)  denote  the  spectra  for  the  angular 
rates  and  linear  acce lerat ions .  When  the  scale  factor  K  =  1, 
these  spectra  model  the  flight  test  data  corresponding  to  what 
is  termed  in  Ref.  2  as  worst-case  "moderate  turbulance."  Ac¬ 


cording  to  the  model  spectra,  the  linear  vibration  in  the  ver¬ 
tical  direction  then  has  an  rrns  acceleration  of  0.09  g,  an  rms 
velocity  of  0.7  m/s ,  and  an  rms  displacement  of  1  m.  Also  for 
K  =  1,  the  angular  vibration  about  the  pitch  axis  has  an  rms 
velocity  of  1  deg/s  and  an  rms  rotation  of  1  deg.  These 
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are  much  larger  than  the  vibration  levels  expected  during  a 

2 

gravity  survey;  the  scale  factor  value  K  =  1/30  reduces  these 
rms  levels  by  a  factor  of  30,  which  is  believed  to  be  a  rea¬ 
sonable  adjustment  for  survey  modeling. 


The  model  spectra  are  plotted  together  with  the  en¬ 
semble  of  available  flight  test  data  for  comparison  in  Fig. 

A 

3.5-1  and  3.5-2.  In  these  figures  the  quantities  2S^(f)  and 
2S^  (f)  (with  K  =  1)  are  plotted  for  easy  comparison  with  the 

a  a 

one-sided  test-data  spectra,  which  were  measured  with  sensors 
at  various  locations  on  the  aircraft. 


n— 509*5 


Figure  3.5-1  Ensemble  of  Flight  Test  Data  for  Pitch, 
Roll,  and  Yaw  Rates  and  Model  Spectrum 
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Figure  f.r)-2  Ensemble  of  Flight  Test  Data  for  Vertical 
Accelerations  and  Model  Spectrum 

These  test  data  display  a  high  degree  of  variability, 
and  they  correspond  to  worst-case  conditions  during  moderate 
turbulance  tor  a  single  aircraft.  Furthermore,  the  data  do 
not  include  lateral  and  longitudinal  linear  accelerations. 
Therefore,  the  vibration  model  is  isotropic  to  allow  a  full 
three-dimensional  representation  of  the  vibration  environment. 
That  is,  the  pitch,  roll,  and  yaw  rates  are  independent  random 
processes,  each  having  t  lie  power  spectrum  S^  (  f )  ,  w’hile  the 
linear  accelerations  along  the  vertical,  lateral,  and  longitudi¬ 
nal  axes  of  the  aircraft  are  also  independent  random  processes, 
each  with  the  power  spectrum  S,Al  (1). 
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3.5.3  G radiometer  Isolation-Syst cm  and  Inertial- 
Platform  Models 


In  projected  airborne  surveys  the  gradiometer  is 
mounted  on  an  inertial  platform,  which  in  turn  is  supported  by 
a  passive  vibration  isolation  system  attached  to  the  airframe 
structure.  Mechanical  stability  requirements  put  a  practical 
lower  limit  of  about  6  Hz  on  the  natural  frequencies  of  the 
passive  isolator.  For  simplicity  the  6  degrees  of  freedom  (3 
translational  and  3  rotational)  are  modeled  as  uncoupled  from 
each  other,  and  the  damping  constant  of  the  resonance  for  each 
mode  is  0.25.  The  transfer  function  Tjg(s)  of  such  an  isolator  is 


TIS<s) 


xout(s) 


=  0.25 


2£lSu)iss  +  w 

~j - 

S  +  t  T  Q® 


IS*ISS  +  WIS 


(3.5-2) 


ujjS  =  2n6  rad/s 

s  =  complex  frequency  variable 

where  xout(s)  anc*  x^n(s)  are  the  Laplace  transforms  of  the 
output  and  input  vibration  quantity,  which  could  be  pitch, 
yaw,  or  roll  rate,  or  vertical,  lateral,  or  longitudinal  line¬ 
ar  acceleration.  The  magnitude  of  the  frequency  response  of 
this  transfer  function,  obtained  by  setting  s  =  i2af,  is  shown 
in  Fig.  3.5-3. 

The  inertial  platform  provides  isolation  from  angular 
vibrations  and  is  described  by  a  transfer  function  Tpj(s)  that 
relates  input  angular  rates  to  the  resulting  output  angular 
rates  of  the  stabilized  gradiometer  with  respect  to  an  inertia] 
reference : 
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MAGNITUDE,  20  log10|  TpL(i  2*  f)|  (dB) 


Figure  3.5-3  Magnitude  of  Transfer  Function  Tlc.(s) 

for  Passive  Isolator  b 


VPL 

s  +  24 


PL^PL^ 

PL^PL5 


CFL  =  0.707 

Wpj  =  2n20  rad/s 

KpL  =  3.16  x  10"3 


(3.5-3) 


This  model  represents  commonly  used  platform  dynam¬ 
ics.  At  the  damped  20  Hz  resonance,  it  attenuates  angular  rates 
by  a  factor  of  520  (50  decibels).  The  magnitude  of  the  fre¬ 
quency  response  of  this  model  is  shown  in  Fig.  3.5-4. 
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Figure  3.5-4  Magnitude  of  Transfer  Function  Tpp(s) 
for  Inertial  Platform 


3.6  ESTIMATED  UPPER  BOUNDS  ON  VI BRAT1 ON- INDUCED  ERRORS 
DURING  AIRBORNE  SURVEYS 


3.6.1  Acceleration  and  Jitter  Spectra 

The  models  for  the  aircraft  vibrations  and  the  gradi- 
ometer  mounting  system  are  now  combined  to  yield  power  spectra 
for  modeling  the  vibration  disturbances  acting  on  a  gradiometer 
during  an  airborne  survey.  The  vertical,  lateral,  and  longi¬ 
tudinal  linear  accelerations  experienced  by  the  gradiometer 
each  have  the  power  spectrum 

saa(  f)  -  I T  j  s  (  i  2tx  f )  |  2  S^f)  (3.6-1) 


while  the  pitch,  roll,  and  yaw  rate  spectra  are 


73 


TpL(i2nf)Tls(i2nf)| 2  s£ft(f) 


(3.6-2) 


These  vibration  spectra  are  shown  in  Figs.  3.6-1  and 
3.6-2.  The  critical  spin- frequency  harmonics  of  the  gradiome- 
ter  are  indicated  by  the  arrows  in  the  0.25  Hz  to  1.0  Hz  range. 
These  are  the  frequencies  at  which  the  gradiometer  exhibits 
linear  vibration  and  linear  jitter  sensitivities. 


3.6.2  A  Bound  on  the  RMS  Error  Due  to  Jitter 


The  practical  significance  of  the  jitter  spectrum  in 
Fig.  3.6-2  is  discussed  first.  As  previously  explained  in 
Section  3.6.3,  jitter  at  any  frequency  induces  error  signals 
in  the  two  gradiometer  outputs.  For  the  jitter  model  dis¬ 
cussed  here,  which  excludes  slid  ion  in  the  inertial  platform, 
the  inline  and  cross-channel  error  spectra  due  to  quadratic 
jitter  errors  (centripetal  gradients)  are  given  by  Eqs .  3.6-15 
and  3.6-2  with  SXJ£(f)  =  Syy(f)  =  SQfi(f): 
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1 1 
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(f)  = 

1  .  1  X 

io-10 

EZ/Hz 

(f)  = 

1.1  >■ 

r-10 

E2/Hz 

m  < 

0  M 

*  *7 

(3.6-3) 


The  fluctuating  white-noise  component  in  each  channel  has  an 
rms  value  of  2.3  *  10  ^  K.  These  quadratic  errors  (due  to 
aircraft  jitter)  are  negligible  because  they  are  more  than  six 
orders  of  magnitude  smaller  than  the  rms  errors  due  to  the 
self  noise  of  the  gradiometer. 

There  are  two  remaining  sources  of  jitter  error: 
linear  jitter  sensitivity  at  sp i n - f requency  harmonics,  and 
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Estimated  Power  Spectrum  for  Pitch,  Roll,  or  Yaw 
Rate  Experienced  by  Airborne  Gradiometer  (No 
Plat  form  Jit  ter) 


Figure  3.6-2 


quadratic  sensitivity  to  stiction  jitter  in  the  inertial  plat¬ 
form.  The  rms  gradient  errors  caused  by  linear  sensitivity 
can  be  estimated  from  the  sensitivities  listed  in  Table  3.4-1. 
The  largest  sensitivity  listed  there  is  6.6x10“*  E/(rad/s), 
while  the  jitter  spectrum  in  Fig.  3.6-2  has  a  value  of  about 
3x10  ( rad/s )2/Hz  at  the  sp in- f requency  harmonics.  It  then 

follows  from  Eq.  3.4-9,  with  all  sensitivities  equal  to  6.6x10 

E/(rad/s)  and  all  spectral  densities  set  equal  to  3x10  (rad/ 

2 

s)  /Hz,  that  the  error  spectrum  in  each  of  the  two  gradiometer 
output  channels  is  bounded  for  |  f  I  <  0.024  Hz  by 

Sfr(f)  <  0.04  E2/Hz  (3.6-4) 

The  corresponding  rms  error  bound  is  0.04  E.  The 
conclusion  to  be  drawn  from  these  calculations  is  that  induced 
errors  due  to  linear  jitter  sensitivity  are  negligible  compared 
to  the  self  noise  of  the  gradiometer,  which  has  an  rms  value 
of  2  E. 

Errors  due  to  i nert i a  1 -pi at  form  jitter  are  a  poten¬ 
tially  significant  source  of  error  during  gravity  surveys. 

Bell  has  proposed  to  measure  the  platform  jitter  with  a  gyro 
and  then  continuously  compute  corrections  to  the  gradiometer 
outputs  in  real  time.  The  feasibility  ol  this  jitter  compen¬ 
sation  is  discussed  in  Appendix  A,  where  it  is  concluded  that 
Bell's  approach  to  compensation  is  viable,  but  the  magnitude 
of  the  residual  error  depends  on  a  number  ol  questions  that 
can  not  be  answered  with  the  presently  available  data. 

3.6.3  A  Bound  on  the  RMS  Error  Due  to 
Linear  Accel e raj  ion 

The  practical  significance  of  the  linear  acceleration 
spectrum  in  fig.  3.6-1  is  now  estimated  by  calculations  similar 


to  those  used  above  tor  jitter.  To  compute  an  upper-bound  for 

the  errors  due  to  linear  accelerations,  the  largest  sensitivity 

in  Table  3.3-1  (4.  4  E/mg)  is  used,  along  with  the  worst-case 

-  S  2 

acceleration  value  S  (0.25)  -  1x10  g  /Hz  in  Fig.  3.6-1;  the 
upper  bound  on  the  error  spectrum  Spp(f)  in  each  output  channel 
of  the  gradiometer  is  then  given  by  Eq .  3.3-14,  with  ail  spec¬ 
tral  densities  equal  to  S  (0.25)  and  all  sensitivities  equal 
to  A. 4  E/mg.  The  result  is,  for  |f|  <  0.02A  Hz, 

Srr(f)  <  5800  E2/Hz  (3.6-5) 

The  corresponding  rms  error  bound  is  17  E.  This  bound 
ing  value  is  eight  times  the  rms  error  due  to  the  grad iomet er ' s 
self  noise,  which  indicates  that  errors  due  to  linear  acceler¬ 
ations  may  be  a  significant  error  source  in  airborne  gravity 
surveys . 

3.7  ERROR  MODEL  FOR  AN  AIRBORNE  SURVEY  SYSTEM 
CONTAINING  THREE  BELL  ('.RADIOMETERS 

3.7.1  1 n troduct ion 

In  projected  gradiometric  gravity  surveys,  three  gradi 
ometers  are  to  be  used  simultaneously  to  recover  all  five  inde¬ 
pendent  components  of  the  gravity  gradient.  This  section  de¬ 
scribes  the  error  model  for  such  an  instrument  triad  consisting 
of  three  Bell  GGl's. 

3.7.2  Recovery  of  t_he_  Gravi  ty  Grad  i  cnt  f  rom 
Grad  ionic  t  er- System  Out  pu  t  s 

Each  of  the  three  grad i omet ers  comprising  l he  instru¬ 
ment  triad  has  two  outputs:  the  ini  inr-cli.inricl  gradient  me. is- 


urement  I  j(t)  and  the  e  ross- charine  1  measurement  f^(t).  These 


six  output  functions  are  arranged  in  a  vector 


f-o(t) 


r  j  ( t ) 


1  c< L ) 


f  j(U 


r(-(t ) 


i  j(t ) 


rc;(l } 


(3.7-1) 


where  the  superscripts  indicate  particular  gradiometers  in  th 


triad,  and  the  subscript  on  l()(t)  indicates  output  signals. 


11  the  gradiometer  system  were  free  of  errors,  then 
the  six  distinct  components  of  the  gravity  gradient 
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"uld  be  computed  I  loin  t  lie  system  output  b  (t  )  as  lollows 
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and,  e.g.,  r xz  =  3gx/3z  =  3gz/3x ,  with  gx  and  g?  being  the  x 
and  z  components  of  the  gravity  vector.  This  transformation 
is  derived  in  Appendix  C.2.  It  corresponds  to  an  x ,  y,  z  rec¬ 
tangular  coordinate  system  with  Lhe  x-axis  pointing  east,  the 
y  axis  pointing  north,  and  the  z  axis  pointing  up.  In  addition, 
the  first  gradiometer  in  the  instrument  triad  is  assumed  to 
have  a  vertical  spin  axis,  the1  second  gradiometer  has  a  south¬ 
pointing  axis,  while  the  third  instrument's  axis  points  east. 

Strictly  vertical  and  horizontal  spin  axes  are  assumed  because 
the  formal  test  data  on  the  baseline  CGI  were  obtained  with 
the  gradiometer  in  these  two  orientations. 

The  measurement  described  in  Eq .  3.7-3  also  applies 
to  gradiometer  triads  in  which  the  spin  axes  are  inclined  at 
5 70  to  the  vertical.  (Such  an  "umbrella"  configuration  is 
p.  ned  for  the  Gravity  Sensor  System  now  being  designed  by 
Bell.)  In  this  case  the  elements  of  T  in  Eq .  3.7-4  are  changed 
to  reflect  the  geometry  of  the  sensors.  Equation  3.7-3  also  ap¬ 
plies  to  gradiometer  triads  consisting  of  Draper  grad i ometers ; 
this  simply  requires  a  redefinition  ol  the  output  vector  E  (t) 
and  appropriate  choice  of  matrix  T  . 

T  is  the  transformation  that  maps  ideal  measurement 

O  r 

-i  to’-s  (error-tree  gradiometer  output  signals)  into  gradient 
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vectors  that  satisfy  Laplace’s  equation  for  the  gravity  poten¬ 
tial.  Since  T  is  a  known  matrix,  it  is  sufficient  to  develop 
an  error  model  for  the  gradioraeter  output  vector  [^(t).  The 
resulting  errors  induced  in  the  gradient  vector  T(t)  can  then 

x 

be  computed  from  Eq.  3.7-3. 

3.7.3  Gradiometer -System  Error  Model 

In  discussing  error  models  it  is  convenient  to  decom¬ 
pose  the  measurement  vector  r  (t)  into  two  component  vectors 

ro(t)  ;  +  V-e{t)  (3.7-5) 

where  Fj^(t)  is  an  error-free  i.deal  measurement  vector  and 
r  ( t )  is  a  vector  containing  the  six  error  signals  in  the  out¬ 
puts  of  the  gradiometer  triad.  The  error  model  for  simulating 
airborne  gradiometric  gravity  surveys  computes  the  spectral 

density  matrix  Sj  p(f)  of  these  errors.  The  structure  of  the 

e  e 

error  model  is  shown  in  Fig.  3.7-1,  where  the  errors  caused  by 
aircraft  vibrations  (linear  accelerations  and  angular  velocities) 
are  modeled  along  with  the  self  noise  of  the  gradiometers . 

3 . 7 . A  A ircraft  Linear  Accelerations 

In  Fig.  3.7-1  the  aircraft  linear  accelerations  at  the 
mounting  site  of  the1  gradiometer  are  characterized  by  a  vector 


;'It  is  recognized  that  Eq .  .3.7-3  does  not  necessarily  describe 
an  optimal  estimator  for  l(t)  when  the  measurements  are  noisy, 
but  rather  a  convenient  transformation  for  error-modeling  pur¬ 
poses.  References  23-25  address  the  estimation  problem. 
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Figure  3.7-1  Airborne  Gradiomet ri c  Survey  Error  Model 
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(3.7-6) 


whose  elements  are  the  power  spectra  of  random  accelerations 
in  the  z,  y,  and  x  directions.  This  ordering  corresponds  to 
up,  across,  and  along  track  axes  for  a  west-to-east  flight 
path.  In  the  error  model  developed  in  Section  3.5.2  these 
three  components  of  acceleration  are  independent  identical 
random  processes  with  the  spectral  density  SA.(1) 
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(3.7-7) 


where  ( f )  is  given  in  equation  3.5-1. 

ci  <i 

3.7.5  Aircraft  Angular  Velocities 

In  similar  fashion  the  aircraft  angular  rates  at  the 
gradiometer  mounting  site  are  represented  by  the  vector 


-«a * f  * 
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(3.7-8) 


whose  elements  are  the  power  spectral  densities  of  the  angular 
velocities  about  the  y,  x,  and  z  axes.  This  ordering  corre¬ 
sponds  respectively  to  pitch,  roll,  and  yaw  axes  for  a  west- 
to-east  flight  path.  These  three  components  of  velocity  are 
modeled  as  independent  indent  ical  random  processes  with  spec¬ 
tral  density  (  f ) 

-QQ ( f  * 

A 

where  S^q(0  is  given  in  Eq .  3.5-1. 

3.7.6  I  so  1  a  t i on  System  and  1 nert i al  PI  at  form 

The  output  of  the  vibration  isolation  and  local -level 
platform  models  in  Fig.  3.7-1  is  the  vibration  vector 


SQQ ( f  ) 
SQQ ( f  * 


(3.7-9) 
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which  contains  the  x,  y,  and  z  components  of  both  the  linear 
acceleration  and  angular  velocity  experienced  by  the  gradi- 
ometer  instrument  triad.  From  Eqs.  3.6-1  and  3.6-2  the 
vibration  vector  can  be  written  out  as 


§aa(f>' 

lTIS<12"fd2  §aa<f> 

_-QQ(f ) . 

|TpL(i2nf)  TIS( i2nf )  | 2  s£ft(f) 

where  the  transfer  functions  Tpp  and  TjS  of  the  inertial  platform 
and  passive  isolation  systems  are  given  in  Eqs.  3.5-2  and  3.5-3. 

3.7.7  Gradiometer-System  Vibration  Sensitivi ties 


The  linear  sensitivities  of  the  gradiometer  to  vibration 

are  represented  in  Fig.  3.7-1  by  a  6  x  6  mean-square  sensitivity 

matrix  S(f)  =  (S,.(f)],  whose  element  is  the  mean-square 

1  J  £  K 

sensitivity  of  the  i  output  of  the  instrument  triad  to  the 
jth  element  of  the  vibration  vector  Sy(l).  The  output  vector 
W  of  white-noise  spectral  densities  is  given  (in  accordance 
with  Eq .  3 . 3- 14  )  as 


(3.7-12) 


4 

w  =  2  £  S(fk)S  (f  ) 

k=0  k  v  k 

where  f.  =  kf  ,  for  k  =  0,1,..., 4  and  f  =  0.28  Hz,  are  the 

K  S  S 

spin- frequency  harmonics  of  the  gradiometer. 

In  the  coordinate  system  used  for  the  gradiometer 
triad,  the  mean-square  sensitivity  matrix  has  the  form 


S(fk)  =  [SLlN(fk)  !  SANG(fk)] 


where 
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2  2 

in  which  and  denote  the  squares  of  the  rms  linear 

sensitivities  defined  in  Sections  3.2.1  and  3.4.1,  HSA  and  VSA 
denote  horizontal  and  vertical  spin  axes,  and  the  symbols  "  ||  " 
and  "x"  respectively  denote  vibration  axes  that  are  parallel 
and  normal  to  the  spin  axes.  The  numerical  values  for  the  rms 
sensitivities  that  have  been  estimated  from  available  test  data 
and  that  are  used  in  the  TASC  error  model  are  listed  in  Tables 
3.3-1  and  3.4-1.  The  values  in  these  tables  are  expressed  in 
E/mg  or  10^  E/(rad/s)  and  must  therefore  be  squared  and  scaled 
properly  to  yield  ^^(f^)  in  E^/g^  — 1  ,c  '  " — J/-'^ 

as  shown  in  Appendix  C.3. 


and  in  E  /(rad/s)' 


The  diagonal  white-noise  spectral  density  matrix  D1AG(W) 
is  constructed  from  the  six  independent  white-noise  densities 
in  W: 


DIAG(W)  = 


Wj  0  0  0 
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0  0  W3  0 
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(3.7-16) 


W5  o 
o  w6 


This  diagonal  form  represents  the  independence  of  the  vibration- 
induced  white  noises  in  the  six  outputs  of  the  gradiometer 
triad . 


Discussion  of  Omitted  Terms  -  The  current  version  of 
the  TASC  error  model  for  airborne  gradiometric  surveys  omits 
measurement  errors  caused  by  platform  jitter.  Active  jitter 
compensation  is  planned  by  Bell.  Although  not  yet  demonstrated 
"on-line,"  preliminary  tests  indicate  that  such  compensation 
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may  be  feasible.  The  analysis  in  Appendix  A  supports  the  con¬ 
cept.  Of  course,  to  demonstrate  the  practicality  of  jitter 
compensation  requires  that  a  compensation  scheme  be  imple¬ 
mented  and  verified  by  laboratory  tests.  As  such  data  become 
available  they  will  be  evaluated.  If  the  magnitudes  of  the 
residual  jitter  errors  are  large  enough,  it  will  then  be  jus¬ 
tified  to  include  platform  jitter  in  the  TASC  error  model. 

The  analysis  of  quadratic  sensitivity  to  aircraft 
jitter  in  Section  3.6.2  showed  that  the  resulting  rms  measure¬ 
ment  errors  are  expected  to  be  more  than  six  orders  of  mag¬ 
nitude  smaller  than  the  self  noise  of  the  gradiometer.  For 
this  reason,  quadratic  sensitivity  to  aircraft  jitter  is  ex¬ 
cluded  from  the  error  model. 


3.7.8  Gradiometer  Self  Noise 


The  final  error  source  in  Fig.  3.7-1  is  the  gradi¬ 
ometer  system  self  noise,  which  is  represented  by  the  spectral 
density  matrix 
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(3.7-17) 


where  S^(f),  j  =  1,2,3,  is  the  spectral  density  matrix  for 
the  self  noise  in  the  two  outputs  of  the  j^  gradiometer  of 
the  instrument  triad.  The  self  noises  of  separate  gradiome- 
ters  are  independent  of  each  other.  The  first  gradiometer  has 
a  vertical  spin  axis,  and  its  spectral  density  matrix  is 


S 


1 

NN 


(  f  ) 


(.  VF.K 

*  1  1 


(  (  ) 


VHR 

1C 


<n 


eVFK 

‘Cl 


(  t  ) 


SVEK 

cc  '  1 


(3.7-18) 


88 


where,  from  the  analysis  in  Sections  2.4.1  and  2.4.2,  the 
trai  densities  for  the  inline  and  cross-channel  noises  are 
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,  0.01  mHz  <  |f  1  <  24  mHz 
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In  similar  manner,  the  two  horizontally  oriented 
dioineters  have  spectral  density  matrices 


SNN<f)  =  SNN(f) 
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All  componen  (  s  of  (he  gradiometer  system  error  model 
have  been  defined.  In  the  next  chapter,  this  model  is  com¬ 
bined  with  models  for  other  gravimetric  sensors  and  used  in 
simulations  to  demonstrate  some  typical  results  for  airborne 
gradiometric  surveys. 
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4. 


GRAD10METR1C  SURVEY  SIMULATION  RESULTS 


4.1  INTRODUCTION 

The  TASC  Multisensor  Survey  Simulation  software  uses 
models  for  the  gravity  field,  survey  geometries,  and  sensor 
errors  to  compute  the  rms  errors  inherent  in  statistically 
optimal  recoveries  of  gravity  quantities.  A  detailed  descrip¬ 
tion  of  the  Simulation  program  and  its  capabilities,  which  in¬ 
clude  trajectory  and  impact  error  analyses,  will  be  given  in  a 
forthcoming  report  (Ref.  25).  This  chapter  describes  typical 
simulation  results  obtained  for  airborne  gradiometric  surveys. 

A  variety  of  gradiometer  error  models  are  used  to  illustrate 
the  flexibility  of  the  software  and  to  indicate  the  usefulness 
of  gradiometer  error  models  for  determining  the  sensitivity  of 
recovery  errors  to  gradiometer  measurement  errors. 

The  overall  structure  of  the  Survey  Simulation  is 
shown  in  Fig.  4.1-1  for  the  case  in  which  airborne  gradiometry 
is  combined  with  either  satellite  altimetry  or  surface  gravi¬ 
metry  to  estimate  the  rms  error  in  point  estimates  of  the  grav¬ 
ity  disturbance  vector.  The  gradiometer-system  error  model 

JL 

described  in  Section  3.7  is  used  in  this  software  to  repre¬ 
sent  the  self  noise  and  sensitivities  to  mechanical  vibration 
of  the  Bell  GGI .  In  addition,  preliminary  error  models  for 
the  Bell  Ball-Bearing  GGI  and  the  Draper  Floated  Gravity 


*The  gradiometer  bandwidth  is  set  equal  to  50  mHz  to  model  a 
survey  sampling  rate  of  one  measurement  every  10  seconds. 
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Figure  A. 1-1  Structure  of  Multisensor  Survey  Simulation 

Gradiometer  are  also  included  for  comparison.  These  prelimin¬ 
ary  error  models  are  subject  to  revision  when  additional  test 
data  are  analyzed. 

A . 2  GRADIOMETER  SURVEY  GEOMETRY 

The  airborne  gradiometer  survey  model  covers  a  rec¬ 
tangular  area  as  shown  in  Fig.  A  .  2  -  1  .  The  survey  aircraft  has 
a  constant  altitude  of  6  km  and  flies  along  parallel  east-west 
tracks  spaced  10  km.  The  measurement  spacing  along  a  track  is 
l.b  km,  which  corresponds  to  one  measurement  every  10  seconds 
at  a  speed  of  TOO  kt  .  Each  track  is  2780  km  long  (5  hour  flight) 
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SUHVf  Y  GLOME  IKY  MODE  L 


•  SIKAIGHI  MNl  fAST  WES1  IKACKS 

■  1  HHEE  AXIS  INSTRUMENT  ATION 

•  I  RACK  SPACING  T  10  km 

■  ML  ASIIRI  MINT  SPACING  a  1  S  km 
110  sec  SAMPLING  at  300  ktl 

•  HEIGHT  OF  SURVEY  h  20.000  h 

•  LOW  FREQUENCY  DRIFT  CHARACTERISTICS 

•  WHITE  NOISE  FLOOR  AT  HIGH  FREQUENCIES 

■  SAMPLING  ERROR:  ALIASING  AND  OMISSION  AT  HIGH  FREQUENCIES 

■  FINITE  TRACK  LENGTH  WAVELENGTHS  LONGER  THAN  ?7R0  km 
(5  hr  flight)  ASSUMED  INACCURATE 

Figure  A. 2-1  Gradiometer  Survey  Geometry  and  Error  Modeling 

4.3  SURVEY  MODELS  FOR  OTHER  SENSORS 

The  satellite  altimetry  model  used  in  the  simulation 
is  chosen  to  represent  GEOS-3  survey  data.  The  track  geometry 
is  shown  in  Fig.  4.3-1,  together  with  numerical  values  for  the 
track  parameters.  The  error  model  for  the  altimetry  is  also 
described  in  this  figure. 

The  gravimeter  survey  model  is  described  in  Fig.  4.3-2. 
Measurements  are  located  on  a  grid  of  equally  spaced  points 
separated  from  each  other  by  150  nautical  miles  (nni)  and  ex¬ 
tending  sufficiently  far  beyond  the  grad i ome ter  survey  region 
to  provide  good  long  wavelength  recovery.  This  spacing  is 
adequate  to  recover  the  lower  frequency  components  of  the 
grav i t  y  field. 
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Figure  A . 3-1  Satellite  Altimeter  Survey  Geometry  and 
Error  Modeling  (GEOS -3  Parameters) 
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The  gravimeter  measurement  errors  are  uncorrelated  with  each 
other  and  have  an  rms  value  of  3  mgal  (this  includes  elevation 
errors  and  reduction  errors). 


4 . 4  GRAVITY  MODELS 

Two  gravity  models  are  used  to  simulate  both  "active" 
and  "moderate"  gravity  fields.  The  active- field  model  is  the 
sum  of  two  third-order  stationary  random  Markov  models  with 
very  strong  high-frequency  content;  the  rms  vertical  deflec¬ 
tions  are  16.8  sec,  and  the  rms  anomaly  is  113  mgal.  Such  a 
model  is  appropriate  for  areas  such  as  mountain  ranges,  ocean 
trenches,  escarpments,  and  other  areas  with  sharp  physical 
relief. 


The  moderate-field  model  is  relatively  low  in  high- 
frequency  energy  and  has  behavior  close  to  the  worldwide  aver¬ 
age.  It  is  a  5-shell  attenuated  white  noise  model  (Ref.  28) 
with  rms  deflections  of  the  vertical  of  6.8  sec  and  rms  anomaly 
of  43  mgal. 


4.5  GRADIOMETER  ERROR  MODELS 

The  gradiometer-system  error  model  for  the  Bell  base¬ 
line  GGI  is  described  in  Section  3.7.  The  error  model  for  the 
Bell  ball-bearing  GGI  is  similar  except  for  the  self-noise, 
which  is  modeled  to  reflect  the  preliminary  data  analysis  de¬ 
scribed  in  Section  2.5.  The  Draper  Floated  Gravity  Gradiometer 
No.  2  is  modeled  in  accordance  with  the  noise  spectrum  presented 
in  Section  2.6,  plus  a  preliminary  analysis  of  vibration  test  data. 
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4.6 


SIMULATION  RESULTS  FOR  THE  GOAL  GRADIOMETER  PLUS 
A  SATELLITE  ALTIMETER 


Before  presenting  simulation  results  for  the  Bell  and 
Draper  gradimeters,  it  is  instructive  to  consider  the  rms  recov¬ 
ery  errors  for  a  fictitious  "goal"  gradiometer ,  which  is  defined 
as  having  a  self  noise  with  a  flat  power  spectrum  and  an  rms 
value  of  1  eotvos  when  a  10  second  moving-window  average  is 
used  to  obtain  one  measurement  every  10  seconds.  The  goal 
gradiometer  is  defined  to  be  insensitive  to  all  environmental 
disturbances.  The  along-track  vertical  deflection  recovery 
errors  for  this  gradiometer  model  and  the  active  gravity  model 
are  described  in  Fig.  4.6-1,  which  shows  the  (two-sided)  power 
spectral  density  of  the  errors  as  a  function  of  wavelength. 

The  lowest  curve  corresponds  to  recovery  errors  when  both  gradi- 
ometry  and  satellite  altimetry  are  optimally  combined.  Also 
shown  are  the  recovery  errors  when  the  altimeter  data  and  the 
gradiometer  data  are  used  separately.  The  effect  of  limited 


Figure  4.6-1  Gravity  Field  Recovery  by  Wavelength 

(Active  Gravity  Model) 
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survey  area  on  long-wavelength  recovery  from  gradiometer  data 
alone  is  also  illustrated.  For  the  survey  simulation  consid¬ 
ered  here,  the  goal  gradiometer  data  provide  most  of  the  grav¬ 
ity  information  at  wavelengths  less  than  1000  nm,  while  the 
altimeter  data  supply  most  of  the  information  for  wavelengths 
longer  than  1000  nm. 


4.7  SIMULATION  RESULTS  FOR  GRAD10METRY  COMBINED  WITH 
GRAVIMETRY  OR  ALTIMETRY 


The  computed  rms  recovery  errors  for  several  gradi¬ 
ometer  surveys  scenarios  are  listed  in  Table  A.  7-1.  These 
results  were  obtained  with  the  smooth  worldwide  gravity  model. 


TABLE  A. 7-1 

MULTI SENSOR  GRADIOMETER  SURVEY  PERFORMANCE 
(WORLDWIDE  GRAVITY  MODEL-RMS  GRAVITY  DISTURBANCE 
FIELD  IS  6.8  SEC,  A3  MGAL)  AIRBORNE  GRADIOMETRY 
COMBINED  WITH  LAND  SURVEY/OCEAN  SURVEY 


H 


GRADIOMETER 

NOR  I  II 
DEFLECTION 
ERROR 

I  set  - rms ) 

EAST 

DEI  LECT  ION 
ERROR 
( sc<  -  rms ) 

VERTICAL 
DISTURBANCE 
I.KKOK 
(mj>,t  1  -  rms  ) 

NONE* 

5. A/1 .A 

T. A/1 . 5 

HH 

GOAL  GRADIOMETER 

0.6/0. 5 

0.5/0. 5 

38/3.3 

BELL  GAS  BEARING  MODEL 

■nMH 

BALL  BEARING  PROTOTYPE 1 

pnjfi 

1.5/1  2 

12.1/7.7 

CSDI.  FLOATED  GRADIOMETER  NO.  2** 

[HBQSi 

0.5/0. 3 

6. A/2. 9 

CSDL  WITH  CURRENT  NOISE  FLOOR  EX- 

0  5/0. A 

0  .  A/0. 3 

3. »/2  .  3 

TENDED  TO  LOW  FREOIIENCIES 

BALL  BEARING  PROTOTYPE  WITH  CURRENT 
NOISE  FLOOR  EXTENDED  TO  LOW 
FREQUENCIES! 

1  3/1.1 

,  -  — 

1.3/11 

8  8/7.5 

•i.o.,  SURFACE  GRAVIMETRY  AI.ONK  OR  CEOS  -  1  SATELLITE  Al.TIMF.TKY  AI.ONI. 
I  BASED  ON  PRF.I.IMINAHY  DATA 

•  BASED  ON  PRELIMINARY  ANALYSIS  OF  VIBRATION  IT  ST  DATA 
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The  rms  errors  are  listed  in  pairs  separated  by  slashes:  values 
to  the  left  of  the  slashes  are  the  errors  for  gradiometry  plus 
gravimetry ,  while  values  to  the  right  are  for  gradiometry  com¬ 
bined  with  satellite  altimetry . 

The  first  row  represents  the  limiting  case  of  no  gradi¬ 
ometry.  The  resulting  errors  are  large  because  the  gravimetry 
and  altimetry  lack  information  about  the  gravity  field  at  me¬ 
dium  and  short  wavelengths.  The  second  row  applies  to  the 
fictitious  goal  gradiometer  system  introduced  in  the  previous 
section.  The  third  row  corresponds  to  the  Bell  GGI  system 
error  model  that  was  developed  in  this  report,  while  the  fourth 
row  was  obtained  with  a  preliminary  error  model  for  the  Bell 
ball-bearing  gradiometer 

Results  for  tht  Draper  (CSDL)  Floated  Gradiometer  No. 

2  are  in  rows  5  and  6.  The  last  two  rows  (6  and  7)  contain 
estimated  recovery  errors  for  idealized  Draper  and  Bell  gradi- 
ometers  that  are  defined  to  be  free  of  long-term  random  drift, 
i.e.,  there  are  no  low-frequency  "red"  components  in  their 
self  noises.  This  limiting  case  may  be  approximated  in  prac¬ 
tice  by  intentionally  rotating  the  gradiometer  platform  so  as 
to  modulate  the  gradients  sensed  by  the  gradiometer.  Additional 
signal  processing  is  then  required  for  demodulating  the  gra¬ 
diometer  measurements  to  recover  the  gradient  components.  An 
analysis  of  intentional  platform  motion  is  given  in  Appendix  D. 

Table  A. 7-2  lists  projected  rms  recovery  errors  for 
the  same  gradiometer  surveys  covered  by  Table  A. 7-1,  except 
that  here  the  act ive  gravity  model  is  used.  The  much  larger 
high-frequency  energy  in  the  active  gravity  model  yields  larger 
rms  recovery  errors  as  compared  with  the  worldwide  gravity  model. 
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TABLE  4.7-2 

MULTISENSOR  GRAD IOMETER  SURVEY  PERFORMANCE  (ACTIVE  GRAVITY 
MODEL-RMS  DISTURBANCE  FIELD  IS  16.8  $EC,  113  MGAL)  AIRBORNE 
GRADIOMETRY  COMBINED  WITH  LAND  SURVEY/OCEAN  SURVEY 


GRADIOMETF.R 

NORTH 

DEFLECTION 

__ERROR 

(sec-rms) 

EAST 

DEFLECTION 

JCRROR 

(scc-rms)  j 

VERTICAL 
DISTURBANCE 
ERROR 
( mga  1  -  mis ) 

NONE* 

16.4/5. A 

| 

16.4/8.0 

110/46.0 

GOAL  GRAD IOMETER 

1.4/1. 3 

1.1/0. 9 

8. 0/7. 5 

BELL  GAS  BEARING  MODEL 

4. 0/2. 2 

1.9/1 .6 

20.9/12  7 

BALL  BEARING  PROTOTYPE^ 

4. 1/2. 5 

2. 3/2.0 

22.3/15.1 

CSDL  FLOATED  GRAD  IOMETER  NO.  2** 

3. 5/1. 9 

1.0/0. 8 

17  .4/9.6 

CSDL  WITH  CURRENT  NOISE  FLOOR  EX¬ 
TENDED  TO  LOW  FREQUENCIES 

1.4/1  .3 

0.9/0. 8 

8. 0/7.0 

BALL  BEARING  PROTOTYPE  WITH  CURRENT 
NOISE  FLOOR  EXTENDED  TO  LOW  FRE¬ 
QUENCIES! 

2. 1/2.0 

2. 1/2.0 

14.2/13.3 

*i.e.,  SURFACE  GRAVIMETRY  ALONF.  OR  GF.OS-3  SATELLITE  ALTIMETRY  ALONE 
t BASED  ON  PRELIMINARY  DATA 

**BASED  ON  PRELIMINARY  ANALYSIS  OF  VIBRATION  TEST  DATA 
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5. 


AIRBORNE  GRAD IOMETER  TEST  PLAN  ELEMENTS 


5.1  STATUS  OF  ONGOING  GRAD IOMETER  DEVELOPMENT 

The  Bell  Aerospace  Division  of  Textron,  Inc.  (under 
contract  with  the  Sperry  System  Management  Division  of  Sperry 
Rand  Corp.  and  sponsored  by  the  Strategic  Systems  Project  Of¬ 
fice  of  the  Navy)  is  currently  developing  an  Advanced  Develop¬ 
ment  Model  (ADM)  gravity  gradiometer  system,  which  is  to  include 
the  following: 


•  Three  ball-bearing  gradiometers 

•  Three  gyros  suited  for  carouselling  ex¬ 
periments  and  jitter  compensation 

•  Inertial  platform 

•  Temperature  controlled  binnacle 

•  vibration  isolators 

•  Associated  electronic  equipment. 

The  package  is  designated  the  Gravity  Sensor  System 
(GSS)  ADM.  Preliminary  to  the  ADM,  Bell  is  testing  two  ball¬ 
bearing  gradiometers,  designated  AF-2  S/N  1  and  S/N  2,  which 
are  development  versions  of  the  gradiometers  to  be  used  in  the 
GSS.  The  delivery  date  for  the  GSS  ADM  is  planned  for  October 
1981,  as  indicated  in  Fig.  5.1-1.  If  the  Navy  development 
proceeds  according  to  plan,  Beil  will  have  the  capability  to 
produce  a  gradiometer  triad  for  airborne  use  soon  thereafter. 


100 


R-JOMt 


I -  GSS  ADM  DEVELOPMENT - ^H| -  AT  SEA  TESTING 

(Already  Scheduled) 


Jan  1900 


_J _ 

Jan  1961 


_l _ 

Jan  1902 


- 1 _ 

Jan  1903 


Figure  5.1-1  Ongoing  GSS  ADM  Development 

Although  actual  flight  testing  of  gradiometers  is 
several  years  away,  and  the  concomitant  details  of  the  test 
scenario  are  unknown,  several  factors  and  issues  are  now  clear. 
This  section  briefly  discusses  these  areas  and  directs  emphasis 
toward  the  "long  lead-time"  items.  The  overall  objective  is 
to  illustrate,  for  planning  use  by  DMA,  the  key  portions  of  a 
coordinated  airborne  gradiometer  test  program.  An  overall 
schedule  which  such  a  program  might  follow  is  presented  in 
Fig.  5.1-2. 

5.2  TEST  PROGRAM  ITEMS 

In  one  sense  the  initial  portion  of  an  airborne  grad- 
iometer  test  program  has  already  begun  with  the  analysis  of 
both  Bell  and  Draper  moving-base  data.  While  a  continuing 
study  of  laboratory  data  is  a  necessary  adjunct  to  a  success¬ 
ful  airborne  gradiometer  program,  this  section  looks  beyond 
the  expected  evolutionary  developments  of  laboratory  testing. 
The  items  below  are  directed  towards  the  epoch  following  de¬ 
livery  of  a  gradiometer  triad  and  its  stabilized  platform  to  a 
test  site  remote  from  the  manufacturer ' s  laboratory. 
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Figure  5.1-2  Elements  of  an  Airborne  Gradiometer  Test  Program 


Ground  Truth  -  The  accuracy,  density,  and  total 
coverage  of  gravity  data  required  to  support  air¬ 
borne  testing  is  yet  to  be  determined.  Ancillary 
issues  include  types  of  data  most  appropriate 
for  the  database  (e.g.,  gravity  vs.  deflections 
of  the  vertical),  the  data  subsets  required  for 
moving-base  checkout  testing  on  the  surface, 
data  needed  for  system  calibration,  and  the  pro¬ 
cessing  techniques  to  be  used  for  upward  continu¬ 
ation  to  aircraft  altitude.  The  acquisition  and 
reduction  of  these  data  will  require  considerable 
lead  time. 

Inertial  Navigat ion  System  (INS)  and  Navigation 
Aid  SelecfTon  -~The  choice  of  aircraft,  test 
si te ,  and  test -program  navigation  accuracy  re¬ 
quirements  will  be  the  chief  constraints  in  se¬ 
lecting  an  INS.  The  use  of  Navigation  Aides 
(Navaid)  will  largely  be  dictated  by  availabil¬ 
ity  and  cost.  Integration  of  a  Global  Positioning 
System  (GPS)  receiver  into  the  system  is  particu¬ 
larly  promising.  In  the  absence  of  GPS,  a  radar 
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altimeter  and  doppler  radar  would  probably  be 
required.  Lead  and  back-up  navigation  systems 
need  to  be  defined,  selected,  and  procured  in 
pace  with  gradiometer  acquisition. 

Data  Recording  and  Monitoring  -  The  equipment 
and  system  required  for  adequate  monitoring  of 
the  gradiometer  triad's  performance  must  be  de¬ 
fined.  Among  the  issues  to  be  resolved  include 
which  variables  are  to  be  recorded,  necessary 
environmental  data  monitors,  the  degree  of  post¬ 
mission  vs.  real-time  data  processing  to  be  done, 
onboard  storage  vs.  telemetry  as  well  as  data 
formats  and  hierarchy.  In  addition  to  allowances 
for  system  procurement  and  setup  time,  a  period 
of  test  equipment  software  development  is  likely 
to  be  needed. 

Ground  Testing  -  Prior  to  operation  in  an  air- 
craft,  testing  and  calibration  of  the  gradiometer 
triad  and  platform  must  be  performed  at  the  air¬ 
field.  The  detail  and  extent  of  this  testing 
will  be  dictated  by  the  philosophy  of  those  di¬ 
recting  the  program  at  the  time.  Ground  testing 
can  range  from  simple  check-out  of  the  gradiometer 
instruments  and  recording  systems  prior  to  flight, 
to.  a  complete  ground  test  program  designed  to 
provide  high  predictability  of  flight  test  per¬ 
formance.  Regardless  of  the  specific  program 
chosen,  well-defined  procedures  and  schedules 
must  be  developed.  Similarly,  the  analysis  ob¬ 
jectives  for  the  ground  test  data  should  be  speci¬ 
fied  early  in  the  program.  Some  software  develop¬ 
ment  is  likely  to  be  required  for  the  data  analysis 

Flight  Test  Operational  Plan  -  Because  of  the 
high  cost  of  operating  test  aircraft,  it  is  de¬ 
sirable  to  maximize  the  information  return  from 
each  flight.  A  pre-flight  checkout  procedure  is 
therefore  necessary  for  the  gradiometer/INS/ 

Navaid  systems,  as  well  as  data  recording  and/ 
or  telemetering  equipment.  In-flight  procedures 
are  required  to  verify  proper  data  acquisition 
and  to  make  data- reasonabl eness  checks.  Post¬ 
flight  calibration  procedures  should  also  be 
well -defined . 

Flight  Test  Data  Ana  1 y sis/Evaluation  -  Compa r i - 
son  of  i  n  -  fl  i  gh  f~g  rati  1 ome  t  e  r  measurements  to 
data  derived  from  ground  truth  will  most  likely 
be  performed  after  the-  flight.  Many  other  facets 
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of  gradiometer  performance  will  probably  have  to 
be  studied  post-mission  too.  Sufficient  time 
and  resources  should  be  allocated  for  this  ac¬ 
tivity.  Experience  has  shown  that  unanticipated 
problems  can  seriously  impact  the  effort  required 
in  this  phase  of  the  test  program.  Close  coor¬ 
dination  of  the  gradiometer  navigation,  data 
acquisition,  and  truth  data  development  are  re¬ 
quired  to  prevent  such  surprises. 

Technical  Improvement/Iteration  Procedures  -  In 
a  development  program  such  as  gradiometry,  it  is 
in  the  interest  of  all  parties  to  assure  that 
needed  adjustments,  design  changes,  and  repairs 
are  affected  from  time  to  time.  It  is  also  nec¬ 
essary  that  such  procedures  occur  with  the  cog¬ 
nizance  of  appropriate  personnel  and  in  an 
orderly  fashion.  Procedures  for  such  actions 
must  be  designed  which  assure  adequate  documen¬ 
tation  and  discipline  for  all  with  access  to  the 
instrumentation,  but  encourage  creativity  and 
rapid  response  to  problems  that  may  develop. 

Test _ Program  to  Survey  Transition  -  Once  confi¬ 

dence  in  Fhe  gradiometer  is  established,  and 
error  estimates  for  its  airborne  performance  are 
obtained,  flights  might  be  extended  to  areas 
where  ground  truth  data  are  not  dense.  Simi¬ 
larly,  the  training  of  new  nonf 1 ight-test  per¬ 
sonnel  might  be  appropriate.  Such  flights  could 
be  interspersed  with  test  program  flights  as  the 
program  nears  completion.  To  effect  an  efficient 
transition  to  survey  procedures,  requirements 
and  goals  of  the  survey  program  should  be  defined. 


5.3  APPROPRIATE  NEXT  STEPS 

The  current  umbrella-angle  design  for  a  triad  of  Bell 
grad  ioineters  (Ref.  3)  is  suitable  for  airborne  survey  opera¬ 
tions  with  a  moderate  augmentation  of  the  instrument's  band¬ 
width.  This  bandwidth  increase  is  not  expected  to  be  difficult. 

A  gradiometer  inertial  platform,  integral  with  the 
gradiometer  triad,  is  under  development  by  Bell.  However,  it 
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is  currently  unclear  to  what  extent  the  Bel ] -designed  platform, 
intended  for  marine  use,  will  have  to  be  modified  for  the  air¬ 
borne  application.  This  is  the  chief  critical-path  item  for 
airborne  gradiometry  as  of  this  date. 

Bell  has  not  yet  implemented  their  gradiometer ’ s  jitter 
compensation  scheme  in  an  on-line  environment.  Although  the 
concept  appears  sound  (see  Appendix  A),  it  should  be  tested 
prior  to  an  extensive  commitment  to  an  airborne  test  program. 
Definitive  testing  would  demonstrate  an  operative  gradiometer 
mounted  on  the  intended  aircraft-compatible  inertial  platform 
while  subjected  to  linear  and  angular  vibration  that  simulates 
the  host  aircraft.  The  testing  should  include  conditions  of 
angular  motion  (such  as  encountered  in  an  aircraft)  of  the 
stabilizing  platform's  base. 

As  stated  earlier,  continued  monitoring  of  the  manu¬ 
facturer's  development  and  laboratory  testing  is  necessary. 

This  effort,  in  addition  to  enabling  DMA  to  maintain  a  high 
level  of  technical  awareness  in  the  field  of  gradiometry,  also 
provides  the  flexibility  to  embark,  without  delay,  on  a  wider 
program  of  airborne  gradiometry  as  the  resources  become  available. 
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6. 


CONCLUSIONS  AND  RECOMMENDATIONS 


The  accuracy  with  which  gravity  can  be  estimated  from 
an  airborne  gradiometric  survey  depends  on  the  gravity  field, 
the  survey  geometry,  the  additional  data  from  other  sensors, 
and  the  gradiometer  measurement  errors.  This  report  has  fo¬ 
cused  on  the  last  of  these  factors  by  describing  error  models 
for  airborne  gradiometers .  These  models  have  been  incorpo¬ 
rated  in  the  TASC  Multisensor  Survey  Simulation  software, 
which  provides  estimates  of  recovery  errors  when  different 
types  of  gravity  survey  data  are  optimally  combined. 

The  error  models  in  this  report  were  chosen  to  fit 
test  data  for  the  baseline  Gravity  Gradiometer  Instrument  No. 

1  t,GGI  No.  1)  developed  by  the  Bell  Aerospace  Division  of 
Textron,  Inc.  This  is  the  gradiometer  for  which  comprehensive 
formal  test  data  have  been  analyzed,  and  its  performance  is  be¬ 
lieved  to  be  indicative  of  the  Bell  bail-bearing  gradiometers 
that  are  currently  being  built  for  field  tests  in  1981.  It  is 
recommended  that  the  error  models  in  this  report  be  updated  to 
fit  test  results  for  the  Ball-Bearing  instrument  when  additional 
data  for  that  gradiometer  are  available. 

Investigative  test  results  and  a  limited  amount  of 
formal  test  data  were  available  on  the  Floated  Gravity  Gradio¬ 
meters  (FGG's)  developed  at  the  Charles  Stark  Draper  Laboratory. 
Analysis  of  these  data  yielded  self-noise  spectra  for  the  FGG 
NO.  2  that  are  included  in  this  report  .  The  FGG  models  for 
v i brat  ion- i nduced  errors,  however,  are  only  preliminary  and 
are  subject  to  change  after  additional  data  are  analyzed. 
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Gradiometer  measurement  errors  are  due  to  systematic 
effects,  self  noise  originating  inside  the  instrument,  and  to 
environmental  disturbances  to  which  the  gradiometer  is  sensi¬ 
tive.  The  dominant  measurement  errors  during  an  airborne  survey 
are  expected  to  be  caused  by  the  self  noise  and  the  mechanical 
vibration  environment  of  the  gradiometers .  Major  emphasis  in 
this  report  was  therefore  placed  on  these  two  error  sources. 


Analysis  of  the  test  data  indicates  that  the  self 
noise  in  the  inline  and  cross  channels  of  the  Bell  baseline 
GGI  is  well  modeled  as  two  random  processes.  Power  spectral 
densities  for  these  processes  were  estimated  from  data  for 
both  vertical  (VSA)  and  horizontal  spin-axis  (HSA)  orienta¬ 
tions  of  the  gradiometer. 


•  At  frequencies  below  0.1  millihertz  (mHz), 
the  spectral  densities  are  inversely  pro¬ 
portional  to  the  frequency  squared.  This 
is  termed  "red"  noise,  and  it  accounts 
for  long-term  drift  in  the  gradiometer 
outputs 

•  At  higher  frequencies  (2  mHz  <  f  <  2U  mHz), 
the  spectral  densities  are  independent 

of  frequency.  This  is  the  "white-noise 
floor,"  and  it  accounts  for  short-term 
variability  in  the  gradiometer  outputs. 

•  The  rms  accumulated  random  drift  projected 
for  a  5-hour  survey  is  4.2  eotvos  (E)  for 
VSA  and  5.3  E  for  HSA.  The  measurement 
averaging  time  for  these  projections  is 

10  seconds. 


Analysis  of  test  data  for  the  Draper  FGG  No.  2  indi¬ 
cates  that  its  self  raise  has  the  following  properties: 


•  The  power  spectra  of  the  FGG's  self  noise 
are  similar  in  shape  but  sma  tier  i n  mag- 
nitude  as  compared  with  those  of  the  Bell 
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CGI  (the  EGG  No.  2  spectral  magnitudes 
are  typically  3  percent  of  the  corre¬ 
sponding  Bell  baseline  GG1  spectral  mag- 
ni tudes ) 

•  The.  projected  rms  accumulated  random 
drift  during  a  5-hour  survey  is  0.6  E 
for  a  10-second  averaging  time. 


Analysis  of  Bell  vibration  data,  together  with  flight 
test  data  on  the  vibrations  of  an  NKC-135A  aircraft,  led  to 
the  following  conclusions  regarding  vibration-induced  measure¬ 
ment  errors  during  airborne  gradiometric  surveys: 


•  Aircraft  angular  jitter  is  expected  to 
induce  small  measurement  errors.  In  con¬ 
trast,  uncompensated  inertial  platform 
jitter  Is  idenTITied  as  a  potentially 
important  error  source.  Compensation 
for  platform  jitter  appears  to  be  feasi¬ 
ble  and  may  reduce  the  j i t ter- induced 
errors.  The  s ucc e ssful  implementation 

of  such  compensaTTon , "however ,  remains 
to  be  demonst rated 

•  Aircraft  linear  vibrations  are  expected 
to  induce  uncorrelated  measurement  errors 
that  are  comparable  in  order  of  magnitude 
to  the  gradiometer ' s  white-noise  floor. 
This  conclusion  is  based  on  available 
flight  test  data. 


Error  models  were  developed  for  airborne  gradiometric 
survey  systems  that  consist  of  three  gradiometers  mounted  on 
one  inertial  platform.  These  error  models  were  incorporated 
into  the  Multisensor  Survey  Simulation  routines  and  used  to 
compute  the  rms  errors  inherent  in  statistically  optimal  re¬ 
coveries  of  point  values  for  the  gravity  vector  when  airborne 
grad  i  nine  t  ry  is  combined  with  long-wavelength  data  from  satel¬ 
lite  altimetry  or  land-based  gravimetry.  The  following  is  a 
sample  of  the1  simul.it  ion  results. 
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•  When  airborne  gradiometry  is  combined 
with  satellite  altimetry,  the  worldwide 
rms  vertical  disturbance  error  is  3.3 
mgal,  and  the  rms  along-track  vertical 
deflection  error  is  1.0  sec  for  the  Bell 
baseline  GGI 

•  For  the  Draper  FGG  No.  2  the  correspond¬ 
ing  rms  vertical  disturbance  error  is  2.9 
mgal,  and  the  rms  along-track  vertical 
deflection  error  is  0.3  sec. 


Future  modifications  of  the  gradiometer  error  models 
are  expected  to  include: 


•  Adjustments  of  parameter  values  to  re¬ 
flect  additional  test  data  from  the  Bell 
Ball-Bearing  GGI 

•  The  addition  of  a  platform  jitter  model 
if  justified  by  future  test  data 

•  Refinements  in  the  aircraft  vibration 
models  if  additional  flight  test  data 
representative  of  survey  conditions  be¬ 
come  available. 


From  the  data  that  were  analyzed  for  this  report  as 
well  as  the  technical  maturity  of  both  the  Bell  and  Draper  in¬ 
struments,  it  appears  likely  that  the  transition  from  labora¬ 
tory  to  various  types  of  field  tests  is  near.  In  this  context 
the  following  recommendations  are  made  tc  assist  in  the  plan¬ 
ning  and  development  of  an  airborne  gradiometric  survey  capa¬ 
bility.  The  focus  is  on  elements  that  should  be  addressed 
soon  to  provide  maximum  flexibility  in  later  airborne  gradi¬ 
ometer  program  activity. 


•  Current  gradiometer  data  analysis  studies 
should  be  continued  to  provide  DMA  with 
ongoing  high-level  technical  insight  into 


gradiometer  development.  Such  studies 
should  include  test  data  processing  as 
well  as  independent  mechanization  studies 
and  should  relate  gradiometer  design 
developments  to  the  airborne  survey 
scenario 

•  Independent  analysis  of  gradiometer  plat¬ 
form  requirements  should  be  initiated 

•  Planning  for  a  gradient  "truth  data  base" 
should  commence  at  the  earliest  possible 
date.  Even  before  the  selection  of  a 
test  site,  the  type,  accuracy,  densifi- 
cation,  and  extent  of  required  surface 
data  must  be  determined.  The  generation 
of  gravity  gradients  for  comparison  with 
gradiometer  outputs  at  aircraft  altitudes 
is  likely  to  require  considerable  time. 

The  necessarily  large  quantities  of  data 
should  therefore  be  matched  with  efficient 
algorithms  (e.g.,  GEOFAST  in  references 

32  and  33). 
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APPENDIX  A 


JITTER  COMPENSATION  TECHNIQUE 
FQR  THE  TSELL  C.RaDJ  OHETEE 


The  purpose  of  this  Appendix  is  to  examine  Bell's 
method  of  gradiometer  compensation  for  angular  jitter  motion 
of  the  inertial  stabilization  platform.  In  the  approach  stud¬ 
ied,  a  gyro  is  used  to  measure  the  jitter  and  to  compute  cor¬ 
rections  continuously  for  the  gradiometer  outputs.  The  com¬ 
pensation  is  to  be  applied  in  real  time. 

The  analysis  herein  finds  this  approach  to  be  viable 
if  the  gradiometer ' s  mechanical  structure  meets  certain  yet- 
to-be-demonstrated  goals.  A  description  of  these  requirements 
is  provided  in  Section  A. 2.  Conclusions  are  given  in  Section 
A. 3. 

A . 1  INTRODUCTION 

Concerns  have  been  expressed  in  the  gradiometer  com¬ 
munity  that  gyro-based  jitter  compensation  could  prove  inef¬ 
fective  because  of  the  wide-bandwidth  multiplications  required 
to  compute  the  jitter-correction  gradient  terms.  The  potential 
problem  is  that  gyro  noise  over  the  bandwidth  from  DC  to  the 
mechanical  response  cutoff  frequency  of  the  gradiometer  plat¬ 
form  would  be  excessive.  The  following  example  is  illustra- 
t  ive . 


Suppose  a  gradiometer  i mounted  on  a  "state-of-the- 
art”  platform  which  exhibits  randomly  occuring  limit -cycle 


behavior  with  an  average  frequency  of  20  Hz.  Consider  sinu¬ 
soidal  platform  motion  with  an  amplitude  of  1  arcsecond  (sec), 
i  .e.  , 


0  =  0  cos  2nft  (A. 1-1) 

m 

where  6^  =  1.0  sec  and  f  =  20  Hz.  The  angular  velocity  of  the 
limit  cycle  is 


8  =  -2nf0  sin  2nft 
m 


(A. 1-2) 


with  a  peak  value  ot 


0  =  2nf0  =  40 n  sec/s  (A. 1-3) 

m  m 

The  corresponding  centripetal  gradient  is 

f  =  6^  =  ^  -  I  0^  cos  Anft  (A. 1-4) 

2  m  2  m 


Only  the  DC  part  of  Eq .  A. 1-4  is  of  concern  because  the  double¬ 
frequency  term  is  beyond  the  grad iometer ' s  operating  bandwidth. 
However,  since  the  platform  jitter  is  not  stationary  (i,e.,  it 
comes  and  goes),  the  "DC"  term  in  Eq .  A. 1-4  appears  as  a  ran¬ 
dom  gradiometer  output. 


Using  the  jitter-rate  amplitude  of  (A. 1-3)  yields  the 
reel i 1 ied  term 


I  1)(.  =  8.0  •  1 02/r 2  sec =  1 86  E 


(A. 1-5) 


This  is  the  measurement  error  that  an  uncompensated 
would  exhibit  in  t  hr*  presence  of  platform  jitter  as 
above.  Note  that  if  the  limit -cycle  amplitude  were 
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gradiometer 
desc  ri bed 
a  factor 

f 

1 


t 


f 


of  ten  larger  (10  arcseconds),  the  effect  on  gradient  error 
would  be  a  factor  of  one  hundred  worse. 


Since  this  single- frequency  error  is  an  impulse  in 
the  frequency  domain,  it  can  be  observed  with  a  very  large 
signal- to-noise  ratio.  Therefore,  compensation  for  jitter- 
induced  error  presents  no  problem  in  the  single- frequency  case. 

However,  several  questions  remain: 


•  What  happens  if  the  limit  cycle  is  more 
realistically  characterized  by  a  spectrum 
containing  harmonics? 

•  What  is  the  net  effect  of  gyro  self-noise 
on  gradiometer  output  error? 

•  What  level  of  gyro  resolution  is  required? 

•  The  jit  ter- induced  gradient  error  origin¬ 
ates  as  an  acceleration  error.  Thus  the 
rectified  portion  of  the  jitter  is  com¬ 
prised  of  all  jitter  frequencies  contained 
within  the  platform's  mechanical  response 
bandwidth.  (Neither  the  gradiometer 
signal  bandwidth  nor  the  accelerometer 
constrainraent  loop  bandwidth  are  impor¬ 
tant.)  How  much  difficulty  do  high  fre¬ 
quency  mechanical  resonances  cause? 


A. 2  JITTER  COMPENSATION  ISSUES 

A. 2.1  Rectification  Effect 

Angular  motion  of  the  gradiometer  instrument  affects 
the  individual  sensed  gradients  via  the  following  relation 
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sr  = 


(»2  +  o2) 

ftx°y 

U)Q_. 


Vy 

(U)2  +  fi2) 
U>Q.. 


U)fi. 


wfl. 


-(O2  *  !12> 


(A. 2-1) 


where  6T  is  the  j i tter- induced  gradient  tensor  measurement 

error,  u>  is  the  jitter  rate  about  the  spin  axis,  and  Q  and  ft 

a  y 

are  jitter  rates  about  two  other  orthogonal  axes. 


Jitter  about  the  spin  axis  does  not  appear  in  the 
Bell  gradiometer ' s  output.  Thus,  terms  of  interest  are  gradi¬ 
ent  errors  of  the  form 

6r  =  Q2,  Q2,  QxQy  (A. 2-2) 

Examination  of  a  single  component  of  jitter  is  instructive. 
Consider  Q  to  be  Gaussian  random  process  and  to  have  the  power 

A 

spectrum  S  (f).  The  gradient  error  spectrum  of  6y  is  given 
by  (Ref.  29  Eq .  19-20) 

Srr(f)  =  2Sxx(f)  *  Sxx(f)  +  6(f)  (A. 2-3) 

2 

where  a  is  the  variance  of  fl  ,  6(f)  is  the  Dirac  delta  func- 

X  X 

tion,  and  *  denotes  convolution.  The  second  term  of  Eq .  A. 2-3 
is  the  source  of  the  rectification.  The  magnitude  of  this 
effect,  if  uncompensated ,  was  illustrated  in  Section  A.l. 


A . 2 . 2  Compensa  t i on 


error 


i. 

ft 


The  jitter  signal  D  is 
(also  a  Gaussian  random 


to  be  measured 
process ) 


by 


a  gyro  with 


119 


m  =  +  £g  ( A . 2-4 ) 

The  estimated  jitter  gradient  is 

6r  =  m2  =  qI  +  2EgQx  +  z2  (A. 2-5) 

The  actual  jitter  gradient  is 

6r  =  ft*  (A. 2-6) 

The  resulting  compensated  gradiometer  output  error  is 

6r£  =  6r  -  6r  =  2£g0x  +  z2  (a. 2-7) 

The  autocorrelation  of  6Tc  is  given  by 

R£|t)  =  4E  { t  g  <  t )  Qx(t)  £g(t-t)  ftx(t-x)} 

+  2E  { t  g ( t )  0x(t)  eg( t-T  )  eg(t-t)} 

+  2E  U  <t)  tit)  t( t-x)  ft  (t-T ) } 

&  S  6  A 

+  E  Ug(t)  £g(t)  eg(t-T>  EgU-T)l  (A. 2-8) 

where  t  is  signal  time  and  t  is  shifted  time.  E{ . . . }  is  the 
ensemble  expectation  operator.  Using  Eq .  4.9  of  Ref.  31,  Eq . 
A. 2-8  can  be  reduced  to 


R«(o  = 4  v°>  *  v>>  Rxx(I>  *  Rgg<’> 

♦ 2  V#»  V°>  +  V>  V>  4  V"  VT) 
* 2  V0)  V°>  +  VT>  V(T)  *  V"  V” 

♦  R2  (0)  +  2R2  (i) 

88  88  (A. 2-9) 
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where  the  correlators  functions  are  defined  according  to 

R  <t)  =  E  { £  ( t )  0  (t-t)}  (A. 2-10) 

gX  g  X 

Since  platform  jitter  and  gyro  errors  are  uncorrelated, 


=  V(I>  =  0 


(A. 2-11) 


Using  Eq.  A. 2-11  to  simplify  Eq .  A. 2-9  gives 


R  (t)  =  4  R  (t)  R  (t)  +  2R2  (t)  +  o  2  o  2 

££  gg  XX  gg  g  X 


(A. 2-12) 


Taking  the  Fourier  transform  of  Eq .  A. 2-12  gives  the  power 
spectral  density  (PSD)  of  the  residual  gradiometer  error  due 
to  jitter, 


S 

€  C 


(f) 


4  S 


gg 


(f) 


S 


xx 


<f>  *  2Sgg(f) 


s 

gg 


(f)  +  agA6(f) 

(A. 2-13) 


The  final  term  of  Eq .  A. 2-13  contributes  only  to  the  zero- 
frequency  component  of  6f‘  ;  the  rms  contribution  is 

6rDC  =  o2  (A. 2-14) 

£  g 

It  is  convenient  to  preserve  the  distinction  between 
the  DC  and  all  other  spectral  components  of  6ye ■  The  DC  effect 
is  treated  first. 


Gyro  errors  for  second  generation  dry-tuned  gyros 
such  as  those  used  the  LN-lb  inertial  navigation  system  have 

-  _  ^ 
"Gyro  vendor  data  indicate  that  error  spectra  as  low  as  10 

2 

(deg/hr)  /Hz  over  this  bandwidth  are  available. 
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1 


been  modeled  by  a  first-order  Markov  process  with  an  rms  value 

_3 

of  10  deg/hr  and  a  correlation  time  of  2  hours  (Ref.  30). 

Bell  has  tested  and  expects  to  use  an  advanced  version  of  this, 
or  a  similar  gyro,  for  jitter  detection.  The  gyro  error  PSD 
(white)  at  low  frequencies  is 

S  =  3.6  x  10'3(sic/s)2/Hz*  (A. 2-15) 

o 

It  is  instructive  to  consider  the  case  in  which  the 

JL 

increased  bandwidth  required  for  jitter  compensation  (100  Hz)” 
is  accomplished  only  at  the  expense  of  accepting  the  error 
spectral  density,  S  ,  over  the  entire  gyro  pass  band.  The 

o 

gyro  error  variance  is  then 
100 

a2  =  f  S  df  =  0.72  (s'Tc/s)2  (A. 2-16) 

o  J  to 

-100 

which  yields  a  DC  gradient  error  of 

6r^C  =  0.02  E  (A. 2-17) 

This  result  indicates  that  the  DC  component  is  not  a  signifi¬ 
cant  source  of  error. 


Turning  to  the  frequency  spectrum  of  61"^  ,  evaluation 

of  the  first  two  terms  of  Eq.  2.2-10  is  required.  Taking  the 

gyro  error  spectrum  to  be  white  out  to  its  cutoff  frequency 

f  ^  ( f . .  =  100  Hz)  gives 
co  co 


(f) 


_IS 


I  f  I  <  f 
I  f  I  >  f 


co 

co 


(A. 2-18) 


*100  Hz  is  used  as  the  effective  cutoff  frequency  for  platform 
mechanical  response. 
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which  is  to  he  inserted  into 


S  (f)  =  4S  (f)  *  S  (f)  +  2S  (f)*S  (f)  (A. 2-19) 

ee  gg  xxv  gg'  ggv 

It  is  convenient  to  rewrite  Eq.  A. 2-16  to  include  f  as  a 
parameter 


S  = 


g  2f . 


(A. 2-20) 


co 


For  the  numerical  values 
volution  in  Eq.  A. 2-19  yields  the 
in  Fig.  A. 2-1.  The  maximum  value 


considered,  the  second  con- 

( 2 ) 

error  spectrum  S'  '(f)  shown 

O  C 

is 


S(2)(0) 

£t 


=  =  3 

co 


10~A  E2/Hz 


(A. 2-21) 


This  result  is  more  than  5  orders  of  magnitude  smaller  than 
the  self-noise  spectral  density  of  the  Bell  GGI  discussed  in 
Chapter  2.  Therefore,  the  second  term  in  Eq .  A. 2-19  is 
negl igi ble . 


Figure  A . 2 -  1 


Frequency  Spectrum 


(2) 


£  £ 


(f) 
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It  remains  to  evaluate  the 
Eq.  A. 2-19.  Although  discussion  of 
until  Section  A. 2.4  it  is  useful  to 
jitter  model  as  illustrated  in  Fig. 


first  convolution  term  of 
jitter  spectra  is  deferred 
anticipate  a  worst-case 
A. 2-2. 


R— 51366 


Sxx  (sec/s)2/Hz 


1 

r44 

f  (Hz) 


-100 


-10  10 


100 


Figure  A. 2-2  Band-Limited  White-Noise  Jitter 
Spectrum  (Schematic) 


The  result  of  convolving  the  jitter  and  gyro  error 
spectra  is 


S<J><0)  =  0.06  E2/Hz  (A. 2-22) 

Although  this  spectral  level  is  somewhat  larger  than  the  one 
in  Eq.  A. 2-21,  it  is  still  tor  small  to  contribute  significant 
gradiometer  error. 

A. 2. 3  Resolution 

Suppose  the  gyro's  resolution  is  no  worse  than  the 
rms  value  of  its  broadband  noise,  given  in  Eq .  A. 2-16  as  0.85 
sec/s.  If  the  resolution  is  the  same  at  all  frequencies,  the 
worst-case  effect  occurs  when  the  jitter  is  wideband  and  of 
density  greater  than  the  resolution  density  over  the  entire 
band.  In  that  case  the  contribution  of  gyro  resolution  is 


119 


nearly  same  as  that  of  gyro  self-noise  and  is  therefore  negli¬ 
gible  in  comparison  to  the  gradiometer ' s  self  noise. 


A . 2 . 4  Gyro- to-Platform  Relative  Motion 

The  gradiometer  body  exhibits  resonances  that  occur 
near  the  mechanical  frequency- response  limits  of  the  structure. 
Because  of  the  modal  motion  associated  with  resonance  condi¬ 
tions,  the  gyro  is  subjected  to  different  angular  vibrations 
than  are  experienced  by  the  gradiometer.  Bell  has  measured 
this  effect  as  a  function  of  frequency  (Ref.  27).  The  depart¬ 
ure  of  gyro-measured  jitter  gradients  from  gradiometer-measured 
jitter  gradients  is  indicated  in  Fig.  A. 2-3.  Note  that  devia¬ 
tion  begins  to  occur  at  10  Hz  and  becomes  serious  at  40  Hz. 


i  to  »  40  eo  no  ioo  toon 

frequency  wt 

Figure  A. 2-3  Bell  Jitter  Compensation  Test  Results 
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Measurement  capability  is  effectively  lost  by  100  Hz.  To  the 
extent  that  these  resonances  are  stable,  some  improvement  will 
be  realized  by  using  shaping  filters  at  the  gyro  output.  How¬ 
ever  only  small  gains  are  to  be  expected  with  this  technique. 

Estimates  of  gradient  errors  (which,  because  of  reso¬ 
nances,  are  "not  seen"  by  the  gyro)  depend  strongly  on  the 
spectral  distribution  of  the  assumed  jitter.  The  following 
two  cases  are  illustrative. 

A . 2 . 5  White  Noise  Platform  Jitter  Model 

A  pessimistic  jitter  model  that  has  a  flat  spectrum 
between  10  Hz  and  100  Hz  can  be  rationalized  as  follows:  Con¬ 
sider  the  fundamental  limit-cycle  frequency  of  20  Hz  to  have  a 
bandwidth  of  ±10  Hz.  Let  the  attenuation  of  the  harmonics  be 
offset  by  the  increased  structural  response  at  higher  frequen¬ 
cies.  Consider  the  case  of  a  structural  response  with  a  cutoff 
frequency  of  100  Hz. 

The  rms  jitter  rate  corresponding  to  Eq .  A. 1-3  is 

ox  =  20/7  7i  —  89  sec/s  (A. 2-23) 

The  flat  PSD  is  obtained  from 

ol  =  2(100  -  10)  S  (A. 2-24) 

A  AA 


which  gives 


Sxx<f>  = 


I  44( sec/s )  /Hz 


1 0  Hz  <  | f  | 
ot  he  rw i se 


100  hz 


(A. 2-23) 


121 


mmmmmmmrn 


If  the  jitter  below  60  Hz  is  compensated,  the  variance  of  the 
uncompensated  portion  is: 


The 


cT  =2 

uncomp 


resulting  DC 


( 100-60 )S  =  3.5  x 

AA 

gradient  measurement 


10"^  ( sec/s 
error  is 


(A. 2-26) 


or  =  82  E  (A. 2-27) 

An  additional  error  results  from  the  spectral  compo¬ 
nents  represented  by  the  convolution  in  Eq.  A. 2-3.  The  result¬ 
ing  error  within  the  gradiometer ' s  passband,  due  to  this  term, 

2 

is  white  noise  with  a  spectral  density  of  173  E  /Hz.  Note 
that  this  value  is  the  error  spectrum  for  the  second  inline 
gradient  term  of  Eq .  A. 2-1.  To  compare  this  figure  with  error 
spectra  for  the  Bell  gradiometer,  it  is  appropriate  to  take 
into  account  the  measurement  provided  by  the  gradiometer ' s 
inline  channel,  Eq .  3.4-12.  When  this  is  done  and  an  equal 
magnitude  of  independent  uncompensated  jitter  is  present  in 
the  y  direction,  the  jitter  error  spectrum  referenced  to  the 
Bell  gradiometer ' s  inline  channel  is  82  E/Hz.  This  is  com¬ 
parable  in  magnitude  to  the  instrument's  self  noise  (see  Fig. 
2.5-1)  and  is  therefore  a  significant  error  source. 

A. 2. 6  Harmonic  Jitter  Model 


Consider  a  platform  limit  cycle  that  is  narrow-band 
at  20  Hz  and,  in  the  worst  case,  is  only  piecewise  continuous. 
The  jitter  spectrum  is  given  by 


S 

xx 


(  f) 


=  A 


2  I  6 ( f - 20n  )  +  6 ( f +20n  )  ] 
ri 


(A. 2-28) 
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where  n=5  corresponds  to  the  100  Hz  cutoff.  The  constant  A  is 
selected  so  that  the  rms  jitter  is  again  given  Eq.  A. 2-23, 


.100  5 


=  2.7  x  10  (sec/s)' 


/i  VV  U  1 

£  -4  [6(f-20n)  +  6(f+20n) ]  df 

-100  n=l  n 


(A. 2-29) 


Neglecting  structural  effects  and  again  assuming  no  gyro  com¬ 
pensation  above  60  Hz  yields  the  uncompensated  jitter  variance 


x( uncomp)  =  2A  £  ?  =  1‘2  X  10  (5Fc/s)‘ 


(A. 2-30) 


This  is  the  appropriate  variance  to  use  in  the  delta- 
function  term  of  Eq.  A. 2-3.  However,  because  the  model  of  Eq . 

A. 2-28  is  comprised  solely  of  impulse  functions,  the  convolution 
of  Eq.  A. 2-3  contributes  additional  zero- frequency  power  to 
the  jitter  gradients.  Since  only  the  zero- frequency  component 
of  the  convolution  falls  within  the  gradiometer ' s  passband , 
only  the  power  at  zero  frequency  is  of  interest.  The  variance 
of  this  uncompensated  jitter  gradient  is 


g( uncorap) 


=  4  A' 


£  K  =  5.2  x  105(£Tc/s)4 
n=3  n 


(A. 2-31) 


Adding  the  v.j^'  term  of  Eq .  A. 2-30  to  the  result  of 
Eq .  A. 2-31  yields  an  rms  value  for  the  total  gradient  bias  of 


a  =  32  E 
g 


(A. 2-32) 
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A. 2. 7  Discussion 

The  values  of  mis  gradiometer  error  given  by  Eq .  A. 2-27 
and  A. 2-31  represent  extremes  for  1  arcseeond  of  platform  jitter 
motion.  The  actual  effect  of  high-frequency  jitter  probably 
falls  somewhere  in  between.  However,  it  is  reasonable  to  con¬ 
clude  that  the  demonstrated  Bell  jitter  compensation  capability 
leaves  room  for  gradiometer  error  of  several  tens  of  eotvos 
(or  more)  for  a  less- than-superb  platform. 

Note  that  some  gradiometer  platform  designs  are  under 
consideration  for  which  the  goal  is  an  inner  element  that  is 
"jitter  free"  when  the  plat  form  is  at  rest.  Although  such  a 
platform  would  provide  a  very  beneficial  gradiometer  environ¬ 
ment,  platform  jitter  errors  may  occur  during  commanded  plat¬ 
form  motion.  During  a  mission  such  as  an  airborne  survey, 
ongoing  small  attitude  changes  ordinarily  associated  with 
flight  could  result  in  essentially  continous  gradiometer 
platform  slewing.  In  such  a  case  jitter  compensation  would 
still  be  requi red . 

A. 3  CONCLUSIONS 

Gyro-based  gradiometer  jitter  compensation  appears  to 
be  viable.  Considerations  such  as  gyro  noise,  bandwidth,  and 
resolution  are  not  serious  issues.  However,  tests  conducted 
to  date  by  Bell  indicate  that  resonances  in  the  gradiometer 
structure  may  be  a  source  of  difficulty.  Since  Bell's  find¬ 
ings  apply  t,o  the  large  gas-bearing  gradiometer  model,  future 
tests  with  the  reduced- s  i  ze  ball-bearing  gradiometer  will  be 
particularly  instructive.  The  key  question  which  remains  to 
be  answered  is  whether  structural  resonances  of  the  gradiometer 
can  be  controlled  at  high  mechanical  frequencies  (>  50  Hz). 
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APPENDIX  B 

MEASUREMENT  EQUATIONS  FOR  AN  IDEALIZED 

- BEIT' GRAVITY  GRADIOMETER - - 


B. 1  INTRODUCTION 

The  Bell  gradiometer  consists  of  four  accelerometers 
mounted  on  a  rotating  platform.  The  axis  of  rotation  is  per¬ 
pendicular  to  the  sensitive  axes  of  the  accelerometers.  The 
orientation  and  reference  frames  to  be  used  are  depicted  in 
Fig.  B.l-1. 


Figure  B.l-1  Coordinate  System  for  Gradiometer  Analysis 

In  this  diagram  the  axis  of  platform  rotation  is  Z. 
The  constant  rate  of  rotation  is  denoted  by  w .  The  four  spe- 
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cific  forces  sensed  by  the  four  accelerometers  are  combined 
into  an  overall  measurement  according  to: 


2  =  fsl  *  fs2  -  £s3  '  fs4 


(B.l-l) 


where 


th 


f  ^  =  specific  force  from  i  accelerometer 


B . 2  MEASUREMENT  EQUATIONS 

B.2.1  Linear  Accelerations 


Generally  , 


(B.2-1) 


where 


I"  h 

g.  =  component  of  gravity  along  1  accelerometer 
sensitive  axis 

a^  =  acceleration  in  the  direction  of  the  i*1*1 
accelerometer  sensitive  axis 


The  basic  process  of  interest  is  mathematically  de 
scribed  by: 


Si  =  So  +  rIj  (B.2-2) 

That  is,  the  gravity  vector  at  the  i1"*1  location  (g.  )  is  de- 
scribed  to  first  order  by  the  sum  of  the  gravity  vector  at  the 
origin  (g  )  plus  the  gravity  gradient  (r)  times  the  position 
vector  ( r .  )  . 

-i 
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The  component  of  gravity  sensed  can  be  expressed  as 
gi  =  ci  gj  (B. 2-3) 

where  is  the  appropriate  matrix  of  elements  to  resolve  . 
Combining  all  of  the  above  yields: 


z  =  cl  <£0+r^l)  '  al  +  c2  (S0+r^2)  "  a2  ‘  c3  (S0+r^3) 

+  a3  '  c4  (S0+r^/4)  +  aA  (B.2-A) 

If  the  platform  is  stabilized  in  inertial  space  (even  if  it 
is  translating) , 


al 

a3 


making 

z  =  c^Fr^  +  c2^-2  ~  c3r I3 
By  direct  calculation  it  can  be 

II  =  -I2  = 


C4r-A 


shown  that 


R  cos  tut 
R  sin  tut 
0 


-3 


-R  sin  mt 
R  cos  utt 
0 


(B.2-5) 


(B. 2-6) 


(B. 2-7 ) ) 


(B.2-8) 
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The  above 


Utilizing 

clr£l 


From  Eqs. 


— 

|  - s i n  uit 

COS  Uit 

0] 

(B.2-9) 

- 

[  -cos  uit 

-sin  uit 

0] 

(B.2-10) 

four  equations  can  be  used  to  reduce  Eq.  B.2-6  to 

z  =  2c^Tr^  -  2c 2^ (B.2-11) 
Eqs.  B.2-7  and  B.2-9  yields 

2  2 

=  R(  -r  cos  uit  sin  uit  +  r  (cos  uit  -  sin  tut) 

1  xx  xy 

+  F  cos  mt  sin  ui 1 1  (B.2-12) 

yy 

B.2-8  and  B.2-10  it  can  be  shown  that 


(B.2-13) 


Thus  , 

2  2 

z  =  4R  1  ( r  -r  )  cos  tut  sin  uit  +  T  (cos  u)t  -  sin  uit)) 
yy  xx  xy 

(B.2-14) 

Now,  the  following  two  identities: 


cos  uit  sin  tut  =  2  sin  2u>t 

2  2 

cos  uit  -  sin  wt  =  cos  2u)t 


(B. 2-15) 


yi  eld 


z  =  2R  (I  -F  )  sin  2uit  4  ART  cos  2u>t 
yy  xx  xy 


(B.2-16)) 


This  equation  describes  the  measurement  obtained  from  the  Bell 
gradioineter  when  its  spin  axis  is  fixed  with  respect  to  inertial 
space . 
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B.2.2  Angular  Velocities 

If  angular  motion  of  the  spin  axis  relative  to  iner¬ 
tial  space  is  present,  the  inertia  forces  generated  by  this 
motion  will  be  sensed  by  the  accelerometers.  For  this  analy¬ 
sis,  let 


f  .  =  c . (g . -a . ) 
si  1  fil  -1 


(B.2-17) 


where 

a.  =  a  +  6  x  r.  +  ft  x  q  x  r.  (B.2-18) 

— 1  — o  —  — 1  -  -  -1 

In  this  equation  Q  and  U  are  the  angular  rates  and  accelera¬ 
tions  of  the  proof  mass  of  the  i1"^  accelerometer.  Using  a 
matrix  notation  for  cross  products, 


fsi  =  ci(£o  +  rIi  -  ao  -  161  II  •  |Q]  V 


(B.2-19) 


where 


in)  = 


[ft]  = 


o 

w 

-Q 

J 

0 

w 

-ft. 


-U) 

0 

ft 

j 

-U) 

0 

ft. 


ft. 


-ft. 


-ft. 


(B.2-20) 


(B.2-21) 


[ft]' 


°x°y 

u)0.. 


V* 

(w2+ft2) 

u»0.. 


u)Q 


u)U 

-<nxH,y> 


(B.2-22) 
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and  ,  Qy  are  the  components  of  angular  velocity  of  the  spin 
axis . 


The  overall  system  measurement  given  by  Eq.  B.2-11 
now  becomes 

z  =  2cj(r  -  [Q 1  -  [Q  ] 2 )  rx  -  2c3(T  -  [ft]  -  [ft]2)  r3 


(B.2-23) 


which  reduces  to 


z  =  2R  (ryy-rxx+ft2-ft2)  sin  2m t  +  4R  (fxy-ftyftx)  cos  2m t 


(B.2-24) 


For  the  special  case  ft^  =  fty  =  0,  this  reduces  to  previous 
result  in  Eq .  B.2-16. 
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APPENDIX  C 


DEFINITIONS  OF  KEY  MATRICES  IN  THE  ERROR  MODEL 
FOR  AN  AIRBORNE  GRADIi'METRIC  SURVEY  SYSTEM 


C.l  INTRODUCTION 

This  appendix  contains  two  major  sections.  In  the 
first  is  a  derivation  of  the  transformation  that  maps  gradio¬ 
meter-system  output  signals  into  the  six  distinct  elements  of 
the  gravity  gradient  tensor.  In  the  last  section  numerical 
values  are  given  for  the  vibration  sensitivity  matrices  in  the 
gradiometer-system  error  model. 


C . 2  RECOVERY  OF  THE  GRAVITY  GRADIENT  FROM 
GRADIOMETER-SYSTEM  OUTPUTS 


In  keeping  with  the  notation  of  Section  3.7,  the  six  out 
put  functions  of  the  gradiometer  system  are  arranged  in  a  vector 


r{(t) 
rc(t> 
r  j(t) 
r2(t) 

r  t ) 


(C.2-1) 
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where  the  superscripts  indicate  the  first,  second,  or  third 
gradiometer  in  the  instrument  triad,  and  the  subscripts  denote 
inline  or  cross-channel  measurements. 


If  the  measurement  system  were  free  of  errors,  then 
the  output  Q  (t)  would  equal  an  ideal  measurement  vector  de¬ 
noted  ['jfl(t),  whose  elements  are  specific  linear  combinations 
of  the  gradient  components  sensed  by  the  individual  instru¬ 
ments  in  the  triad.  The  measurement  vector  Fj^(t)  depends 
both  on  the  gravity  field  and  the  orientations  of  the  three 
grad i omet ers .  To  compute  this  measurement  vector,  let  an  x, 
y,  z  rectangular  coordinate  system  be  chosen  with  the  x  axis 
pointing  east,  the  y  axis  pointing  north,  and  the  z  axis  point¬ 
ing  up.  The  gravity  gradient  at  the  center  of  the  instrument 
triad  can  then  be  represented  by  the  matrix 


XX 

xy 

xz 

yx 

1  yy 

r  yz 

zx 

rzy 

fzz 

(C.2-2) 


where,  e.g.,  =  dg^/dz  dg^/dx,  with  gx  and  gz  being  the  x 

and  z  components  of  the  gravity  vector.  The  six  distinct  com- 

-J- 

ponents  of  this  matrix  can  be  arranged  in  a  column  vector 


r  ( t ) 


XX 


xy 


xz 


I 

yy 

r 

yz 

i 

7.7. 


(t) 
(t  ) 

(r ) 
(t  ) 
(t  ) 
(t  ) 


(C.2-3) 


Note  that  only  five  of  the  gradient  elements  are  independent. 

The  inline  gradients  arc  related  by  I  +1  +1  =0. 

-  xx  yy  zz 
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The  particular  orientations  of  the  three  gradiometers 
in  the  triad  then  uniquely  determine  a  t rans toi mat  ion  matrix  F 
that  maps  the  gradient  vector  F(t)  into  the  ideal  measurement 
vector  r IM( t ) 

=  T~(L>  (C.2-A) 

In  the  system  error  model  of  Section  3.7,  the  first 
gradiometer  in  the  instrument  triad  is  assumed  to  have  a  ver¬ 
tical  spin  axis,  the  second  gradiometer  has  a  south-pointing 
axis,  while  the  third  instrument's  axis  points  east.  The  ideal 
measurement  vector  is  then 


'vvO  -  rxx(t> 

- 1 - 


r .. .( t ) 


rzz<‘>  -  'xx(t) 


rxz(t> 


rzz<[>  -  rvv(t) 


ryz(t> 


and  the  transformation  matrix  T  is 


T  - 


"  h 


0 

0 

0 


0 

1 

0 

0 

0 

0 


0 

0 

0 

1 

0 

0 


0 

0 

0 

I 

"2 

0 


0 

0 

0 

0 

0 

1 


0 

0 

Is 

0 


(C.2-5) 


(C. 2-6) 
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In  practice,  the  outputs  of  the  instrument  triad  (rep¬ 
resented  in  the  ideal  case  by  r^(t))  are  the  observed  quan¬ 
tities,  and  the  gradient  P(t)  is  to  be  computed  from  them. 
However,  the  t ransformat ion  matrix  T  can  not  be  inverted  be¬ 
cause  it  is  singular;  the  sum  of  its  first,  fourth,  and  sixth 
columns  is  the  zero  vector.  Therefore,  to  each  ideal  measure¬ 
ment  vector  in  the  range  space  of  T,  there  exist  an  infinite 

number  of  gradient  vectors  P  that  satisfy  P-^  =  TP.  However, 
only  one  of  these  r  is  admissable,  namely  that  for  which  La¬ 
place’s  equation  for  the  gravity  potential  is  satisfied. 


r  +  r  +  I  =  0 

xx  yy  zz 


(C.2-7) 


Solving  this  equation  for  and  substituting  the 


result  for  f  in  the  expression  for  P[M(t)  yields 


IM' 


r-IM(t>  = 


yy( 1 }  '  rxx(l_> 


r  ( t  ) 

xy 


•|yy(t)  '  2fxx(t) 
- 2 


I  ( t  ) 
xz 


-r  .(t )  -  2i  (t ) 

xx  yy 

" . .  T ^  - 

r  ( t ) 
yz v 


(C. 2-8) 


Since  I  ( l )  no  longer  appears  in  1  JM(t),  the  ideal 
measurement  vector  can  be  computed  from  a  reduced  gradient 
vector  P^(l),  which  is  obtained  by  deleting  the  r  (t)  com- 
porien  t  t  rom  I  (  t  )  : 


]  14 


(C  .2-9) 


The  corresponding  reduced  transformation  matrix  Tp 
satisfies  fj^ft)  =  T^r^(t)  and  is  given  as 

0  0  ^0 

10  0  0 
0  0  -!*  0 

(C. 2-10) 

0  1  0  (, 

0  0-10 

0  0  0  1 

The  reduced  transformation  matrix  Tp  is  also  singular 
(because  it  is  not  square),  but  it  has  full  column  rank.  There¬ 
fore,  each  ideal  measurement  that  satisfies  Lj^(t)  =  T^L^(b) 

corresponds  to  exactly  one  reduced  gradient  vector  fp(t). 

This  unique  vector  can  be  computed  as 

rR( t )  =  Tp  rJM(t )  (C.2-11) 

where  Tp  is  the  pseudo  inverse  of  Tp 

Tp  =  (Tj  Tp)*1  Tj  (C.2-12) 

As  a  result,  the  full  gradient  vector 


I  35 


r  ( t ) 


'H(t) 


r 

zz 


<t) 


(C. 2-13) 


satisfies 


I  ( t ) 


To  ~-IM 


(t) 


where  the  transformation  matrix  Tq  is  given  by 


(T^  T  ) ' ^ 

T^ 

1 - 

o 

0  2/3  0 

2/3  0 

(C. 2-14) 


(C.2-15) 


C.3  G RADIOMETER  SYSTEM  VIBRATION  SENSITIVITIES 

This  section  gives  estimated  numerical  values  for 
both  the  linear  and  angular  vibration  sensitivity  matrices, 
Siin^r)  and  ^anc  (  *k )  ’  as  in  Section  3.7.7.  These 

estimates  are  based  on  test  data  as  summarized  in  Tables  3.3-1 
and  3.4-1 . 


Seiner)  c°nsidered  first.  This  6x3  matrix  of 
mean-square  sensitivities  to  linear  vibration  depends  on  the 
vibration  frequency  f^  =  kf  ,  k=0,l,2,3,  where  the  spin  fre¬ 
quency  f  ~  0.28  Hz.  Equation  3.7-3  is  used  to  associate  a 


specific  entry  in  Table  3.3-1  with  each  element  of  ^^(f^)- 
The  units  of  SIIM(f,,)  are  EVg^ 


El  N 


therefore,  the  entries  in 

i  it  i\ 

able  3.3-1  are  mult  iplied  by  1000  mg/g  and  then  squared, 
h  i  :>  y  i  e  i  ds 


1  3b 


SUN  <IIVSA.<0> 


6.8  x  10 


P2,  2 
E  /g 


SLIN  (IIVSA,  f| )  =  1.7  x  10* 


SLIN  (HVSA.f2) 


=  6.2  x  10' 


'l  .  2 
E  /g 


'2.  2 
E  /g 


(C. 3-1 > 


SL1N  (HVSA<f3)  =  2.8  X  10 


r2  ,  2 

E  /g 


SLIM  UVSA’fo> 


2.1  x  105  E2/g2 


SL1N  <1VSA’fl>  = 
SLIN  (IVSA,^)  = 


1.8  x  107  E2/g2 


2.4  x  103  E2/ g2 


(C.3-2) 


SL1N  (J-VSA>  f 3) 


1.9  x  107  E2/g2 


S21n  (1IHSA,  fk)  =  S2in  { II VSA  ,  f  k  )  for  k  =  0,1. 2, 3  (C.3-3) 


SLIN  (-LHSA,f0) 


2.1  x  105  E2/ g2 


SLIN  (lHSA’fi> 


SkIN  (1.HSA , f 2 ) 


6.6  x  103  E2/g2 


s  ?  9 

3.2  x  10J  EVgZ 


(C.3-4) 


SL1N  (2-HSA-f3> 


1.4  x  10^  E2/g2 


The  same  procedure  is  used  to  determine  numerical 
values  for  the  elements  of  f ^ )  in  Eq .  3.7-15,  except  that 

Table  3.4-1  is  used,  and  its  scale  factor  of  10^  must  be  taken 
into  account.  The  results  for  k=l,2,3,4  are 
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(l|VSA,fA) 

=  5 . 8  x 
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UVSA,f1) 

=  9 . 4  x 
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E2/( rad/s )2 
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E2/( rad/s )2 
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1  0'^ 

E2/( rad/s )2 
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)  for  k=l , 2 , 3 , 9 
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APPENDIX  D 

INTENTIONAL  GRAD IOMETER  PLATFORM  MOTION 


D.l  INTRODUCTION 

Most  studies  of  gradiometer/inertial  platform  mech¬ 
anizations  visualize  a  triad  of  gradiometers  attached  to  the 
inner  element  of  a  local-level  gimballed  platform.  The  ge¬ 
ometry  of  the  triad  involves  a  symmetric  distribution  of  the 
three  instruments  about  the  vertical  with  each  gradiometer 
inclined  at  the  same  "umbrella"  angle.  This  "conventionally- 
mechanized"  gradiometer  is  illustrated  in  Fig.  D.l-1.  Figure 
D.l-1  also  defines  the  gradiometer  body  axes  u,  v,  and  w  and 
the  local-level  coordinates  x,  y,  z. 

In  the  conventional  mechanization,  the  platform's 
inner  element  is  held  fixed  about  the  two  horizontal  axes  and 

JL. 

in  azimuth.  However,  as  illustrated  in  the  following,  speci¬ 
fic  advantages  can  be  gained  by  superimposing  platform  rota¬ 
tion  motions  onto  the  nominal  local-level  platform  orientation. 
Such  motions  have  the  potential  for  decreasing  the  effective 
noise  of  the  gradiometer  and  reducing  the  number  of  gradiome¬ 
ters  needed  for  completely  measuring  the  gradient  tensor. 
Several  variations  exist.  The  simplest  is  constant  rate  ro¬ 
tation  or  "carousel 1 ing"  about  the  vertical  axis.  Carousel- 
ling  is  discussed  in  detail  in  Sections  D.2  and  D.3.  The  other 
rotation  concepts  have  similar,  although  more  mathematically 


'>Azimuth  wander  mechanizations  would  also  be  suitable  for 
gradiometers . 


L39 


P-42&0U 


PLATFORM 
WWf»  LUMIN  7 


Figure  D.l-1  Gradiometer  Triad  Element 

involved,  descriptions.  Compound  mechan izat i ons  involving  tum¬ 
ble  or  rocking  motion  of  the  plat  form  (in  addition  to  carousel- 
ling)  are  treated  in  Section  D.4.  Section  D.5  concludes  with 
a  description  of  step- index i ng  ,  a  discrete  form  of  carousel  1 ing . 

I) .  2  OUTPUT  EQUATIONS  FOR  A  SINGLE  CAROUSEL- 
M ECU AN  I  ZED  CHAD  IOME  TER 

I) .  2  .  1  Coordinate  Transformations 

Three  orthogonal  coordinate  frames  are  of  interest: 
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(1)  The  u,  v,  w  gradiometer  body  axes  of  Fig. 
D. 1-1  . 


(2)  A  frame  (x^,  ,  z)  associated  with  the 

(vertical)  carousel  axic  (distinct  from  the 
gradiometer ' s  internal  spin  or  weights  axis). 
The  (x^,  y^ ,  z)  frame  rotates  synchronously 
at  the  carousel  rate. 


(3)  The  frames  (u,  v,  w)  and  (x^ ,  y^ ,  z)  are 

both  rotating.  Let  (x,  y,  z)  be  the  fixed 
local-level  frame  which  corresponds  to  (Xj  , 
y^ ,  z)  at  a  specified  time. 

Relations  among  the  coordinates  are  summarized  as  follows: 


Transformation  Between  Gradiometer  Body  F rame 
and  Carousel  Frame  (See  Figure  I).  2-1) 


Figure  D.2-1  Gradiometer  Body  Axes  to  Carousel  Axes 


0  . 


The  transformat  ion  is  given  by 
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0 
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So 
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(D.2-1  ) 


where  the  abbreviations  C  and  S  are  used  to  indicate  sine  and 
cosine . 


Transformation  Between  Carousel  Frame  and  Fixed  Frame 
Yl 

x 

Y  /  z  is  common  (out  of  paper) 

X  /  B  =  ui  t  ,  u>  =  carousel  rate 

\  /  _  P  P 


C2  yl 
2 


Figure  D.2-2  Gradiometer  Carousel  Axes  to  Fixed  Frame  Axes 


The  transformation  is 
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(D. 2-3) 


The  gradients  in  the  fixed  frame  (x,  y,  z)  are  related 
to  the  gradients  in  the  grad i ome ter ' s  body  frame  (u,  v,  w)  by 

I  =  r2  |  r° 

i  ()  i.Q  i  2 


(D.2-4) 


where 


C2  =  (C®)1 


and  the  gradient  matrix  is  represented  by 
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(D.2-6) 


The  superscript  T  indicates  matrix  transpose. 

D .  2 . 2  Bell  Gradiometer 

A  single  rotating  accelerometer  gradiometer  provides 
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the  two  measurements 


m, D  =  P  -  V 
IB  uu  vv 


m..d  =  r 

2B  uv 


(D.2-7) 
(D. 2-8) 


When  the  expression  indicated  in  (D.2-4)  is  evaluated  and  sub¬ 
stituted  into  (D.2-7)  and  (D.2-8),  the  following  output  equa¬ 
tions  of  a  single  gradiometer  are  obtained  after  taking  advan¬ 
tage  of  several  trigonometric  identities. 

m,K  =  (C2B  -  C2a  S2B)T  +  (S2B  -  C2orC2B)r  -  S2af 
IB  xx  yy  zz 


-  (1+C a)  S2Br  +  S2«SBr  +  S2aCBT 

xy  xz  yz 


(D.2-9) 


*In  practice,  the  gradiometer  outputs  are  modulated  by  the  in¬ 
strument  block's  spin  frequency  (~l/6  Hz).  Tn  this  develop¬ 
ment  the  grad i en t - s i gna 1  and  carousel -carrier  frequencies 
being  considered  are  much  lower  than  the  gradiometer  spin 
rate.  Hence,  it  is  appropriate  to  avoid  the  additional 
complexity  of  the  extra  spin  frequency  terms. 
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«2B  "  2  CuS2BI  xx  '  2  Casm'yy  *  CaC2BI'xy  -  SaCBP^  +  SaSBryz 

(D.2-10) 

It  is  useful  to  express  the  two  measurements  in  terms  of  their 
components  at  zero  frequency  and  quadrature  components  at  and  2uip. 

n'lB  =  n,0B  *  mCA  +  mSl  *  mC2  +  mS2  (D.2-11) 

_  2B  ^  2B  ^  2B  ^  2B  ^  2B  /rk  „  .  „  . 

,U2B  '  m0  +  n,Cl  +  ,nSl  +  mC2  +  mS2  (D.2-12) 


Each  frequency  component  is  given  by 


1  B 
m() 

3  ...  .2 

-  -  ^  1 zz  s 1 n  « 

(D.2-13) 

IB 

mci 

=  T  s  i  ri  2a  cos  w  t 

yz  p 

(D.2-14) 

IB 

msi 

-  1  sin  2a  sin  uj  t 

XZ  p 

(D.2-15) 

IB 

IUC2 

=  ^(r  -  r  )(i  +  cos2  «)cos  2w  t 

2  xx  yy  p 

(D. 2-16) 

1  B 
mS2 

-  f ' x y f  1  +  cos2  a)  sin  2w^t 

(D.2-17) 

2B 

m0 

-  0 

(D. 2-18) 

2B 

mci 

=  -T  sin  a  cos  uj  t 

xz  p 

(D. 2-19) 

2B 

mSl 

=  1  sin  »  sin  m  t 

yz  p 

(D.2-20) 

2  B 
"'C2 

=  1  cos  a  cos  2ui  t 

xy  p 

(D.2-21 ) 

2B 

"*S2 

-  ->(l  v  -  1  )  cos  a  sin  2uj  t 

2  xx  yy  p 

(D.2-22) 

''here  i.aplace's  ec|u.it  ion 
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r  +i  +i  -  o 

xx  yy  zz 


(0.2-23) 


has  been  used  to  simplify  (D.2-13). 


D . 2 . 3  Draper  Gradiometer 


A  single  floated  gradiometer  provides  the  two  measure¬ 


ments 


mlr.  =  r 

ID  vw 

2D  uw 


(D. 2-2A) 
(D.2-25) 


Following  the  same  procedure  as  in  the  previous  section 

m1n  =  SaCaS2Br  +  SaCaSBCBI  -  S2oSBr  +  SaCoSBCBP 
ID  xx  xy  ?.x  xy 

+  SaCaC2Br  -  S2aCBI  +  C2aSBI  +  C2aCBP 

yy  yz  xz  yz 


-  SaCaP 


zz 


(D.2-26) 


m0r.  =  SaSBCBr  -  SaS2BP  +  SaC2BP  -  So  SBC  BP 
2D  xx  yx  xy  yy 


+  CaCBT  -  CaSBT 

xz  yz 


(D.2-27) 


The  individual  frequency  components,  defined  in  the  same  fashion 
as  (D.2-11)  and  (D.2-12),  are 
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T  cos  2a  sin  w  t 
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(D.2-28) 

(D.2-29) 

(0.2-30) 

(  D .  2  -  3  1  ) 


1A5 


* 

11) 

111 S  2  r 

.1  1  sin  2«  sin  2w  t 

2  xy  p 

(D. 2-32) 

m21)  = 

0 

(D. 2-33) 
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r  cos  <*  cos  Ul  t 

xz  p 

(D. 2-34) 

21) 

msi 

-T  ins  a  sin  in  t 

yz  p 

(D. 2-35) 

21) 

n,C2 

1  sin  ii  cos  2ui  t 

xy  p 

(D. 2-36) 

2D 

mS2  ' 

J(l  -  1  )  sin  «  sin  2w  t 

2  xx  zz  p 

(D.2-37) 

Again  t  hr  Laplace  relation  (I). 2-23)  has  been  used  in  obtaining 
(D.2-2H)  . 

1) .  3  DISCUSSION 

D . 3 .  1  Si  gnu  1  Keeove  rv 

From  equal  ions  (I). 2-13)  through  (I). 2-17)  it  is  evi¬ 
dent  that  the  grad  i  <>me  t  e  t  '  s  inline  channel  (n\j)  provides  in¬ 
dependent  ampl i t ude-modu 1  at ed  outputs  containing  all  of  the 
information  required  to  specify  the  gravity  gradient  tensor. 

The  ^radiometer's  cross  channel  provides  only  four  independent 
elements  of  tensor  information;  they  are  (with  the  exception 
of  I  t  tin  *,«,me  as  provided  by  the  inline  channel.  However, 
the  cross  channel  outputs  have  a  different  scale  factor  depend¬ 
ency  on  the  amine  I  la  angle,  <v  .  As  a  result,  for  certain  values 
of  n,  tie  cn  c.  channel  measurements  cense  to  be  completely 
r  i  >  Iiiiii  I  an  t  *  1  h  the  lullin’  measurements.  Such  a  situation  is 

•li  •  u  •  •  so.  I  i  ,  Set  t  i  on  I) .  A  similar  case  holds  for  the 

two  out  | .  *  i  *  I  i  he  Inapt  i  g  i  ad  i  olio- 1  e  r  .  Note  that  the  output 

a  I  .1  I  on  se  I  <  <!  Hell  g  I  ad  i  <  line  ter  |  (  I)  .  2  -  1  3  ) 


1  a  f  i 


A 


•  t  p  I  a  I  I  t  a  i  t  t  mic  -  I  hi 


through  (D. 2-22)  |  provide  the  k<iiik-  information  as  those  for  a 
carouseled  Draper  gradiometer  ((D.2-28)  through  (D.2-32)). 
Only  the  scale  factors  differ. 


It  is  interesting  to  note  the  effects  of  orienting  a 
carousel  axis  in  a  direction  other  than  vertical.  Although 
the  equations  associated  with  the  additional  angular  degree  of 
freedom  are  not  reproduced  here,  the  following  findings  are 
str  ad. 


•  The  zero- f requency  (i.e.,  unmodulated) 
element  is  the  gravity  gradient  defined 
by  the  carousel  axis 

•  Sufficient  information  is  available  from 

a  single  gradiometer  instrument  to  specify 
the  complete  gradient  tensor  in  any  coor¬ 
dinate  frame  (subject  to  the  demodulation 
constraints  discussed  in  Section  D.d.M 

•  Simultaneous  carousel  ling  of  three  gradi- 
ometers  on  a  common  platform  poses  no 
theoretical  difficulty.  Each  gradiometer 
provides  a  set  of  independent  outputs  of 
the  form  D.2-9  and  I). 2-10  or  D.2-26  and 
D.2-27. 


D . 3 . 2  Number  of  G radiometers  Kequ i red  for  Tensor  Recovery 
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If  the  carousel  frequency,  is  sufficiently  high, 

then  the  gradient  signals  of  Eqs.  (D.2-9)  and  (D.2-10)  or  (D.2-26) 
and  (D.2-27)  can  be  recovered  by  conventional  demodulation  tech¬ 
niques.  The  five  independently  recoverable  signals  provided 
by  a  carouselled  gradiometer,  together  with  Laplace's  equation, 
completely  specify  the  gravity  gradient  tensor.  Thus,  in  prin¬ 
ciple,  only  a  single  instrument  (instead  of  a  triad)  is  required 
to  provide  a  full  set  of  gradient  measurements. 


“Criteria  for  selection  of  m  are  discussed  at  further  length 
in  Section  D . 3 . 6 .  * ’ 
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It  is  interesting  to  observe  the  simplified  form  which 
the  demodulated  output  equations  of  a  Bell  gradiometer  take 
when  the  umbrella  angle  is  horizontal  (a  =  90°).  They  are 
summarized  below. 

0 

r 

xz 

r 

yz 

0 

0 

The  scale  (actors  of  the  output  terms  aiso  simplify  and  only 
two  redundant  terms  remain  in  the  cross  channel  signal.  Such 
a  configuration  could  ultimately  be  attractive  in  a  carousel- 
mechanized  system. 

I) .  5 .  1  Not  se  Reduc  t  i  on 

The  frequency  content  of  gradiometer  self-noise  is 

treated  extensively  in  the  main  section  of  this  report.  For 

missions  requiring  gradient  recovery  at  frequencies  below  about 

0.1  mHz  t he  instrument  noise  at  low  frequencies  (red  noise) 

can  impose  significant  performance  limitations.  The  potential 

advantage  of  carousel ing  in  reducing  low-frequency  noise  can 

be  seen  if  the  carousel  r.itc,  w  .  is  selected  to  lie  above  the 

P 

red-noise  frequency  range  of  the  grad i ome t e r ' s  error  spectrum, 

■A  s  i  m  i  l /i  r  degenerate  form  occurs  tor  the  Draper  gradiometer 
equations.  However,  the  carousel  orientation  having  the  axis 
<>f  weights  horizontal  is  of  less  interest  for  the  Draper  in¬ 
strument  because  the  I  term  becomes  unobservable. 
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(e.g.,  Wp  >  0.6  mrad/'s )  .  When  t  he  signal  is  lowpass  demodu- 
lated,  the  only  gradiometer  errors  which  remain  in  the  output 

.t, 

are  those  originally  within  a  narrow  band  about  Thus,  from 

the  carouselled  gradiometer  output  equations,^  it  is  apparent 

that  the  red  noise  can  be  suppressed  in  all  of  t he  gradient 

elements  except  F 

r  zz 

No  direct  improvement  in  F  occurs  because  it  is 
r  zz 

unmodulated  by  the  carousel  action.  Note  that,  while  low- 

frequency  gradiometer  errors  are  not  eliminated  from  1 '  ,  some 

improvement  occurs  by  virtue  of  the  improvement  in  the  (r  - 

T  )  term.  This  results  because  (F  -  I  ),  in  combination 
yy  xx  yy 

with  Laplace's  equation,  is  used  to  recover  F  .  However,  the 

inability  of  carouselling  to  eliminate  the  majority  of  the  red 

noise  from  T  (and,  as  a  result,  from  the  other  two  inline 
zz 

gradients)  is  a  deficiency.  Extensions  which  allow  the  I 
term  to  be  modulated  also  are  discussed  in  Section  0.4. 

Note  that  these  comments  still  apply  if  more  than  one 
gradiometer  is  carouselled  on  the  same  platform.  Regardless 
of  the  orientation  of  any  additional  grad iometers ,  there  is 
still  no  motion  of  any  instruments  in  the  z  direction  (carousel 
axis).  Thus  there  is  no  basis  for  modulation  of  F  .  01  course, 
the  addition  of  instruments  is  always  advantageous  due  to  the 
reduction  of  noise  inherent  with  the  averaging  of  separate 
instrument  outputs. 


'"This  conceptual  discussion  omits  important  gradiometer  signal 
bandwidth  considerations.  These  are  taken  up  in  the  next 
sec  t ion . 


t  Equations  0.2-13  through  D.2-22  ( Be- 1  1  )  and  0.2-28  through 
0.2-32  ( Orape r ) . 


■» 
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I) .  3 .  A  ('a  rouse  1  Speed  arid  Out  put  Bandw  i  d  th  Considerations 

Two  benefits  from  carousel  ling  a  grad iome ter  have 
been  identified,  namely  the  capability  to  recover  the  gravity 
gradient  tensor  with  a  single  gradiometer  instrument  (instead 
of  a  triad)  and  low-frequency  noise  reduction.  This  section 
examines  some  ol  the  different  carousel -speed  and  operating- 
bandwidth  constraints  associated  with  each. 


To  utilize  a  single  gradiometer  for  recovery  of  all 
ol  the  tensor  elements,  the  carousel  frequency,  ,  must  be 
somewhat  larger  than  the  highest  frequency  contained  in  the 
gradient  signal,  w  .  Otherwise  information  loss  results. 

This  is  illustrated  in  fig.  I).  3-1  for  the  case  where  the  band¬ 


width  of  the  gradient  is  slightly  less  than  Note  that 

Fig.  If. 3-1  depicts  the  process  of  amplitude  modulation  and 
demodulation  with  ui^  as  the  carrier.  The  low-frequency  noise 
suppression  inherent  in  this  process  is  also  illustrated. 

Note  also  that  il  the  carousel  frequency  is  less  than  wcq ,  the 
distinct  lobes  shown  in  Fig.  11. 3-1  overlap.  In  this  situation, 
if  the  lowpass  filter  were  set  to  pass  all  frequencies  up  to 
uj(  ,  error  would  result  from  the  sideband  overlap.  On  the 
other  hand,  if  the  filter  bandwidth  were  reduced  to  omit  the 
overlap  region  from  the  output  ,  the  high-frequency  portion  of 
the  signal  w  i  be  lost. 


The  va I ue  of  w  depends  on  the  nature  (e.g. ,  rough 

or  smooth)  o f  the  local  gravity  field  and  the  speed  of  the 

g rad i ome t e r ' s  host  vehicle.  A  rule  of  thumb  is  to  consider 

the  cutoff  wavelength  of  the  earth’s  field  to  be  10  km.  Then, 

t  lie  minimum  carousel  frequency  which  meets  the  criterion  that 

m  ui  is  given  by 
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Which  Avoids  Signal  Loss 
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where  uj  is  given 
P 

vehicle  such  as  an 


in  rad/ s  and  v  is  in  km/s .  For  a  slow- speed 
ocean  vessel  moving  at  2  m/s  (3.9  kts),  the 


minimum  carousel  frequency  evaluates  to  0.72  cycle  per  hour. 


For  a  rapidly  moving  vehicle  such  as  an  aircraft  trav 
eling  at  130  m/s  (291  kts),  the  minimum  value  of  the  carousel 
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rate  works  out  to  0.9  cycles  per  minute.  As  seen  later  this 
value  is  significantly  beyond  practical  limits. 


If  multiple  gradiometers  (e.g.,  a  triad)  are  carou- 
seled,  the  value  of  required  to  effect  low  frequency  noise 
reduction  is  much  lower.  This  can  be  illustrated  by  onsider- 
ing  a  single  quadrature  component  of  the  output  of  oi  carou¬ 
sel  led  grad iome t e r  ,  e.g., 


Sj(t)  -  F(t)  cos  u>pt  vD.3-2) 

where  !'(t)  is  the  gradient  element  measured  by  that  component. 
The  output  of  a  neighboring  gradiometer,  aligned  at  an  angle  0 
with  respect  to  the  first  is  given  by 


S„(t)  =  r  ( t  )  cos  (u>  t  +  0)  (D.3-3) 

P 

With  both  S-^(t)  and  $2(0  available,  the  following  process  il¬ 
lustrates  the  recovery  of  T(t)  at  frequencies  well  beyond 

If  the  first  signal  ( Eq .  D.3-2)  is  divided  by  cos  w^t 
at  points  other  than  the  zeros  of  the  cosine  function,  the 
piecewise  function 


'In  practice,  for  an  aircraft  flying  at  high  altitude,  a  lower 
value  of  would  be  permissible,  since  the  upward  continued 

gradient  field  would  cut  off  at  a  longer  wavelength  than  10  km. 
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Sj(t) 

cos  w  t 
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U1  t 
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/  nn ,  n=0 , 1 , 2 . 


(I).  3-0  ) 


results.  A  similar  division  is  carried  out  for  the  second 
signal  (Eq.  D. 3-3)  . 


s2(t) 


S®  / i \  -  _ 

z  '  '  cos(ui  t  +  <|>  ) 


(uipt+0)  /  nn  ,  n=0,l,2... 


(D.3-5) 


Of  course,  in  an  actual  implementation,  the  gaps  would  be  of 
finite  width  due  to  noise,  quantization,  etc.  The  multiple 
gradiometer  demodulation  process  is  shown  in  Fig.  D.3-2  where 
the  gap  width  is,  for  illustration  purposes,  indicated  as  90°. 
Also,  for  simplicity,  the  figure  illustrates  the  case  where 
the  phase  angle,  $ ,  is  also  90  degrees.  In  practice,  the  right- 
angle  gradiometer  arrangement  would  be  adverse  to  demodulation 
of  the  2u>p  signals.  It  is  use  here  only  for  i  1  1  us t ra t  i vi*  con¬ 
venience.  In  fact  a  gradiometer  triad  as  illustrated  in  Fig. 
D.l-1  is  the  appropriate  geometry. 


For  purposes  of  examining  the  frequency  effects  of  the 
demodulating  process  described  in  the  previous  paragraph,  it 
is  useful  to  generalize  Eqs.  I). 3-0  and  1).3-T>  as  finite  pulse- 
width  samplings  of  the  modulated  signal.  Analysis,  not  pro¬ 
vided  here,  shows  that  the  piecewise  demodulated  signals  are 
subject  to  edge  errors  which  degrade  the  recovered  gradient 
signal.  This  effect  becomes  severe  only  at  signal  frequencies 
substantially  beyond  the  carousel  frequency  and  can  be  miti¬ 
gated  somewhat  by  the  use  of  a  tapered  window  to  perform  the 
pu 1 se  samp  ling. 


An  alternative  technique  for  demodulating  S j ( t  )  and 
S^(t)  is  to  formulate  the  demodulated  signal  as  the  output  of 
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a  Kalman  filter  estimator.  In  this  case  the  observation  equa¬ 
tions  (time  varying)  contain  the  modulation  frequency  and  requi¬ 
site  statistical  models  for  the  gradient  signal  and  gradiometer 
noise  processes.  As  with  the  conventional  filter  technique,  the 
Kalman  approach  has  a  high  error  sensitivity  to  signal  content 
at  frequencies  which  are  considerably  beyond  the  carousel  fre¬ 
quency.  This  motivates  as  high  a  carousel  rate  as  possible. 

D . 3 . 5  Error  Sources  Limiting  Carousel  Rates 

JL. 

The  chief  limitation  on  carousel  rates  is  the  cen¬ 
trifugal  gradient  field  which  the  carousel  motion  induces. 

The  dynamic  range  of  current  gradiometer  instruments  is  of  the 

order  of  10^  E.  The  gradiometer  bias  induced  by  w  (measured 

P 

in  rad/s)  is  given  by 

re  =  wp  x  I0<)  E  (0-3.6) 

For  the  case  of  a  bias  offset  of  2*10^  E  (i.e.,  twice  the  in¬ 
strument's  dynamic  range),  the  corresponding  carousel  rate  is 
2.6  cycles/hour.  To  maintain  a  bias  stability  error  of  less 
than  1  E,  the  bias  offset  setting  is  required  to  be  stable  to 
50  parts  per  million  (ppm).  Similarly,  to  maintain  an  inertial 
gradient  field  constant  to  1  K,  the  carousel  speed  control  is 
required  to  hold  a  constant  rate  within  25  ppm.  These  figures 
suggest  that  implementing  continuous  carousel  speeds  higher 
than  about  3  cycles/hour  would  prove  formidable. 

If  the  2.6  cycles/hour  figure  is  considered  to  be  a 
maximum  feasible  carousel  rate,  it  is  of  interest  to  estimate 


"An  approach  for  removing  this  limilat  ion  is  discussed  in 
Sec  t ion  D . 5 . 
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the  highest  vehicle  velocity  consistent  with  single  gradiome- 
ter  (vs  a  triad)  operation  and  recovery  of  the  complete  gra¬ 
dient  tensor.  Solving  Eq .  D.3-1  for  v  yields  the  value  7.1 
m/s .  This  suggests  that  the  use  of  a  single  carouselled  gra- 
diometer  instead  of  a  full  triad  on  board  an  aircraft  is  well 
beyond  the  state  of  the  art. 

D . A  GENERALIZATIONS  OF  CAROUSEL  MOTION 

In  Section  D.3  it  was  observed  that  carouselling  pro¬ 
vides  no  modulation  of  the  gradient  tensor  element  defined  by 
the  carousel  axis.  As  a  result  low-frequency  gradiometer  noise 
continues  to  be  a  major  error  source  for  inline  gradients  (for 
the  gradient  tensor  expressed  in  the  reference  frame  of  the 
carousel  motion).  For  obliquely  defined  frames,  all  of  the 
gradient  elements  can  contain  significant  amounts  of  the  gra- 
dioineter's  red  noise.  The  techniques  discussed  in  this  sec¬ 
tion  are  directed  toward  implementing  intentional  platform 
motion  which  provides  a  suitable  modulating  signal  on  the 
carouse  1  -  ax i s  gradient. 

D.4.1  Tumble  Motion 

If,  in  addition  to  carouselling,  constant  -  rate  rot a- 
t ion  of  the  gradiometer  platform  is  implemented  about  an  axis 
orthogonal  to  the  carousel  axis,  the  desired  additional  modu¬ 
lation  is  effected.  The  generalization  of  the  carousel  rate 
output  equations  of  Section  I) .  2  to  include  t  he  effects  of  tum¬ 
ble  mot  ion  proceeds  similarly  to  the  derivation  of  that  section. 
For  the  umhrel la  angle,  a ,  given  by 

a  -■  u»t  t  (D.4-1) 
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where  u)^  is  the  tumble  rate,  the  outputs  equations  for  the 
inline  channel  of  a  carousel  and  tumble  mechanized  Bell  gra- 
diometer  are 
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where  the  quadrature  terms  for  each  frequency  sum  to  provide 
the  total  inline  output  in  the  same  fashion  as  F.q .  D.2-11.  As 
before  the  m  superscript  1  identifies  the  inline  channel  output 
and  B  designates  that  the  measurement  is  associated  with  t he 


157 


Bell  gradiometer.  The  subscripts  indicate  the  specific  quad¬ 
rature  component  involved. 

For  the  cross  channel  a  different  set  of  frequency 
components  is  obtained. 
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Note  that  while  the  inline-channel  outputs  alone  are 
sufficient  to  completely  specily  the  entire  gradient  tensor  if 
uij  and  u>  are  not  harmonically  related,  the  cross-channel  terms 
are  not  .  This  is  not  surprising  since  it  also  occurs  for  carou¬ 
sel  motion  alone.  However  all  of  t he  distinct  frequency  signals 
yield  a  richness  of  redundant  informat  ion  which  provides,  in 
add i t  i on  to  the  desired  modulation  of  F  ,  the  extra  capability 
for  noise  reduction  by  averaging. 
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Unfortunately  the  tumble  mechanization  is  not  prac¬ 
tical  for  current  gradiometer  instruments  because  of  the  time 
varying  stresses  induced  by  the  earth's  gravity  field.  How¬ 
ever,  key  elements  of  it  are  retained  in  the  platform  rock 
mechanization  discussed  in  the  next  section. 

D . 4 . 2  Platform  Rock  Motion 

To  reduce  the  adverse  gravi ty- induced  effects  associ¬ 
ated  with  tumbling  the  gradiometer,  consider  low-amplitude 
platform  rocking  about  the  nominal  umbrella  angle,  a  .  The 
motion  is  given  by 


a(t)  -  a  +  a  sin  ut  t  (D.4-22) 

or  r 

where  a^.  is  the  rocking  amplitude  and  is  the  rocking  rate. 

The  rocking  motion  introduces  frequency  modulation  of  the  gradi¬ 
ometer  outputs  as  a  result  of  terms  of  the  form  sin  (a^  sin  wrt  )  . 
Expansion  of  the  inline-channel  equations  for  the  Bell  gradi¬ 
ometer  and  the  use  of  suitable  trigonometric  identities  pro¬ 
vides  the  following  frequency  and  phase-separable  terms.  The 
notation  is  defined  in  the  same  fashion  as  earlier. 
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where  J  (2a  )  is  the  Bessel  function  of  the  first  kind,  order 

n  r 

n.  A  similar  set  of  equations  occurs  when  the  Draper  gradi- 
ometer  is  rock  mechanized. 


Note  that  there  is  sufficient  information,  separated 
in  frequency,  to  recover  all  elements  of  the  gradient  tensor 
with  the  low  frequency  noise  suppressed.  Note,  too,  that  be¬ 
cause  the  rock  amplitude  a ^  is  small  the  coefficient  of  each 
harmonic  term  (i.e.,  the  scale  factor)  is  also  small.  Th i s 
mot i vat  es 


As  large  a  rock  amplitude  as  possible, 
consistent  with  minimizing  g-induced 
gradiometer  errors 

Retention  of  as  many  sideband  terms  as 
possible  in  the  demodulation  process. 


■As  with  the  tumbling  mechanization,  w  and  w  must  not  be 
harmonically  related.  !  1 
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I).1)  STF.P  INDKXINc;  MECHANIZATION 

It  has  already  been  observed  that  one  disadvantage  of 
the  carousel  mechanization  is  the  effective  rate  limitation  im 
posed  by  the  associated  centrifugal  gradient  field.  A  suppres 
sion  technique  is  to  alternately  rotate  and  stop  the  gradiom- 
eter.  During  the  rotation,  while  the  instrument  is  slewed 
through  a  portion  of  a  revolution,  it's  output  is  disregarded. 
Once  the  slew  motion  is  completed,  the  gradiometer  remains  at 
rest  until  the  transients  resulting  from  the  slew  have  settled 
The  instrument  output  is  then  sampled,  after  which  another 
slew  cycle  begins.  This  mode  of  operation  is  illustrated  in 
Fig.  D.r)-1  for  the  case  of  A r)  degree  slew  steps  and  equal 
length  dwell  periods.  In  the  figure  the  last  portion  cj  f  the 
dwell  cycle  is  used  for  sampling. 

I"  A  I  M)HM 
UAH 


I  (Mf 

F  i  gu  rt-  I ) .  >  -  1  C.  r  ad  i  <  >me  l  <•  r  I’l.it  lorm  Motion  in  Step 

Indexing  Mri  li.iinz.il  inn 


1  b2 


rm 


Because  the  sampling  occurs  only  when  the  ^radiometer 
is  motionless,  no  biases  are  introduced  to  t he  ^radiometer  in 
the  course  of  the  step  indexing  motion.  Similarly,  there  is 
no  accuracy  requirement  for  .  (However,  the  sampling  period 
must  be  maintained  and  known  accurately.) 

The  upper  limitation  on  the  effective  carousel  rate 
of  this  mechanization  is  dictated  by  the  settling  time  re¬ 
quired  for  a  given  slew  speed.  Preliminary  tests  have  been 
made  by  Bell  which  indicate  good  gradiometer  performance  with 
45-deg  slew  steps,  a  slew  rate  of  11  deg/s,  and  a  dwell  time 
of  5  minutes.  During  the  dwell  period  approximately  15  sec¬ 
onds  is  required  for  transients  to  settle;  the  remaining  4-3/4 
minutes  are  used  for  sampling.  The  effective  carousel  rate  of 
this  motion  is  approximately  1.5  cycles  per  hour.  Further 
Bell  tests  indicate  that,  during  the  dwell  interval,  transient 
effects  seem  to  pass  after  about  15  seconds.  If  another  15 
seconds  were  allowed  for  sampling,  an  equivalent  carousel  fre¬ 
quency  of  0.25  cycles  per  minute  is  obtained.  Because  of  these 
promising  indications  it  is  likely  that  the  first  implementa¬ 
tions  of  carousel  motion  for  gradiometers  will  be  in  the  step 
indexing  mode. 
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